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MAGNITUD MAXIMA EN ZONAS PUENTES PARA ESTUDIOS DE 
PELIGROSIDAD SISMICA DEL ECUADOR 



Roberto Aguiar Falconf y Celia Castro Torres 



Centro de Investigaciones Cientificas 
Escuela Politecnica del Ejercito 
raquiar(S)espe.edu.ec 



RESUMEN 

Se presentan las magnitudes maximas esperadas en cada una de las veinte zonas fuentes que 
se han considerado para posterlores estudlos de pellgrosidad sismica del Ecuador. Las magnitudes 
maximas se han hallado tenlendo en cuenta la maxima magnitud registrada Instrumentalmente entre 
1 540 y 2009 y en base a una relaclon entre la longitud de ruptura y la magnitud. 

Prevlamente, se obtuvo una relaclon entre la longitud de ruptura y la magnitud, trabajando con 
las curvas de Intensidades sismicas en la escala de Mercalll de 19 sismos que causaron demaslado 
daho en las estructuras, en el Ecuador. Esta relaclon se obtuvo con la hipotesis de que la longitud de la 
superficle de replicas de un sismo, es la longitud de ruptura y es aproximadamente Igual a la maxima 
longitud de la curva de la Isosista de Intensidad VI 11. Se comprobo la valldez de esta hipotesis con los 
sismos de Baeza 1987, Macas 1995 y Bahia de Caraquez 1998. 



ABSTRACT 

The maximum expected magnitudes In each of the twenty considered source areas are 
presented, for further study of seismic hazard of Ecuador. The maximum magnitudes were found taking 
Into account the maximum Instrumentally recorded magnitude between 1540 and 2009, based on a 
relationship between the rupture length and magnitude. 

Previously, a relationship between the rupture length and magnitude was found, working with 
seismic Intensity curves In 19 seisms Mercalll that have caused heavy damage to the structures. In 
Ecuador. This relationship was obtained with the assumption that the length of the surface replica of an 
earthquake. Is the rupture length and Is approximately equal to the maximum length of the Isoselsmal 
curve of Intensity VIII. It was proved the validity of this hypothesis with the earthquakes of Baeza in 1987, 
Macas in 1995 and Bahfa de Caraquez in 1998. 
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1. INTRODUCCION 



En base a la tectonica, geologia y sismicidad registrada, se determinan las zonas fuentes, que 
son la base de los estudlos de pellgrosidad sismlca. Las zonas fuentes son volumenes de lltosfera en los 
cuales se producen slsmos de las mlsmas fallas geologicas o de la mlsma zona de subducclon. Para el 
Ecuador, Aguiar y Castro (2009) ban definido 1 1 zonas fuentes asoclados con el fallamlento local y 9 
zonas fuentes asocladas a la subducclon, las mlsmas que se presentan en las figuras 1 y 2. 

En la figura 3 se presentan cuatro cortes en las latitudes 1° S, 2.5° S, 4°S y 6°S., que ayudan a 
comprender la forma de las zonas fuentes en el piano vertical. En ellas tamblen se determine el espesor 
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Figura 1 Zonas fuentes asocladas al fallamlento local. 
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Ahora bien, en cada zona fuente se debe encontrar la magnitud maxinna que en ella se espera y 
la tasa anual de ocurrencia de slsmos mayores o Iguales a una deternninada magnitud. En este articulo 
solo se obtlene la magnitud maxima esperada en cada zona sismica. 



2. MAGNITUD MAXIMA 



Existen diferentes formulas empirlcas que relaclonan la longitud de rotura de la falla L, con la 
magnitud maxima Mmax- Una de las primeras relaclones fue propuesta por Idda (1959) para fallas 
inversas 

^»..=5.47 + 1.01ogL (1) 

M_=6.04 + 0.791ogL (2) 

La ecuaclon ( 1 ) es para slsmos profundos y la ecuaclon ( 2 ) para slsmos superficlales e 
Intermedlos. Por otra parte, Grases (1985), en base al estudio estadistico de 58 slsmos ocurridos a nivel 
mondial y lo deducldo con base en la geologia para fallas transcurrentes se tiene: 

4^ =6.04 + 0.704 log L (3) 

Donde la longitud de ruptura debe ser mayor a 35 km., y menor a 420 km. La maxima magnitud 
de los slsmos empleados en el estudio fue Ms = 8. 
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En la evaluacion de la pellgrosidad sismlca de Colombia para la normativa NSR-98, utlllzaron las 
ecuaclones propuestas por Ambrasseys para determlnar la magnitud maxima, estas son: 



L = g*' Limite Inferior 

L = Mejor Ajuste ( 4 ) 

j (1.615M„ -8.58) 7- • V C 

L = e " Limite Superior 
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Se aprecia que la relacion entre la Magnitud Maxima M^ax y el logarltmo de la longitud de rotura 
es de tlpo lineal, de la forma: 



=^log^ + fi 

Donde A y B son constantes que se obtlenen por regreslon lineal. Es convenlente que los 
valores de A y 6 se obtengan con sismos de una determlnada zona geografica como lo propone Acharya 
(1979). Por otra parte, Slemmons (1977) propone calcular la magnitud maxima en base al tlpo de falla, 
es asi como propone tres ecuaclones una para falla normal, otra para falla Inversa y otra para falla 
transcurrente, todo esto fue obtenido con sismos registrados en el Occldente de los Estados Unidos. 

Cuando en una zona fuente, no se tiene un tlpo de falla definido es convenlente utilizer la 
metodologfa propuesta por Acharya. Para America del Sur, Acharya propone la ecuaclon ( 5 ) para 
encontrar la magnitud maxima, la misma que fue obtenida a partir de 31 eventos con magnitud superior 
a 7. Por lo tanto, para magnitudes Inferlores su apllcaclon es Inclerta. 

M_=2.30 + 2.831ogL (5) 

Por otra parte, Dorbath et al (1990) proponen la sigulente ecuaclon para encontrar la magnitud 
M ^ a partir de la longitud de ruptura L . 

=1.62 log L+ 4.44 (6) 

Lo Interesante del trabajo de Dorbath et a! (1990) es la forma como determlnan la longitud de 
ruptura L. Lo hace midlendo el eje mayor de la curve de Isosista de Intensidad Vlll en la escala de 
Mercalll, para esto medlante el estudio de las replicas de varies sismos demuestra que la longitud mayor 
del area de las replicas es aproximadamente Igual a la longitud del eje mayor de la curve de Isosistas de 
Intensidad Vlll. 

Gallo y Tavera (2002) Incorporan un factor de ajuste ^ a la ecuaclon de Dorbath et a! (1990) 
para poder encontrar la magnitud a partir de Isosistas de diferentes Intensidades. 



= 1.62 log 



f-1 



+ 4.44 



(7) 



Cuando se trabaja con la Isosista de Intensidad Vlll el valor de K = \, para Intensidad Vll, se 
tiene K = 2. Para Intensidad VI, el valor de ^ = 4. Lo Importante es notar que se puede hallar la 
longitud de ruptura a partir de las curvas Isosistas de Intensidades. SI se trabaja con la curva de 
Intensidad Vlll se halla directamente la longitud de ruptura, si se trabaja con la curva de Intensidad Vll se 
divide para dos la longitud mayor de la curva de intensidad y si se trabaja con la curva de intensidad VI 
se divide para cuatro, con lo que se halla la longitud de ruptura, en forma aproximada. 
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3. ANALISIS DE ALGUNOS SISMOS REGISTRADOS EN ECUADOR 

Con el proposito de comprobar la propuesta de Dorbath et al (1990) y de Gallo y Tavera (2002) 
en el sentido de que la longitud mayor de la curva de Isosistas de Intensidad VIII es aproximadamente 
Igual a la longitud de la zona donde se producen las replicas y esta es a su vez es la longitud de ruptura 
del sismo. Se estudlan los sismos del 4 de agosto de 1998, conocldo con el nombre del sismo de Bahfa 
de Caraquez, Aguiar ef a/(1998). El sismo del 2 de octubre de 1995, conocldo como el sismo de Macas, 
Aguiar (2000) y el sismo del 3 de Marzo de 1 987, conocldo como el sismo de Baeza. Se destaca que las 
Isosistas de estos sismos fueron realizadas por Egred (1988, 1999). 




Figura 4 Replicas del sismo de Bahia de Caraquez y curva de Intensidad VII. 



En la figura 4 se presenta el epicentre del sismo del 4 de agosto de 1998 y las replicas 
registradas hasta diciembre de ese mismo aho. La longitud de la superficie de replicas es 75.48 km., 
Tambien se ha dibujado la curva de Intensidad VII, en ella se tiene que la longitud mayor es 150.96 km. 
Como es la curva de intensidad VII, para hallar la longitud de ruptura se divide para 2, Callo y Tavera 
(2002) con lo que se obtiene 96.65 km. 

En la figura 5 se aprecia el epicentre y las replicas del sismo del 2 de octubre de 1995 que tuvo 
una magnitud M ^ = 6.5 . La longitud de la superficie de las replicas es 65.72 km., y la longitud mayor 
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de la curva de intensidad VIII es 62.38 km. Se tiene una muy buena aproximaclon entre estos dos 
parametros. 

En la figura 6 se aprecia el epicentre del sismo de Baeza y las pocas replicas del mismo, 
lamentablemente solo se tienen 7 replicas registradas no porque no haya habldo sismos despues del 
terremoto del 5 de marzo de 1987, sino por que simplemente no se tienen registros. Con esta acotaclon 
la longitud de las pocas replicas es 45.38 km. La longitud de la Isosita de Intensidad VIII es 73.03 km. 




- 5 ' 1 1 1 ^ 

79 78 77 -7G 75 



Figura 5 Replicas del sismo de Macas y curva de Intensidad VIII. 



Si bien es cierto, son solamente tres sismos los que se ban analizado, sin embargo de ello hay 
una buena aproximaclon entre la longitud mayor de las replicas y la longitud inferida de la curva de 
isosista. Resultados similares los encontro Gallo y Tavera (2002) para el Peru. 



Tabla 1 Longitud de superficie de replicas e inferida de curva de isosista. 



Sismo 


Longitud de superficie 
de replicas (km) 


Longitud inferida de 
curva de isosista (km) 


Bahia de Caraquez (1 998) 


75.48 


96.65 


Macas (1995) 


65.72 


62.38 


Baeza (1987) 


45.38 


73.03 
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4. RELACION ENTRE LONGITUD DE RUPTURA Y MAGNITUD 

Una vez que se comprobo que las propuestas de Dorbath et al (1990) y de Gallo y Tavera 
(2002), dan resultados satisfactorlos para los sismos registrados en el Ecuador, se procedio a encontrar 
una relaclon entre magnltud maxima M ^ y la longitud de ruptura L. Pero antes de eso, los sismos de 

magnitud y M ^ se pasaron a magnltud M^, con las sigulentes ecuaclones. Cahuarl (2006) 

=0.9588 +0.458 

=0.7044 +1.702 



La primera ecuaclon es para sismos con magnltud entre 5.1 y 6.8. La segunda ecuaclon en 
camblo se aplica para sismos con magnltud M ^ comprendidos entre 4.5 y 6.8. 



2 

1 

0 
-1 

-2 
-3 
4 

-81 -80 -79 -78 -77 -76 -75 

Figura 6 Replicas del sismo de Baeza y curva de Intensidad VIII. 




La mayor parte de las isosistas de los sismos que se consideran en el estudio fueron realizadas 
por Egred (1988, 1999, 2009) pero tambien hay isosistas que constan en el estudio de peligrosidad 
sismica desarrollado para el Proyecto Hidroelectrico Coca-Codo Sinclair. Pasquare G., (1988). 
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En la figura 7 se presentan dos de estas isosistas y corresponden al sismo de 1797 que destruyo 
la antigua cludad de Rlobamba y a la del sismo de 1995 que llevo al colapso al puente colgante sobre el 
no Upano. En Castro (2009) se encuentran todas las Isosistas utlllzadas en el estudlo. 

En la tabla 2 se muestran los 19 sismos con los cuales se encontro una relaclon entre la longitud 
de ruptura L con la magnitud M^. Se Indica tamblen la longitud maxima de la curva de Intensidades, 

utlllzada y en la ultima columna de la tabla se muestra la longitud de ruptura utlllzando la propuesta de 
Gallo y Tavera (2002). 

En la figura 8 se presenta la magnitud y longitud de ruptura de los 19 sismos de la tabla 2, se 
Indica tamblen la recta hallada en este estudlo y la recta que se obtlene con la propuesta de Dorbath et 
a/ (1990). Se nota que la propuesta de Dorbath et a! (1990) reporta valores mas altos de magnitud pero 
son comparables sus valores con la hallada en este estudlo. La ecuaclon encontrada en este trabajo es: 

= 4 . 519 + 1.358 logL (9) 




5. MAGNITUD MAXIMA ESPERADA EN CADA ZONA FUENTE 

Egred (2009), en base a las Isosistas de 129 sismos, obtuvo un mapa de pellgrosidad sismica del 
Ecuador pero de Intensidades Sismicas en la escala de Mercalll. Este mapa se presenta en la figura 9. 

Pero en la figura 9, sobre el mapa de Intensidades maximas se han colocado las zonas fuentes 
asocladas a fallamlento local, descritas en el apartado 1, con el proposito de encontrar la magnitud 
maxima esperada en cada zona fuente, para ello se obtuvo la longitud de ruptura per medio de la curva 
de Intensidades y luego se encontro la magnitud maxima esperada en la fuente empleando las 
ecuaclones ( 9 ) y ( 6 ). Los resultados obtenidos se muestran en la tabla 3. En la ultima columna de esta 
tabla se ha colocado la magnitud maxima que ha sido registrada Instrumentalmente desde 1540 a 2009. 
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Tabla 2 Sismos considerados para relacionar longitud de ruptura con magnitud maxima. 



# 


SISMOS 


FECHA 


M 


Mw 


LONGITUD MAXIMA EN LAS 
INTENSIDADES 


LONGITUD DE 
ROTURA 


VIII 


VII 


VI 


L (km) 


1 


RIOBAMBA 


04-02-1797 


7.7 


7.70 


158.67 






158.67 


2 


IBARRA 


16-08-1868 


7.17 


7.17 


102.65 






102.65 


3 


ESMERALDAS 


31/01/1 906 


8.6 


8.60 


228.05 






228.05 


4 


PASTO CALLE 


1 5/09/1 944 


5.7 


5.92 






23.13 


5.78 


5 


PELILEO 


05/08/1 948 


6.75 


6.93 


52.36 






52.36 


6 


AMBATO 


05/08/1 949 


6.17 


6.37 


41.00 






41.00 


7 


NAPO 


11/05/1 955 


6.8 


6.98 


15.00 






15.00 


8 


BAHIA 


16/01/1 956 


7.3 


7.30 


78.82 






78.82 


9 


ESMERALDAS 


10/01/1 958 


7.8 


7.80 




105.71 




52.86 


10 


PASA 


30/07/1 960 


5.75 


5.97 




19.56 




9.78 


11 


ALAUSi 


08/04/1 961 


5.79 


6.01 




31.30 




15.65 


12 


BAHIA 


1 9/05/1 964 


5.4 


5.64 


76.19 






76.19 


13 


NABON 


22/03/1 966 


6 


6.21 


33.45 






33.45 


14 


PASTO CALLE 


06/10/1 976 


5.7 


5.92 


30.41 






30.41 


15 


ESMERALDAS 


09/04/1 978 


6.1 


6.31 




85.38 




42.69 


16 


BAEZA 


03/05/1 987 


7 


7.00 


73.03 






73.03 


17 


BAHIA DE 
CARAQUES 


08/04/1 998 


7.1 


7.10 




150.96 




75.48 


18 


MACAS 


03/10/1 995 


6.9 


6.90 


62.38 






62.38 


19 


PUJILI 


28/03/1 996 


5.7 


5.72 






39.39 


9.85 




Figura 8 Relacion entre magnitud y longitud de ruptura L . 
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En la tabla 3 se aprecia que las magnitudes maximas obtenidas, con la ecuaclon hallada en este 
estudio son mas parecldas a la magnitud maxima registrada Instrumentalmente, son mas parecldas en 
relaclon a las que se hallan con la ecuaclon de Dorbath. 



Tabla 3 Magnitudes maximas esperadas en cada zona fuente. 



ZONA 


mu 


L 

(km) 


{Este estudio) 


(Dorbath) 


(Registrada) 


ZONA 1 


VII 


209.24 


7.7 


8.2 


7.8 


ZONA 2 


VII 


61.05 


6.9 


7.3 


6.5 


ZONA 3 


VII 


80.12 


7.1 


7.5 


5.8 


ZONA 4 


VII 


107.47 


7.3 


7.7 


lA 


ZONA 5 


VIII 


186.96 


7.6 


8.1 


7.8 


ZONA 6 


VII 


99.31 


7.2 


7.7 


8.1 


ZONA 7 


VI 


67.87 


7.0 


7.4 


6.9 


ZONAS 


VI 


20.95 


6.3 


6.6 


6.8 


ZONA 9 


VI 


76.65 


7.1 


7.5 


7.3 


ZONA 10 


VI 


65.67 


7.0 


7.4 


7.0 


ZONA 1 1 


VI 


48.42 


6.8 


7.2 


6.4 



Las magnitudes maximas que deben considerarse en un estudio de peligrosidad sismica del 
Ecuador, deben ser mayores a las registradas instrumentalmente. Este criterio sirve para el case en que 
la magnitud maxima hallada en base a la longitud de ruptura sea menor que la magnitud maxima 
registrada instrumentalmente. Para los sismos de subduccion se trabajara con las magnitudes maximas 
registradas instrumentalmente. En base a estos resultados en la tabla 4 se presentan las magnitudes 
maximas que deben considerarse en cada zona fuente. 



Tabla 4 Magnitudes maximas esperadas en cada una de las zonas fuentes. 



Fallamiento Local 


Subduccion 


Zona 


M 

W 


Zona 


M 

W 


1 


7.8 


12 


8.6 


2 


6.9 


13 


7.5 


3 


7.1 


14 


6.9 


4 


7.4 


15 


7.1 


5 


7.8 


16 


6.4 


6 


8.1 


17 


6.4 


7 


7.0 


18 


7.1 


8 


6.8 


19 


7.5 


9 


7.3 


20 


6.6 


10 


7.0 






11 


6.8 
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MAPA DE INTENSIDADES MAXIMAS 




Escuela Pdiltenica Nactond - Instituto Geofisico ■ 



EKita - Intwwlsdv MSK 



Figura 9 Mapa de intensidades maximas registradas en el Ecuador. 



6. CONCLUSIONES 



En base al estudio de las isosistas de 19 sismos registrados en el Ecuador se ha determinado 
una relacion entre la Magnitud y la Longitud de ruptura, considerando que la longitud de la curva de la 
isosista de intensidad VIII es aproximadamente igual a la longitud maxima de las replicas del sismo y por 
ende de la longitud de ruptura. Previamente se vio la validez de esta hipotesis al encontrar estas 
longitudes en los sismos de Baeza de 1987, de Macas de 1995 y de Bahia de Caraquez de 1998. 

Posteriormente se aplico esta ecuacion, para determinar la magnitud maxima esperada en las 
zonas fuentes, trabajando con un mapa de peligrosidad sismica del Ecuador de Intensidades en la 
escala de Mercalli. La longitud maxima encontrada de esta manera, en cada zona fuente es muy 
parecida a la magnitud maxima registrada instrumentalmente entre 1540 y 2009. 

Finalmente para cada una de las zonas fuentes de fallamiento local se determine la magnitud 
maxima como el valor mayor entre la magnitud maxima registrada instrumentalmente y la magnitud 
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inferida a partir de la longitud de ruptura. Para las zonas fuentes asocladas a la subducclon, la magnltud 
maxima esperada, se considero Igual a la magnltud maxima registrada Instrumentalmente en el periodo 
ya senalado. 
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RESUMEN 

Se presentan 20 curvas de isosistas correspondientes a 1 9 sismos ocurridos en el Ecuador y 1 
slsmo ocurrido en Colombia, los mismos que causaron danos severos en las estructuras de nuestro pais. 

En base a un modelo de ley de atenuaclon se precede a elaborar en base a las curvas Isosistas 
una ley de atenuaclon acorde a las exigenclas geotectonicas de nuestro pais la cual sera vallda para las 
zonas de fallamlento local y para roca, la misma que sera evaluada en base a otras leyes elaboradas en 
el resto de paises de Suramerica. 



ABSTRACT 

The 20 Isoselsmal curves for 19 earthquakes occurred In Ecuador and 1 earthquake occurred In 
Colombia, the same that caused severe damage In the structures of our country. 

Based on a model law of attenuation Is to develop on the basis of a law Isoselsmal curves 
attenuation geotectonic commensurate with the demands of our country which will be valid for the zones 
of faulting and local rock, which will be evaluated on the basis of other laws made In other countries of 
South America. 

1. INTRODUCCION 

Las leyes de atenuaclon son expreslones semlempiricas que relaclonan magnltud-dlstancla- 
Intensldad sfsmica. Entendlendose por esta ultima palabra a la aceleraclon, velocldad, desplazamlento e 
Intensidad proplamente dicha de eventos sismicos; estas relaclones se obtlenen de los dates que existen 
sobre los parametros menclonados. Aguiar (1982). 

Una hipotesis, con la cual se obtlenen estas leyes de atenuaclon, es suponer que los 
movimlentos del suelo se van atenuando conforme la distancia epicentral se Incrementa. En termlnos 
generales esto es clerto, aunque hay cases, en que a una misma distancia epicentral se tengan varlas 
aceleraclones del suelo o que a mayor distancia epicentral se tenga una mayor aceleraclon, son cases 
que se dan pero por lo general a mayor distancia menores Intensidades de los movimlentos del suelo. 
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Figura 1 Aceleraciones maximas registradas en el sismo de 1 979 en Imperial Valley. 

En la figura 1, se presentan las aceleraciones maxinnas del sismo de 1979 en Imperial Valley, 
California, registrado en diferentes estaciones. A 1 km., del epicentre, se aprecia en la figura 1, que la 
aceleracion maxima de la componente horizontal vario entreO.4 y 0.8 g., siendo g., la aceleracion de la 
gravedad; a 10 Km., del epicentre estuvo entre 0.15 y 0.60 g., a 100 km., del epicentre vario entre 0.01 y 
0.04 g. Miranda (2002). 

En la figura 1, se observe la curve de valores medios, que vendria a representar la ley de 
atenuacion para el sismo de Imperiai Valiey de 1979 y una muy buena cantidad de registros se 
encuentran entre el valor medio ± a . Siendo a ia desviacion estandar. Per io tanto, cuando se hablen 
de leyes de atenuacion, por la incertidumbre que ellas tienen es necesario incluir un termino que 
contenga la desviacion estandar. 

Se puede pensar que obtener una ley de atenuacion de movimiento de suelo, es bastante 
sencilio que todo Io que se requiere es tener suficiente cantidad de registros sismicos y apiicar tecnicas 
de regresion iineai para encontrar la ley, pero no es asi. En efecto, Ordaz et al (1989) a partir de 64 
registros de ia red sismica de Guerrero con magnitudes entre 3 y 8.1, obtuvieron ia siguiente ley de 
atenuacion siguiendo ia metodoiogia propuesta por Joyner y Boore (1981). 
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Donde es la aceleracion maxinna del suelo en cm/s^; E{ ) es el valor esperado; R es la 
minlnna distancia al area de la falla; M es la magnltud. La desviaclon estandar del logarltmo comun de 
la es 0.25. Ahora blen, Ordaz (1999) compara las aceleraclones maxlnnas registradas (cuadros) 

durante el slsmo del 21 de octubre de 1 995 que tuvo una magnltud de 7.2., a la Izquierda de la figura 2 y 
la obtenida con la ecuaclon (1 ) y se aprecia que todos los valores registrados tienen aceleraclones 
mayores a las predichas; se ha dibujado la curva de valores medlos y las curvas correspondlentes a los 
percentiles 16 y 84, que corresponden a menos y mas una desviaclon estandar. 




Figura 2 Aceleraclones maximas registradas y estimadas con ecuaclon 1. Ordaz (1999). 



A la derecha de la figura 2 se muestra las aceleraclones maximas registradas (cuadros) durante 
el slsmo del 25 de febrero de 1996, de magnltud 7.1 y las predichas con la ecuaclon (1). Ahora los 
valores registrados son menores a los predichos con la ecuaclon (1) para valores de i? < 400 km. y son 
mayores a las que se hallan con la ecuaclon (1 ) mas a . 

Todo esto demuestra que hay mucha Incertidumbre en la formulaclon de leyes de atenuaclon, de 
pronto hay otras variables que deben ser Incorporadas en la ley. Pero per otro lado tamblen se debe 
pensar en utilizer valores mas altos de la desviaclon estandar. 

2. CARENCIA DE REGISTROS 

En un buena parte de Raises, con alta pellgrosidad sismica no se tienen acelerogramas de 
sismos fuertes, con lo que se complica aun mas la obtenclon de leyes de atenuaclon. Entonces no queda 
otra opcion, que determiner leyes de atenuaclon de Intensidades en la escala de Mercalll a 

sablendas que esta Intensidad tienen tamblen Incertidumbre, como per ejempio la escala no considera el 
tlpo de suelo en la evaluaclon del grade de daho de una estructura y es muy blen conocldo que el tlpo de 
suelo si Influye en el comportamlento de la estructura. 

En forma general las curvas Isosistas si dan una Idea de la atenuaclon del movimlento del suelo 
durante un slsmo, si IndIcan per donde se atenuo mas rapidamente el movimlento y per donde se 
propago mas lentamente este movimlento, de tal manera que si es bastante conflable encontrar leyes de 
atenuaclon de Intensidades a partir de las curvas de Isosistas. 

Lo dificll es como pasar esas Intensidades a aceleraclones maximas de movimlento del suelo. La 
forma de hacerlo es medlante una ecuaclon de la sigulente forma: 
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Jmm+7 



( 2 , 



Algunos valores de o y y se indican en la tabla 1 y han side obtenidos en base a dates de 
Intensidades, y aceleraciones maximas del suelo, 



Tabla 1 Valores de a y y 



AUTOR 


a 


Y 


Trifunac y Brady (1975) 


0.300 


-0.014 


Bolt 


0.313 


0.340 


Murphy, O' Brien (1977) 


0.250 


0.250 


Lomnitz (1974) 


0.333 


-0.500 


Saragoni ef a/(1982) 


0.345 


-0.222 


Saegesser y Baumgartner 


0.26 


0.19 


(1978) 


0.31 


-0.4 


Medvedev-Sponheuer (1969) 


0.30 


-0.1 



En las cinco primeras ecuaciones la Intensidad corresponde a (Mercalli Modificada), en las 
dos ultimas es (Medvedev, Sponheuer y Karnik). Estrictamente la ecuacion de Medvedev- 

Sponheuer (1969) es lasiguiente: 

0.3 - 0.4 < log A_ < 0.30 7^,^ - 0. 1 ( 3 ) 



Esta ultima ecuacion fue utilizada en la evaluacion de la peligrosidad sismica del proyecto 
hidroelectrico Coca-Codo-Sinclair, que esta en construccion actualmente en el Ecuador. Pasquare 
(1988). 

3. LEYES DE ATENUACION MACROSl'SMICAS 



En la literatura existen una gran cantidad de modelos de leyes de atenuacion para Intensidades, 
entre los que se destacan los propuestos por Sponheuer (I960); Shebalin (1968); Howel y Schultz 
(1975), entre otros. 



El modelo de Sponheuer (1960), basado en el trabajo de Kovesligethy (1907) fue utilizado en el 
proyecto Coca-Codo-Sinclair y es el siguiente: 



^ 7 ) A 



h 



(4) 



y n y 



Donde 7^ es la intensidad macrosismica en el epicentre; 7^ es la intensidad de la isosista con 



radio medio R . ; D. es la distancia hipocentral, D. = ^Rf + ; /z es la profundidad hipocentral; 

k, bson coeficientes que deben determinarse en base a las isosistas macrosismicas. 
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El modelo de Shebalin( 1968) fue derivado a partir de la atenuaclon de la energia disipada 
desde un foco puntual. La ecuaclon propuesta es la sigulente. 

A + B\nR + C R ( 5 ) 

Todas las variables ban side definidas, faltando IndIcar queA, S, C, son coeficlentes a 
deternninar en base las curvas Isosistas. 

Otra ley basada en conceptos energeticos es la de Blake-Sheballn y es adecuada para 
atenuaclones fisicas pequenas. Su ecuaclon es la sigulente. Sheballn (1968) 



=-^o~rlog (6) 

En este caso el coeficlente y se deternninara con las curvas Isosistas disponibles. 

En los modelos presentados no se Incluye la Magnitud, como parte de la ley de atenuaclon, 
slendo este un parametro fundamental, de ahi la Importancia del modelo propuesto por Howel y Schultz 
(1 975) que se Indica a continuaclon. 

L= AM +B\aR+C R + D^-^ (7) 

Donde M es la magnitud y A, B, C, D son los coeficlentes a determlnar con las Isosistas 
disponibles. Existen algunas varlantes a la ecuaclon ( 7 ) como se observe en la table 2 en que se 
presentan algunas leyes de atenuaclon de Intensidades para pafses de America del Sur. 

Las Intensidades se transforman en aceleraclones maximas del suelo, medlante la ecuaclon ( 2 ) 
y utlllzando alguna de los coeficlentes Indicados en la table 1. Para Venezuela, por ejempio, se trabajo 
con la propuesta de Murphy y O’ Brien (1977). Para el Ecuador con la propuesta de SaragonI et al 
(1982). Con esta acotaclon en la table 3 se presentan las leyes de atenuaclon de aceleraclones 
obtenidas a partir de Isosistas. 

En la figure 3 se presentan las curvas de atenuaclon de la aceleraclon maxima del suelo, 
encontradas con las ecuaclones de la table 3. Para Venezuela se trabajo con la ley propuesta por 
Quijada et al (1993), que por clerto se denomlna ley de INTEVEP. Se obtuvo la atenuaclon para un 
sismo de magnitud M ^ = 6, para el caso de Venezuela este sismo corresponde a M ^ = 5.92. Es 
necesarla esta aclaraclon debldo a que en la propuesta de Quijada et al (1 993) se trabaja con magnitud 
M ^ y en las restantes ecuaclones con M Tamblen se debe aclarar que Quijada et al (1993) trabaja 

con la distancia al hIpocentroD, en camblo en las otras ecuaclones se trabaja con la distancia al 
epicentre /? . En la obtenclon de las curvas que se encuentran en la figura 3 se considero que D = R . 

La ecuaclon propuesta por SaragonI et al (1982) para Chile y Argentina presenta valores muy 
bajos de aceleraclon maxima del suelo. La ecuaclon de Casaverde (1980) para el Peru es la que reporta 
valores mas altos y las tres restantes tienen valores comprendidos entre estas dos ecuaclones. Se 
considero que la desviaclon estandar es nula en las curvas de la figura 3, es declr se trata de curvas de 
valores medlos. 
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Tabla 2 Leyes de atenuacion de Intensidades macrosismicas. 



Autor 


Pais de aplicacion 


Ley de atenuacion 


G rases (1980) 


Venezuela 


/ = /(,+ 2.60 - 0.004 li? - 2.221og(i? + 20) 
Oriente 

/ = /g + 1 .20 - 0.0063i? + 1.68 \og(R + 20) 
Centro y Occidente 


Gershanik y Gajardo 
(1983) 


Venezuela 


I = 6.44 + 0.92 M, - 3.331og(D + lO) 
f^iogA/, =0.70 


Gomez y Salcedo 


Colombia 


I = Io~ 0.019(i? - 15)- 0.76 ln(i? - 15) 
a = 0.29 Para h< 20 km 
I = (2.3M, -6.8)*l.li?^°“ 

<7 = 0.51 Para h< 60 km 

/ = /g *1.26*e“°“"^ 

r = 0.89 (correlacion) Para h> 60 km 


Aguiar (1982) 


Ecuador 


/ = /g - 1.3 - 0.0045 i? - 0.61 log R 
a = 0.6 Para la Sierra 
/ = /q - 0.7 - 0.0087 R + 0.26 log R 
a = 0.6 Para la Costa 


Aguiar (1988) 


Ecuador 


I = 8.645 + 1 .256 M, - 2.226 ln(i? + 40) 
<7 = 0.5 


Saragoni (1980) 


Peru y Chile 


/ = 8.19 + 0.94 M - 1.761n(i? + 40) 



Tabla 3 Leyes de atenuacion de aceleraciones maxinnas del suelo 



Autor 


Pais de aplicacion 


Ley de atenuacion 


Gershanik y Gajardo 
(1983) 


Venezuela 


log = 4.28 + 0.53 M, - 0.83 ln(D + lO) 

(7 = 0.70 


Quijadaefa/ (1993) 


Venezuela 


In Tl^ax = 5.40 + 0.36 M ^ -0.86 ln(D + lO) 
(7 = 0.66 


Aguiar (1988) 


Ecuador 


In = 6.35 + 0.99 M, - 1 .76 \n{R + 40) 

(7 = 0.6 


Saragoni et al (1982) 


Chile-Argentina 


In = 8.54 + 0.57M, - 1.73 ln(i? + 60) 


Saragoni etal (1982) 


Peru 


In A„^, = 8.18 + 0.68M, - 1.63 ln(i? + 60) 


Casaverde (1980) 


Peru 


In A_=4.23 + 0.8M,-ln(i? + 25) 
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4. LEYES DE ATENUACION INSTRUMENTALES 

En el apartado 1 se presento una ley de atenuaclon encontrada per Ordaz et at (1989) para 
Mexico, en base a registros sismlcos. Varlas ecuaclones slmllares ban sido obtenidas por Me Guire 
(1978), ecuaclon 8; Esteva y Vlllaverde (1973), ecuaclon (9); que fueron utlllzadas en el Proyecto 
HIdroelectrIco Coca-Codo Sinclair. Estas ecuaclones son las sigulentes: 

In =3.60 + 0.89 M- 1.17 In D (8) 

<7 = 0.62 

In =8.63 + 0.8M-21n(D + 40) (9) 

<7 = 0.64 

Pero para evaluar la pellgrosidad sismlca de Ecuador. Aguiar (1982) utilizo la propuesta de 
Mc.GuIre (1976) sigulente: 

In A^^^ =6.16 + 0.64M-1.301n(i? + 25) (10) 

En la tabla 4, se muestran varlas leyes de atenuaclon encontradas en base a registros sismlcos y 
que toman en consideraclon el tlpo del suelo, el tlpo de falla, el coeficlente de atenuaclon Inelastico, entre 
otros factores. Son ecuaclones mas elaboradas en relaclon a las presentadas en las ecuaclones 8 a 10. 
Para Venezuela, Bendito (2001 ) y Bendito et at (2001 ) determlnaron en base a las ecuaclones de la tabla 
4, ecuaclones para roca, estas se presentan en la tabla 5, las ecuaclones de Sadigh et at (1993) y de 
Idris (1993), son las Indicadas en tabla 4 para una falla obllcua. En la figura 4 se presenta la varlaclon de 
la aceleraclon maxima del suelo en roca, para una magnitud de 6, en la que se ha Incluldo la ecuaclon 
propuesta por Quijada et at (1 993). 

Como se Indico, la ecuaclon de Quijada et at (1993) fue obtenida en base a curvas de Isosistas 
de macrosismos. En camblo las leyes de Boore et at (1993), Campbell (1993), Sadigh ef a/(1993) e Idris 
(1993) fueron halladas en base a registros de acelerogramas, sin embargo de ello se aprecia en la figura 
4 que la ecuaclon de Quijada et at (1993) presenta valores medlos a los hallados con las otras 
ecuaclones para distanclas epicentrales mayores a 20 Km., para distanclas menores a esta cantidad la 
aceleraclon encontrada es baja, con relaclon a las reportadas por las otras ecuaclones. 
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Figura 3 Leyes de atenuacion para varies paises de Sur America. 



Tabla 4 Leyes de atenuacion obtenidas con registros sismicos. 



Autor 


Ley de atenuacion 


Boore etal (1993) 


\ogY = b^+b 2 {M ^-6) + b^{M„-6f +b^D + b^\ogD + bf^Gs +b-,Gc 
Y es la aceleracion maxima horizontal en unidades de g. 
b^ =-0.136 Para fallas transcurrentes; b^ =-0.051 Para fallas inversas 
= 0.229; &3 =0; b,= 0; b, = -0.778; b^ = 0.162; b, = 0.251; 
h = 5.57 ; Gg,G^ dependen del tipo de suelo. 


Campbell (1993) 


InF = +0.683M + J3i tanh[0.647(M -4.7)]-1.01nr-a/? + 

0.27F + [{j3^ - 0.105 In R)]s + (5^ tanh(0.620 ) 


r = + [0.0586 *exp(0.683M)f 

or = 

Y es la aceleracion maxima horizontal en unidades de g. 

M = M ^ para M <6. M = M ^ para M > 6 
F = 0 para fallas transcurrentes; F = l para fallas inversas. 

5=0 para suelos cuaternarios; 5 = 1 para depositos terciarios y para roca 
a es un coeficiente de atenuacion anelastico (km'^). 
es la profundidad del lecho rocoso. 
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Para roca: J3„ = -3.15; j3^=j3^=j3^= 0; J3^ = 0.0150; jB, = -0.000995 


Sadigh et al 

(1993) 


In F = C, + + C, \n[R + + C,M J] 

Y es la aceleraclon maxima horizontal en unidades de g. 

Para Falla Transcurrente y roca con M ^ <6.5 

In F = -0.624 + 1 .0 M „ - 2. 1 ln[i? + exp(l .29649 + 0.250M „ )] 

Para Falla Transcurrente y roca con M ^ >6.5 

In F = -1 .274 + 1 . IM „ - 2. 1 \n[R + exp(-0.4845 1 + 0.524M „ ] 

Para Fallas Inversas se multipllcan las Fallas Transcurrentes por 1 .2 
Para Fallas obllcuas se multipllcan las Fallas Transcurrentes por 1.09 


IdrIss (1993) 


In F = [oTg + exp(orj + a 2 M )] + 

- exp(;^i + )]ln(i? + 20) + 0.2 F 

F es la aceleraclon maxima horizontal en unidades de g. 

M = para M <6. M = M ^ para M > 6 

F = 0 para fallas transcurrentes; F = l para Inversas; F = 0.5 obllcuas 
Para fallas transcurrentes y roca con M < 6.5 

In F = [- 0. 1 5 + exp(2.26 1 - 0.083M )] - exp(l .602 - 0. 1 42M ) ln(i? + 20) 
Para fallas trasncurrentes y roca con M > 6 

InF = [- 0.050 + exp(3.477 - 0.284M)]- exp(2.475 - 0.286M)ln(i? + 20) 



En la figura 4 tambien se aprecia una muy buena correlacion entre las ecuaclones de Campbell 
(1993) y de Sadigh et al (1993). Los resultados hallados con la ecuaclon propuesta por Idriss (1993) 
tambien son semejantes. En todas las curvas de la figura 4 no se ha Incluldo la desviaclon estandar. 



Tabla 5 Leyes de atenuaclon, para roca utlllzadas en Venezuela por Bendito (2001) 



Autor 


Ley de atenuacion 


Boo re et al 
(1993) 


log F = -0. 105 + 0.229(M „ - 6) - 0.7781og Vi?" + 5.57 "+ 0. 1 62 
F es la aceleraclon maxima horizontal en unidades de g. 

5.0<M„<7.0 

R < 100 km 


Campbell 

(1993) 


In F = -3. 15 + 0.683M -l.Oln (r)- (0.0150-0.000995 M)i?- 0.105 Ini? 
F es la aceleraclon maxima horizontal en unidades de g. 


r = ^i?" + [0.0586 *exp(0.683M)f 
M = M ^ para M <6. M = M ^ para M > 6 
M>4J R< 300 km 
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Figura 4 Comparacion de leyes de atenuacion obtenidas instrumentalmente. 

Se destaca que todas las leyes de atenuacion, hasta ahora Indicadas, son apllcables a zonas 
fuentes de fallamlento local, no para sismos de subducclon. 

Para sismos de subducclon, Bendito (2001), utilize la ecuaclon de Crouse (1991)) que fue 
hallada en la zona de Cascadia en Oregon ya que segun este autor estas son parecldas a la subducclon 
de Colombia. La forma de la ecuaclon general de Crouce (1 991 ) es la sigulente: 

In^max = Pi+ Pi + Pi Ml + p, ln[R + p^exp{p^Mj+ h] ( 11 ) 

Donde M ^ es la magnitud de memento; R es la distancia epicentral en km.; h\a profundidad 
focal en km; es la aceleraclon maxima del suelo en gals. Para determlnar la aceleraclon horizontal 
maxima del suelo se tiene: 

In = 6.36 + 1 .76 - 2.73 ln(7? + 1 .58 exp(0.608 ))+ 0.009 16 h 

En la figura 5 se muestra la varlaclon de la aceleraclon horizontal del suelo, cuando se aplica la 
ecuaclon de Crouse (1991) para profundidades focales de 50, 100 y 150 km. Para un sismo de magnitud 
M ^ = 6 . Tamblen se presenta la ley de atenuacion de Campbell (1993) para fallamlento cortical, con el 

proposito de ver la necesidad de contar por lo menos con dos tipos de ley de atenuacion, una para ser 
apileada en las zonas fuentes asocladas con la subducclon y otra para las zonas fuentes de fallamlento 
local. 
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Figura 5 Leyes de atenuacion de Crouse para subduccion y de Campbell para sismos corticales. 

5. DETERMINACION DE LAS LEYES DE ATENUACION DE INTENSIDADES. 

En base a 20 curvas Isosista correspondlentes a 19 sismos ocurridos en el Ecuador y 1 en 
Colombia se determlna, previo a un modelo base de ley de atenuacion se precede a determlnar la ley 
que corresponda a un modelo de atenuacion de Intensidades acordes a las condiclones sfsmicas de 
nuestro pafs. 

I = Cj + (^2^ C + i?o) (12) 

Donde, 1: Intensidad esperada a una Intensidad hipocentral D; Ci, C2, C3: Se Inflere del ajuste 
con los dates zonales correspondlentes.; Rq: Parametro empfrico que define el campo cercado. La 
ecuaclon 12 sirve para Isosistas cerradas, en caso de las Isosistas ablertas se hara una relaclon de 
angulos. 



La gran mayorfa de las Isosistas utlllzadas en este estudio fueron tomadas de varies aytores que 
ban realizado mapas de Isosistas en el Ecuador uno de los cuales es (1988, 1999, 2009) quien recoplla 
la Informaclon existente acerca de los principales terremotos en el Ecuador, asi como tamblen se 
utlllzaron Isosistas que constan en el estudio de pellgrosidad sismica desarrollado para el Proyecto 
HIdroelectrIco Coca-Codo Sinclair. Pasquare G., (1988). A continuaclon en las sigulentes graficas 
presentamos las diferentes Isosistas utlllzadas para este estudio: 
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81*W 80*W 79*W 78*W 




Figura 12 Isosistas del Sismo de Pasa 
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En base a los valores de Ri y los hipocentros de cada uno de los sismos se deternninan los 
valores de D . tal come se muestra en la tabla 6 de este apartado. 



Tabla 6 Valores de D para cada Intensidad de isosistas 



# 


SISMOS 


FECHA 


Mw. 


D (Km). 


mm 


IX 


VIII 


VII 


VI 


V 


IV 


mm 


II 


1 


RIOBAMBA 


04-02-1797 


7.70 


25.25 


49.14 


88.69 


136.20 


182.40 


230.90 


276.76 


392.23 




2 


OTAVALO 


16-08-1868 


7.17 


44.79 


81.23 


105.31 


233.45 


271.51 










3 


ESMERALDAS 


31/01/1906 


8.60 




29.48 


84.76 


153.57 


215.94 


291.26 








4 


PASTO CALLE 


15/09/1944 


5.92 








34.03 


38.47 


45.15 








5 


PELILEO 


05/08/1948 


6.93 


11.46 


14.80 


20.95 


28.94 


37.50 


53.90 


77.24 






6 


AMBATO 


05/08/1949 


6.37 




77.79 


93.89 


143.76 




334.23 








7 


NAPO 


11/05/1955 


6.98 






40.06 


66.71 


180.72 










8 


BAHIA 


16/01/1956 


7.30 






32.60 


53.00 


107.26 


214.95 


326.78 






9 


ESMERALDAS 


10/01/1958 


7.80 








101.97 


146.89 


219.59 


284.28 






10 


PASA 


30/07/1960 


5.97 








22.14 


24.28 


30.99 


55.86 


88.02 




11 


ALAUSi 


08/04/1961 


6.35 








26.63 


43.15 


57.47 


83.58 






12 


BAHIA 


19/05/1964 


6.40 






44.39 


58.70 


85.40 


109.02 


176.52 






13 


NABON 


22/03/1966 


6.21 






44.03 


49.39 


76.51 


92.27 


174.11 






14 


PASTO CALLE 


06/10/1976 


5.92 






12.64 


17.82 


24.02 


43.76 


74.62 


103.96 




15 


ESMERALDAS 


09/04/1976 


6.31 








15.96 


105.61 


173.95 


229.08 


278.52 




16 


TUMACO 


12/12/1979 


7.70 




27.49 


68.94 


163.99 


250.60 










17 


BAEZA 


03/05/1987 


7.00 




16.61 


28.70 


56.78 


80.47 


106.85 


132.32 


200.64 




18 


BAHIA DE 

carAquez 


08/04/1998 


7.10 






350.53 


354.02 


369.14 


381.07 


422.51 


430.89 




19 


MAGAS 


03/10/1995 


6.90 






160.93 


202.61 


222.90 


247.40 


257.29 






B!l 


PUJILI 


28/03/1996 


5.72 








33.57 


42.07 


60.69 


83.82 


118.03 


174.42 



Se determlnan los valores que exige la ecuaclon modelo para obtener finalmente la formula de la 
ley de atenuaclon. En Castro (2009), se presenta este calculo de forma mas detallada. Como el objetivo 
es mlnimizar error en el calculo, planteamos la sigulente ecuaclon: 

e(c,C„C,)=f^ {I,-[C, + C,M + cMd + K)]Y ( 13 ) 

/=! 



En la cual n equivale al numero de puntos que se anallzan y e es el error. En base a esta 
ecuaclon se obtlenen los valores de los dates zonales correspondlentes, encontrados a partir de un 
sistema de ecuaclones llneales que se describen con mayor detenimlento en Castro (2009). Obtenlendo 
de esta forma los sigulentes resultados: 

Ci= 3.3577 
Cs= 1.0013 
Ca= -0.8856 



Reemplazando los resultados en la ecuaclon 12 que es nuestra ecuaclon modelo tenemos que: 
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/, = 3.3577 + 1.0013Mw-0.8856Ln(D + 10) (14) 

La ecuacion 14 representa la variacion de la Intensidad en funclon de la magnltud, para 
deternninar la variacion de la aceleraclon maxlnna del suelo en funclon de la magnltud se reemplazan los 
valores en la ecuacion 2 y se toman los valores de de a y de y utlllzados por Saragonl(1982) los cuales 
se detallan en la tabla 1, los mismos que han sido utlllzados en los estudlos de Aguiar (1989) y Palacio 
(1987). 



En base a esta nueva Informaclon se obtuvo la nueva ley de atenuaclon en base a dates 
sismicos de nuestro pais, tal como se define a continuaclon: 

A _ 2.197+0.807Mrv-0.714(ln(D+10)) 

^MAX ~ ^ (15) 

6. SELECCION DE LEYES DE ATENUACION. 

Una vez estudladas las diferentes leyes de atenuaclon para zonas de fallamlento local y zonas 
de subducclon se realize un anallsis comparative de las leyes estudladas en America del sur y la ley 
obtenida en este estudlo. 

En la figure 6 se detallan las curves de las leyes de atenuaclon en funclon de las aceleraclones 
maximas del suelo con una magnltud Mb=7. En esta figure se particularizan las curves obtenidas 
medlante ecuaclones obtenidas en base a Isosistas las mismas que son comparadas con el modelo 
alcanzado para este estudlo. 

En la figure 7, las curves que se detallan representan a las leyes obtenidas Instrumentalmente, a 
excepcion de la curve determlnada por Aguiar Castro (2009), puesto que este estudlo fue realizado por 
medio de Isosistas, sin embargo notamos que la tendencia de esta curve se asemeja a la propuesta por 
Boore et al (1993), la cual fue utlllzada en Venezuela por Bendito (2001) para las condiclones de suelo 
roca. 



Se determine que debldo a que las comparaclones de las curves de Aguiar, Castro (2009) y las 
determlnadas medlante estudlo de Isosistas en America del Sur cumplen con una tendencia y se podrfa 
declr que son simllares, en base a este anallsis la ley obtenida en este estudlo es apllcable para la 
determlnaclon de la Pellgrosidad Sismica. 

En base a la comparaclon de la Figure 23 se aprecia que las leyes obtenidas por Boore et al 
(1993) para roca de forma Instrumental y la ley obtenida en este proyecto tienen de Igual manera una 
tendencia similar. Hay que recalcar que la Ley de Boore et al (1 993) y la ley de Aguiar Castro (2009), son 
leyes proplas para un perfodo estructural Igual a 0 el mismo que corresponde a suelo roca. 

Para leyes de atenuaclon en el case de Subducclon se trabajara con la ecuacion propuesta por 
Crouce (1991), debldo a que se emplea para zonas de Subducclon en Colombia, las cuales tienen un 
simllltud con las zonas de Subducclon de Ecuador. 
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Figura 22 Leyes de Atenuacion para America del Sur. 




Figura 23 Comparacion de curvas obtenidas instrumentalmente.. 
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7. Conclusiones 

• La carencia de registros sismicos dificulta la obtencion de las leyes atenuaclon, lo que obllga la 
utlllzaclon de las Intensidades en la escala de Mercalll, mlsma que no es suflclentemente 
conflable. 

• Para algunos casos autores que deternninan leyes de atenuaclon en America del Sur sugleren 
varies modelos de ley de atenuaclon en funclon de aceleraclon del suelo debldo a que las 
condiclones sismicas son distintas en cada region. 

• Se elaboro una ley de atenuaclon para nuestro pais en base a 20 Isosistas que representan 
sismos que causaron gran dano en las estructuras del Ecuador, la mlsma que fue comparada con 
otras leyes elaboradas en America del Sur y como resultado se obtuvo que la ley de este estudio 
tenfa una similar tendencia a estas. 

• Se determine la ley de Crouce(1991) para la zona de subducclon de Ecuador debldo a que esta 
fue utlllzada en para la zona de subducclon de Colombia, la mlsma que es similar a la zona de 
subducclon de Ecuador. 
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RESUMEN 

El acero come material alternative en la construcclon de edificlos ha side considerado en los 
ultimos ahos en el Ecuador, la necesidad de conocer las ventajas economicas y estructurales respecto al 
hormigon armado son el principal objetivo de este trabajo que se presenta. 

El alto grade de sismicldad que presenta el Ecuador hace necesarlo que el diseho sismo resistente sea 
primordial en los edificlos de varies pIsos, a traves de un anallsis matriclal de estructuras se obtlene las 
respectivas fuerzas estaticas equivalentes a las del sismo rare del CEC2000. 

Se realize dos tipos de anallsis el primero en funclon de la derive de piso y el segundo en funclon de la 
demanda capacldad de los elementos estructurales de acero. 

El anallsis del sistema estructural de hormigon armado en funclon de la derive de pIso es el que se 
propuso en el trabajo que se presenta. 

La comparaclon estructural entre ambos sistemas estructurales es el resultado final del anallsis 
estructural de este trabajo. 



ABSTRACT 

Steel as an alternative material In the building constructions have been considered In the last 
years In Ecuador, the necessity of knowing the economic and structural advantages about concrete are 
the principal goal of this article. 

The high degree of seismic activity that In Ecuador occurs makes necessary that the seismic design Is 
primordial In buildings. A matrix structural analysis In order to get the static forces of CEC2000 rare 
seismic Is applicable In this work. 

In this article there are two types of analysis, the first Is about the drift and the second Is about the 
Demand Capacity of the structural elements of steel. 

The structural concrete system has a drift analysis which Is the propose of this work. 

The structural comparison between both systems Is the final result of this work which Is presented. 

1. INTRODUCCION 

La reclente gran acogida de los sistemas estructurales de acero en el Ecuador, ha hecho que 
se presente clertas Incognitas comparativas con su par de hormigon armado, el cual por muchos ahos ha 
side la unica alternativa constructiva. En los pafses Industrlallzados, come los Estados Unldos, la 
construcclon de edificlos de acero se viene dando desde el aho de 1920, con la construcclon de grandes 
edificlos metallcos, este debldo a la producclon de acero en este pafs, lo que ha dado come resultado 
que practicamente una gran mayorfa de los edificlos de este material sean preferldos. 

Los criterlos estructurales que se utlllzan para los edificlos de hormigon armado tamblen son utlllzados 
para los edificlos de acero, especlalmente conceptos de mecanica de materlales, Ingenleria estructural y 
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como una parte fundamental y muy importante la tecnologfa de soldadura. Las consideraclones de 
diseno sfsmico Incluyen la utlllzaclon del espectro Inelastico de dlsefio del CEC2000, el cual permite 
obtener las aceleraclones de diseno para posterlormente con un anallsis modal espaclal de la estructura 
poder hallar las fuerzas estaticas equivalentes del sismo. En los ultimos afios en el Ecuador se esta 
Innovando con la elecclon del acero como material alternative para la construcclon de edificlos, el 
mejorar las condiclones de construcclon, abaratar los costos, recuperar la Inversion de forma pronta, son 
las razones fundamentales por las cuales se esta prefirlendo al acero en la construcclon de edificlos de 
varies pisos. 

2. MARCO TEORICO 



Se considera un sistema estructural de echo pisos, el mismo que se disena en acero y en 
hormigon armado sin muros de corte. La primera planta consta de una altura de 4.0 m, y los sigulentes 
pisos de 3.5 m, se consideran cuatro vanos de 5 m cada uno en los dos pianos ortogonales verticales. 




Figura 1.- Esquema del sistema estructural de echo pisos 



2.1 Obtencion de las Fuerzas Si'smicas estaticas equivalentes del sismo 



Para la obtencion de las fuerzas sfsmicas estaticas equivalentes se realize un anallsis modal 
espaclal, obteniendose previamente la matriz de rigidez de la estructura, a traves del programa riaxinfi. 
(Aguiar, 2007, 2008), en el cual se introduce como dates las propiedades geometricas de columnas y 
vigas, ademas de su longitud, el modulo de elasticidad del material, el numero de pisos y el numero de 
nudos restringidos. 

Se calcula la matriz de rigidez en coordenadas de piso, a traves de la formula, (Aguiar, 2004): 







;=1 



Ci3 



Donde: Ke es la matriz de rigidez de la estructura en coordenadas de piso; A es la matriz de 
compatibilidad de deformaciones; KL es la matriz de rigidez lateral de portico; NP es el numero de pisos 
de la estructura. 
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A traves del programa modalespacialSGDLnew.m (Aguiar, 2008) se realiza un analisis modal espacial 
con 3 grades de libertad por planta, y considerando la losa como piso rigido (Aguiar 2008). En este 
programa se introduce como dates la matriz de rigidez de la estructura (KL), la altura acumulada de cada 
piso (alt), la distancia desde el centre de masas hacia cada portico (r), el peso de la estructura (pesoD) y 
la carga viva de la misma (pesoL), con estos dates el programa calcula las fuerzas sismicas equivalentes 
del sismo rare del CEC2000, ademas genera un control P-Delta, el cortante basal, la deriva de piso y el 
indice de estabilidad de piso, con sus respectivos analisis de reforzamiento de la estructura de ser el 
caso. 

2.2 Sistema Estructural de Acero 

Se realiza el diseno del Sistema estructural de acero segun 
los perfiles estructurales (ANSI/AISC 360-05, 2005). 

El acero seleccionado para el sistema estructural es ASTM A992 con 
ksi y una resistencia ultima de 65 ksi (ANSI/AISC 360-05, 2005). 

Las columnas y vigas son las siguientes: 

Tabla 1 Cuadro de perfiles estructurales del sistema estructural 



Piso 


Coiumnas 


Vigas 


1 


W14X159 


W 16X50 


2 


W14X159 


W 16X50 


3 


W14X132 


W 16X45 


4 


W14X132 


W 16X45 


5 


W14X132 


W 16X45 


6 


W 14X82 


W 16X40 


7 


W 14X82 


W 16X40 


8 


W 14X82 


W 16X40 



la relacion demanda capacidad de 
una resistencia a la fluencia de 50 



Para el sistema de piso se utiliza los perfiles W8X10, los mismos que son direccionados de forma 
perpendicular a las costillas del deck (AISC, LRFD Manual 2005). 

El sismo puede darse ya sea en sentido X o en sentido Y, para equilibrar las inercias de las columnas se 
establece la siguiente posicion: 




Figura 2.- Distribucion de los perfiles Wide Shape utilizados como columnas 
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En cada eje se encuentran inercias fuertes y debiles haciendo que el sistema estructural se 
encuentre equilibrado para que cuando ocurra un sismo en cualquier direccion este sea resistido per las 
inercias fuertes en su mayoria en cada uno de las direcciones antes mencionadas, cabe recalcar que el 
valor del factor de reduccion de las fuerzas sismicas es de 6, (Aguiar 2007). Los valores de las fuerzas 
stsmicas equivalentes que reporta el programa modalespacialSGDLnew.nn son: 



Tabla 2.- Fuerzas laterales de sismo en sentido X, (toneladas) 



Piso 


Portico 1 


Portico 2 


Portico 3 


Portico 4 


Portico 5 


1 


5,68 


1,34 


3,54 


5,92 


2,63 


2 


3,21 


4,19 


3,69 


3,11 


3,91 


3 


4,02 


3,65 


3,31 


4,16 


3,55 


4 


2,57 


1,23 


1,68 


2,77 


1,49 


5 


1,28 


4,31 


5,95 


1,34 


4,95 


6 


4,61 


3,44 


3,26 


4,85 


3,48 


7 


7,03 


6,26 


5,59 


7,28 


6,07 


8 


14,59 


14,21 


10,19 


15,93 


12,44 



Tabla 3.- Fuerzas laterales de sismo en sentido Y, (toneladas) 



Piso 


Portico A 


Portico B 


Portico C 


Portico D 


Portico E 


1 


2,50 


5,52 


2,06 


4,83 


4,50 


2 


3,98 


3,26 


4,07 


3,42 


3,53 


3 


3,52 


4,10 


3,38 


3,98 


3,91 


4 


1,43 


2,59 


1,32 


2,33 


2,13 


5 


5,35 


3,32 


6,25 


1,64 


2,62 


6 


3,39 


4,66 


3,21 


4,37 


4,23 


7 


6,05 


7,23 


5,80 


6,98 


6,80 


8 


12,16 


15,46 


11,25 


14,84 


14,09 



Tabla 4.- Resumen de los dates obtenidos a traves del programa modalespacialSGDLnew.m. 



Piso 


Sentido X 


Sentido Y 


Deriva de 
piso 


Desplazamientos 
inelasticos en 
centre de masas 
(m) 


Deriva de 
piso 


Desplazamientos 
inelasticos en 
centre de masas 
(m) 


1 


0,0128 


0,0512 


0,0127 


0,0506 


2 


0,0183 


0,1154 


0,0182 


0,1144 


3 


0,0194 


0,1831 


0,0192 


0,1818 


4 


0,0184 


0,2474 


0,0183 


0,2457 


5 


0,0174 


0,3082 


0,0173 


0,3061 


6 


0,0214 


0,3832 


0,0209 


0,3792 
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7 


0,0189 


0,4494 


0,0183 


0,4433 


8 


0,0128 


0,4942 


0,0124 


0,4867 



En lo que respecta al diagrama de momentos debido a cargas gravitacionales este es 
simetrico, ya que se trata de una estructura simetrica en los dos ejes coordenados, come lo indica la 
figura 4. 



Es evidente que a pesar de tener una estructura simetrica en ambos ejes coordenados, la 
disposiclon de los perflles estructurales hace que el diagrama de momentos de la figura 5 no sea 
simetrico, ya que se encuentran Inerclas fuertes y deblles en el mismo portico anallzado. 



DIrecclon 
del sismo 




±— w 



I 





Figura 3.- DIrecclon Y del sismo sobre el portico central 



Como se observe en la figura 3 en la direccion del sismo existen 3 inerclas debiles y 2 inerclas 
fuertes dando como resultado el diagrama de la figura 5. 




Figura 4.- Diagrama de momentos causados por cargas gravitacionales del portico central 
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Figura 5.- Diagrama de momentos causados per cargas de sismo del portico central de la estructura 

Con las solicitaciones de momento y de cortante se empieza el analisis y el diseno del sistema 
estructural de acero, se tiene como base las Specifications for Structurai Steei Buiidings, of American 
Institute of Steel Construction, (ANSI/AISC 360-05, 2005). 

Se reaiiza un anaiisis de esbeitez bajo condiciones sismicas (ANSI/AISC 341-05, 2005), de ios 
eiementos estructuraies utiiizados como vigas y coiumnas. 







Figura 6.- Diagrama Momento-Soporte lateral adecuado 



Cuando el perfil se encuentra dentro de la zona compacta se puede alcanzar la plastificacion 
del acero, en la figura 6, (AISC, LRFD Manual, 2005), se observe las zones compacta, no compacta y 
esbelta para un perfil estructural de acero (ANSI/AISC 360-05, 2005). Un perfil estructural es compacto 
cuando: 
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L'C = 1,75 « rv " '"Z" 



(23 



Donde: Lb: Es la longitud sin soporte lateral adecuado; Lp es el punto limite de zona compacta; 
ry es el radio de giro del perfil estructural; E es el modulo de elasticidad del acero; Fy es el Ifmite de 
fluencia del acero. 

Con la plastificacion del acero se consigue un aprovechamiento total de las propiedades de 
ductilidad de este material (Vinnakota, 2006), y el memento de diseno sera: 

Mp = fy « Z (3) 

Donde: Mp es el memento plastico de diseno; Z es el modulo de seccion plastico. 



El disenar en acero en base del memento de plastificacion es el principio de diseno del LRFD 
(Load Resistence and Factor Design, ANSI/AISC 360-05, 2005). 



Se realize un analisis de demanda capacidad de vigas, considerando las fuerzas axiales de 
diseno y de solicitacion (ANSI/AISC 360-05, 2005), se precede: 



Pr 

— > 0,2 

Pc 

Pr S {Mrx Mrv\ 

-F-" -F — -] < 1.0 

Pc 9 \Mcx Mev/ - 



C^3 



Pr 

— < 0,2 
Pc 

Pr fMrx Mrv\ 

-F -F — ^1 < 1,0 



(5) 



2Pc \Mcx Mcy- 

Donde: Pr es la carga axial ultima; Pc es la carga axial resistente; Mrx es el memento ultimo 
alrededor del eje X; Mry es el memento ultimo alrededor del eje Y; Mex es el memento resistente de la 
seccion alrededor del eje X; Mcy es el memento resistente de la seccion alrededor del eje Y. 



Con las ecuaciones (4) y (5), se verifica si las solicitaciones de memento son capaces de ser 
resistidas por las vigas. De identica manera se verifica para las columnas del sistema estructural, en 
donde tambien se realize un analisis de segundo orden y se utilize las siguientes ecuaciones: 

1 

S2 = > 1 (6) 



1 - 



iPu 



ID 



Donde: B2 es el factor de amplificacion de memento P-A; EPu es la sumatoria de todas las 
fuerzas axiales ultimas del piso donde se encuentra la columna analizada; FI es la sumatoria de las 
fuerzas sfsmicas del piso donde se encuentra la columna analizada; ID es la deflexion lateral de 
entrepiso. 



Sl = 



Cm 




C"3 
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Donde: B1 es el factor de amplificacion de momento P-5; Cm es el factor de memento 
equivalente; Pu es la carga axial ultima de la columna de anallsis; Pe es la carga de pandeo elastico. 

Las ecuaclones (6) y (7), (ANSI/AISC 360-05, 2005), permiten realizar un anallsis de segundo 
orden, y obtenlendose el momento mayorado (8), se puede realizar un dlsefio plastico de las columnas 
con la ecuaclon (3), ya que todos los perflles escogidos para tal proposito son compactos. 

Mu = SI « Mnt -t- S2 « Mti (8j 

Donde: Mu es el momento ultimo; Mnt es el momento causado por fuerzas gravltaclonales; Mtl es el 
momento causado por fuerzas de sismo. 

Las vigas de piso se las disenan tomando en cuenta los recortes realizados a causa de 
equillbrar el nivel entre la viga principal y la viga de pIso (AISC, LRFD Manual, 2005), como se observe 
en la figure 7. 




Figura 7.- Recorte de la viga de piso 



2.2.1 Disefio de conexiones 

Se consideran dos tipos de conexiones, Conexiones de Cortante Simple y Conexiones de 
Momento, (ANSI/AISC 360-05, 2005). 

Las conexiones de cortante simple son aquellas que se utiliza para la conexion de la viga de 
piso y la viga principal y tan solo transmite el cortante debido a cargas gravltaclonales. 

Las conexiones de momento son aquellas que transmiten ademas del cortante el momento 
ocasionado por las cargas de sismo. 

Como se indica en la figura 8, (AISC, LRFD Manual, 2005), cuando existe una mayor 
restriccion a la rotacion de la conexion se realize una conexion de momento, mientras que a una mayor 
libertad de rotacion se establece una conexion de cortante simple. 
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Conexion de Memento FR 




Figura 8.- Grafico Momento-Rotacion para conexiones de elementos de acero 

En el presente trabajo se realiza la conexion de cortante simple con cordones de soldadura de 
acuerdo a los codigos de diseno, (ANSI/AWS D1.1, 2006), este tlpo de conexion se la realiza entre la 
viga principal y la viga de piso, el medio de union es un angulo estructural como se aprecia en la figura 9. 




Figura 9.- Conexion de cortante simple 

De igual manera se utilize soldadura para las conexiones de memento ademas de la utilizacion 
de pernos, como se observe en la figura 10, esta conexion se la realiza entre las vigas principales y las 
columnas del sistema estructural. La presencia del sismo hace que la conexion totalmente restringida 
sea necesaria, la geometrfa que se observe en la figura 10 es definida para condiciones sismo 
resistentes (ANSI/AISC341-05, 2005). 
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Los pernos que se utilizan en el presente trabajo son pernos de alta resistencia ASTM A490, 
los mismos que tienen tratamiento termico y una resistencia a la tension de 150 ksi (Vinnakota 2006). 

Para desalojar las articulaciones plasticas debidas al sismo se realizan los cortes a las vigas 
principales, como se indica en la figura 1 1 , (ANSI/AISC 358-05) 




Figura 11.- Vista de planta del detalle del recorte de los patines de la viga principal 



Los valores de la figura 1 1 se los presenta en la tabla 4: 
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Tabla 4.- Valores de los recortes realizados a los patines de las vigas principales 



Geometn'a 

(pig) 


pisos 1-2 


Pisos 3-5 


Pisos 6-8 


A 


4,50 


4,50 


4,50 


B 


12,50 


12,50 


12,50 


C 


1,60 


1,25 


1,25 


R 


13,00 


16,25 


16,25 



El sistema de piso de la presente estructura esta constituldo de un sistema metallco conocldo 
come deck y la losa de hormigon, todo este conjunto acompanado de una malla metalica, come se Indica 
en la figura 12, con esto se consigue que todo el sistema de losa funclone como un solo elemento, 
considerandolo como pIso rigido necesarlo para el anallsis sismico de la presente estructura, (Me 
Cormac, 2005). 



Malla 




Deck 



2.2.2 Disefio del Sistema estructural de acero segun la deriva de piso. 

Para el diseno del sistema estructural de acero se escogleron los sigulentes perflles 
estructurales de acero: 



Tabla 5.- Cuadro de perflles estructurales del sistema estructural 



Piso 


Columnas 


Vigas 


1 


W 14X370 


W18X192 


2 


W 14X370 


W18X192 


3 


W1 4X311 


W18X158 


4 


W1 4X311 


W18X158 


5 


W1 4X311 


W18X158 


6 


W14X193 


W 18X97 


7 


W14X193 


W 18X97 


8 


W14X193 


W 18X97 



En este diseno se observe que la deriva de piso es inversamente proporcional a la demanda 
capacidad a causa del sismo, los perflles estructurales indicados en la tabla 5, son extremadamente 
pesados, con ellos se obtiene una deriva de piso de 1,4%, que es menor a 1,5% recomendado para 
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estructuras de hormigon armado (Aguiar 2007). El disenar en base a la demanda capacidad o deriva de 
piso, dependera del fin ocupacional del sistema estructural, muchas veces el equipo que se encuentra en 
el interior de un edificio es mas costoso que la misma estructura, es por ello que un diseno en funcion de 
la deriva de piso es recomendable para cuando se tenga este proposito ocupacional. 

2.3 Diseno dei sistema estructurai de Hormigon Armado 

El hormigon elegido para el diseno del sistema estructural es de 210kg/cm^, varillas de 4200 
kg/cm^ de resistencia a la fluencia de diversos diametros. Se siguio lo prescrito en el ACI (2005). 



Tabia 6.- Cuadro de longitudes de vigas y columna s (cm) 



Piso 


Columnas 


Vigas 


1 


75X75 


40X70 


2 


75X75 


40X70 


3 


70X70 


35X60 


4 


70X70 


35X60 


5 


70X70 


35X60 


6 


65X65 


30X50 


7 


65X65 


30X50 


8 


65X65 


30X50 



A diferencia del sistema estructural de acero, en el cual se obtuvieron diagramas de momento 
anti simetricos, en esta estructura existen diagramas simetricos, ya que cada una de las columnas es de 
geometrfa cuadrada, como se indica en la figura 13. 
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Para el diseno del slstema estructural de hornnigon armado se considero una deriva de piso del 
1,25%, con los programas rlaxlnfl.m y modalespaclalSGDLnew, (Aguiar, 2008) se establecleron los 
sigulentes resultados, que se destacan en las tablas 7 y 8. 



Tabla 7.- Fuerzas laterales de sismo para la estructura en sentido X e Y (toneladas) 



Piso 


Portico 1 


Portico 2 


Portico 3 


Portico 4 


Portico 5 


1 


7,42 


7,42 


7,42 


7,42 


7,42 


2 


8,16 


8,16 


8,16 


8,16 


8,16 


3 


6,17 


6,17 


6,17 


6,17 


6,17 


4 


4,35 


4,35 


4,35 


4,35 


4,35 


5 


4,66 


4,66 


4,66 


4,66 


4,66 


6 


5,48 


5,48 


5,48 


5,48 


5,48 


7 


10,04 


10,04 


10,04 


10,04 


10,04 


8 


22,99 


22,99 


22,99 


22,99 


22,99 



Tabla 8.- Resumen de los datos obtenidos a traves del programa modalespacialSGDLnew.m 



Piso 


Sentido X e Y 


Deriva de 
piso 


Desplazamientos 
inelasticos en centre 
de masas (m) 


1 


0,0055 


0,0222 


2 


0,0082 


0,0510 


3 


0,0099 


0,0857 


4 


0,0109 


0,1237 


5 


0,0107 


0,1612 


6 


0,0123 


0,2041 


7 


0,0125 


0,2479 


8 


0,0097 


0,2818 



2.3.1 Diseno de vigas y coiumnas 

A partlr de las sollcitaclones de momento y cortante obtenidas, se Inicia el diseno de vigas y 
coiumnas del slstema estructural de hormigon armado. A partlr de las consideraclones del bloque de 
WIthney, (Romo 2008) se realize el diseno respectivo. 
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Se debe considerar los parametros del bloque de Withney indicados en la figura 14. A traves de 
la ecuaclon 9 se establece la cantidad de acero de refuerzo que se requiere de acuerdo a la geometrfa 
de la viga a disenar, (Romo 2008). 



k / 2 -Mu \ 

,45 = — « 1- : 

Fy \ « k- « c j 


m 


k = 0,85 - fc - b -d 




Donde: Fy es el limite de fluencia del acero de refuerzo; Mu es el memento ultimo; d es la 
distancia desde el centrolde de la armadura de acero a la fibra superior de la viga; b es el ancho de la 
viga; p es el factor de resistencla. 


A traves de las ecuaclones 10, 11 y 12 se calculan las cuantfas balanceada minima y maxima, 
con esto se controla que la cantidad de acero calculada con la ecuaclon 9 este dentro de los rangos 
permitidos. 


fc ( 0,003 \ 

ei) = 0,S5 « ^ 

1^4- 0,003/ 


(103 


Donde: (3 es el coeficlente para vigas rectangulares; f'c es 
hormigon; Es es el modulo de elasticldad del acero de refuerzo. 


la resistencla a la compreslon del 


14 

= — 

Fy 


(11) 


p-TTiax = 0,667 - pb 


(12) 



Donde: pb es la cuantia balanceada 



Para comparar si los esfuerzos causados por las sollcitaclones de memento son soportados 
por las dimenslones de la vIga, se precede a calcular los sigulentes parametros que se detallan en las 
ecuaclones 13, 14, 15 donde todos ellos provlenen de la figura 14. 
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= ,4s « fv 



a = 




(13) 

CH3 



J.W = r « ijc 




(15) 



Donde: T es la fuerza de traccion del acero de refuerzo; a es la altura del bloque de WIthney; M 
es el momento resistente de la secclon de la viga. 

Como el diseno de la estructura de hormigon armado, al 
metodo LRFD, entonces se establece la sigulente relaclon: 

Mu < 0 « M 

Donde: Mu es el momento ultimo; M es el momento resistente de 
a la flexion. 

Se establece la consideraclon que establece el CEC2000, en lo correspondlente a la colocaclon 
de los estrlbos de refuerzo. Para el diseno de columnas se establece un anallsis de segundo orden, con 
el fin de anallzar los efectos P-delta de la estructura. Con los mementos ultimos obtenidos se realize un 
anallsis de flexocompreslon uniaxial a traves de la ecuaclon 17 (Romo 2008). 



Mu =y:M-ux- -t-Muy- (17) 

Donde: Mux es el momento ultimo alrededor del eje X; Muy es el momento ultimo alrededor del 

eje Y. 



Igual que su similar de acero utilize el 
(le) 

la secclon; p es el factor de resistencia 



Se establecen las condiclones de apoyo de las columnas a traves de los nomogramas para 
porticos no arriostrados (Romo 2008) que se Indica en la figure 15. 
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CO 

100.0 
50.0 — 
300 — 
20 . 0 - 



9.0- 



8 . 0 - 

7.0 
0 . 0 - 

5 . 0 - 

4 . 0 - 

3.0 

2.0 



1.0 — 



0—1 



Op 

=^ 20.0 

-— 10.0 



-50 

-40 



-30 



-20 



. 1 5 



-10 



'B 

03 

100 0 

— 50.0 

— 30 0 

— 20.0 



— 10 0 
F90 
-80 
70 
-60 

-50 
-40 

30 
20 



— 1 0 



1—0 



Figura 15.- Nomograma para porticos no arriostrados 

Con el grafico de diseno de columnas (Romo, 2008), que se indica en la figura 16, se establece 
la cantidad de acero de refuerzo que las columnas de hormigon armado necesitan para soportar las 
sollcitaclones de diseno. 
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Figura16.- Grafico de diseno para columnas de hormigon armado 



De identica manera que las vigas en lo referente al cortante, se establece lo que dispone el 
CEC2000. Para el analisis de la conexion viga columna se realize un analisis de code vertical, horizontal, 
adherencia y confinamiento (Aguiar-Barbat, 1997), ademas de un chequeo de la teoria columna fuerte- 
viga debil. La losa disehada en el presente trabajo es nervada de 25cm de espesor y se utilize un 
armado minimo requerido para las solicitaciones de cargas gravitacionales. 

3. Resultados 

En lo referente a los desplazamientos inelasticos se puede observer que la estructura de acero 
posee desplazamientos mayores en alrededor de un 50% mas que su par de hormigon, como se puede 
observer en la table 9, esto predice en decidirse en una estructura de acero cuando los fines 
ocupacionales del edificio sean los mismos para ambos materiales. 

Cabe destacar que las derives de piso son distintas para ambos sistemas estructurales, esto 
se debe considerar, por condiciones de rigidez en el edificio de hormigon armado y flexibilidad en el 
edificio de acero, es decir se debe conocer las bondades que presta cade uno de estos materiales a fin 
de determiner con que derive de piso se debe disehar cada una de estas estructuras. 
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Tabla 9.- Comparacion de los desplazamientos inelasticos de los sistemas estructurales 



Acero 


Hormigon 

Armado 


Desplazamientos 
inelasticos en 
centre de masas 
(m) 

(Deriva 2,14%) 


Desplazamientos 
inelasticos en 
centre de masas 
(m) 

(Deriva 1,25%) 


0.0512 


0.0222 


0.1154 


0.0510 


0.1831 


0.0857 


0.2474 


0.1237 


0.3082 


0.1612 


0.3832 


0.2041 


0.4494 


0.2479 


0.4942 


0.2818 



Con respecto a los tiempos de ejecucion se destacar que la estructura de acero se la termlna 
en la mltad de tlempo que la estructura de hornnigon arnnado, por razones de fraguado del hornnigon. 

La mano de obra requerlda para la ejecucion de la estructura de acero es mas costosa, que su 
par de hormigon armado, como se observe en la tabla 10 los soldadores deben ser callficados con 
procedimlentos de soldadura WPS (Welding Procedure Specification) (ANSI/AWS D1.1, 2006), y esto 
conlleva a un nivel de preparaclon alto por parte de este personal, ademas se debe contratar un 
Inspector de callficaclon de soldadores CWI (Certified Welding Inspector), por su parte la ejecucion de la 
estructura de hormigon armado relativamente es barata, ya que el nivel de preparaclon de los obreros no 
exige que sea alto. 

Tabla 10.- Comparacion del costo de la mano de obra 



Hormigon Armado 


Acero 


Personal 


Costo 


Personal 


Costo 


Jefe de obreros 


Medio 


Inspector de Soldadura AWS- CWI 


Alto 


Obreros 


Medio 


Soldadores callficados AWS D1 ,1 


Alto 


Equipo 


Medio 


Equipo 


Alto 



La construccion del edificio de acero es mas costoso que el edificio de hormigon armado, esto se 
debe a que el Ecuador no es productor de acero y por ende se lo debe importer, pero este inconveniente 
se lo compensa con la rapidez de ejecucion de la obra lo que hace que la recuperacion de la inversion 
sea mas pronta que la estructura de hormigon armado. 

4. Conclusiones 



En el sistema estructural de acero se establece que las conexiones de viga a columna es la 
parte mas compleja a disehar, ya que se deben considerar conceptos de soldadura, en los cuales 
interviene clasificacion de procesos y electrodes. 
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El recorte realizado a las vigas principales del slstema estructural de acero es de suma 
Importancla, ya que a traves de ellas se desalojan las articulaclones plasticas que genera un 
deternninado sismo. 

Se puede realizar un diseno en el cual la derlva de piso sea mayor a la recomendada por el 
CEC 2000, ya que graclas a la gran ductllldad y flexibllldad del acero los desplazamlentos Inelasticos de 
cada pIso pueden ser soportados sin un mayor dano de la estructura. 

El diseno plastico realizado en este trabajo ha demostrado la gran utllldad que presta el diseho 
por resistencia y factores de carga LRFD, que ha entrado en reemplazo del ya tradiclonal diseho por 
esfuerzos admisibles, el cual consideraba un anallsis elastico el mismo que no utlllzaba del todo las 
bondades caracteristicas del acero. 

El comportamlento de la losa del slstema estructural de acero debe ser de tal forma que se 
comporte como un solo material, es declr el deck metallco y el concrete deben ser capaces de resistir el 
cortante generado por las cargas gravltaclonales, y esto se lo consigue con la colocaclon de conectores 
de cortante, los cuales se disehan de acuerdo a los codigos de diseho de AISC. 

La estructura de hormigon armado es mas rigida, y por ende no es tan flexible como su par de 
acero, esto conlleva a disehar especificamente con derivas de piso menores a lo que establece el CEC 
2000 . 



Las conexiones de la estructura de hormigon armado tienen que ser revisadas de tal forma 
que las fuerzas laterales no ocaslones dahos severos a la estructura, controles de cortantes verticales, 
horizontales, adherenclas del hormigon, son algunos de los controles que se deben efectuar a dichas 
conexiones, ademas de la verlficaclon columna fuerte-viga debll. 

El colocar acero de refuerzo en forma excesiva puede ser perjudiclal ya que no podrfa fluir de 
forma esperada ante un sismo, haclendo que las conexiones se vuelvan fraglles ya que el sismo solo 
estarfa slendo soportado por el hormigon a causa de la gran cantidad de acero que se demora en entrar 
en fluencla. 

Se debe calcular la cuantia balanceada para establecer la cuantia maxima de diseho, con el fin 
de evltar las consecuenclas descritas en el parrafo anterior. 
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RESUMEN 

Bella Rica, Namblja y Portovelo son distrltos mlneros donde se explota oro con la ayuda de 
mercurlo para la amalgannaclon. Los residuales liquidos del proceso de extracclon mezclados con 
partfculas de amalgama son ellmlnados a los cuerpos de agua. El objetivo de esta Investigaclon fue 
medlr la contamlnaclon actual per mercurlo en aguas y sedimentos de los rios, que reciben efiuentes de 
la mlneria de oro en los tres distrltos. Las concentraclones mas elevadas de mercurlo encontradas en el 
agua de los rios y piscinas de los tres distrltos mlneros fueron: en las Chanchas del rio Calera 
(Portovelo) 1812,5 p,g/L, en Bella RIca-parte baja (Ponce Enriquez) 47,88 pg/L y 29,84 |o,g/L en la 
Quebrada Calixto (Namblja). Per otra parte, los mayores contenidos de mercurlo en los sedimentos 
fueron: 440,76 pg/kg, 398,21 pg/kg, 256,88 pg/kg en Ponce Enriquez, Portovelo y Namblja, 
respectivamente. Las elevadas concentraclones de mercurlo encontradas en las aguas de los rios 
pueden ser atrlbuldas al vertido de las aguas residuales del proceso de extracclon de oro o tamblen a la 
llxivlaclon del metal desde los sedimentos por las crecldas de los rios. Los sedimentos en camblo 
acumulan mercurlo porque sus partfculas solldas activas, llgan e Inmovlllzan al mercurlo en sus 
estructuras a medida que el material particulado vlaja con las aguas de los rfos y ademas debldo al 
hidrofobicldad del mercurlo. De estos resultados, se puede declr que las operaclones mlneras que 
extraen oro son las responsables de las elevadas concentraclones de mercurlo en los medlos 
estudlados. 



ABSTRACT 



Bella Rica, Namblja, and Portovelo are mining areas where gold Is extracted using mercury for 
amalgam formation. Waste liquids from the extraction process mixed with amalgam particles are 
eliminated to water bodies. The objective of this study was to measure actual mercury contamination In 
waters and sediments of rivers that receive effluents from gold mining In the three areas. The highest 
mercury concentrations found In waters of rivers and ponds of the three mining sites were: 1812,5 pg/L at 
the Chanchas of the Calera River (Portovelo), 47, 88 pg/L In Bella RIca-lower part (Ponce Enriquez), and 
29,84 pg/L In Quebrada Calixto (Namblja). Conversely, the highest levels of mercury In sediments were: 
440,76 mg/kg, 398,21 pg/kg, 256,88 mg/kg In Ponce Enriquez, Portovelo y Namblja, respectively. The 
highest concentrations of mercury found In waters of rivers can be attributed to the release of wastewater 
from gold extraction process or due to metal leaching from sediments associated to river overflows. On 
the contrary, sediments accumulate mercury because their active solid particles link and Immobilize 
mercury within their structures as the particle material travels with river waters and besides due to 
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mercury hydrophobicity. From these results, it can be concluded that mining operations that extract gold 
are responsible for the high concentration in the investigated media. 

1 . INTRODUCCION 

En la mineria del oro, luego de que el metal es concentrado gravimetricamente, es amalgamado 
con mercurio. La amalgamacion se utiliza en la pequeha mineria primaria y en la pequeha mineria aluvial 
por ser una tecnica economica y que necesita muy pocos instrumentos de trabajo. Entre los metodos que 
se emplean para la amalgamacion se incluye, el de circuito abierto, proceso en el que todo el material 
aurffero se pone en contacto con mercurio en un flujo continue. Con esta tecnica no es posible recuperar 
todo el mercurio y una parte de este, en forma de mercurio elemental o en forma de amalgama escape 
con las agues residuales, contaminando los alrededores: El de circuito cerrado en cambio es una tecnica 
en la que solo porciones del material es puesto en contacto con el mercurio, por ejempio los 
concentrados del material que sale de los canalones o aquellos que estan a continuacion de los molinos. 
Luego de former las “bolas de amalgama”, estas son quemadas en retortas artesanales. Las retortas 
recuperan el 20% mercurio en forma de condensado y el resto se evapora o emite al ambiente 
(Donoghue, 1998). La quema de la amalgama produce varies especies de mercurio como el vapor de 
mercurio (Hg°), HgO, Hg^"" y el mercurio condensado en partfculas (Lacerda & Salomons, 1998). 

El mercurio elemental liberado a la atmosfera durante el proceso de amalgamacion y 
vaporizacion del oro se oxide a Hg^"^ a traves de reacciones en mediadas por el ozono. En las areas 
agrfcolas donde se realizan quemas de materia vegetal y bosques, se incrementa la concentracion del 
ozono en la parte baja de la atmosfera y consecuentemente puede acelerarse la oxidacion del mercurio 
mediante la siguiente reaccion Hg° +o^+2H +o^+h^o (Hall, 1995). El ion mercurio Hg^’" es 

rapidamente lavado por las frecuentes Iluvias en los bosques tropicales y penetra a los ambientes 
acuaticos y se absorbe en los sedimentos y en la materia particulada desde donde es metabolizado por 
los organismos a metilmercurio y dimetilmercurio (A. Sigel & H. Sigel, 1997). Estos compuestos son 
toxicos y se absorben con facilidad a traves de la piel, son solubles y entran a la red trofica acuatica, 
concentrandose en los predadores como los peces (USGS, 2000). El impacto en la poblacion que ingiere 
especies acuaticas con elevadas concentraciones de metilmercurio o dimetilmercurio es letal y pueden 
sufrir dahos series en el higado y en los rihones (PNUMA, 2005). La exposicion cronica a mercurio 
elemental, puede generar deficit en el desarrollo neurologico y de comportamiento (Ozuah, 2001 ; Tirado 
et al., 2000), que puede incluir dahos sutiles en la memoria visual, atencion y velocidad en las respuestas 
visuales, auditivas y psicomotoras (Counter et al., 2002, Ellingsen et al., 2001), perdida reversible de la 
capacidad para distinguir colores (Cavalier! y Gobba, 1998) y ademas de inflamaciones severas de la piel 
(Zimmer et al., 1997; Boyd et al., 2000). En un estudio de Millar (2002), se afirma que las personas que 
habitan en los alrededores de un foco de contaminacion con mercurio estan expuestas a traves de dos 
rutas de exposicion: I) inhalacion e ingestion de mercurio elemental derramado; e ii) ingestion de 
mercurio ionico, luego de la oxidacion del mercurio elemental. Por otra parte, la exposicion al mercurio 
puede causar efectos neurologicos graves en organismos superiores, como se comprobo en Minamata, 
Japon en 1950-1952. En este lugar, las aves mostraron signos de enfermedades neurologicas, incluso 
convulsiones y ataques; en los gatos cuya dieta contenia mucho pescado, se observe movimientos 
altamente erraticos (PNUMA, 2005). 

Por la alta toxicidad del mercurio, la Agenda de Proteccion Ambiental de los Estados Unidos 
(USEPA, 1974) fija como limite maximo alcanzable en agua de consume humane 2 pg/L de mercurio y el 
Texto Unificado de Legislacion Ambiental Secundaria (TULAS, 2003) establece como Ifmite maximo 
permisible 1 pg/L. Adicionalmente, en las aguas residuales que se vierten en cuerpos de agua dulce, el 
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contenido de mercurio no debe ser mayor a 5 |xg/L (TULAS, 2003). Por otra parte, la Organizacion para 
la Agricultura y Allmentaclon y la Organizacion Mondial de la Salud ha fijado como Ingesta Semanal 
Provisional Tolerable (ISPT) 0,33 pg metllmercurlo/ kg peso corporal (FAO/OMS, 1999). 

En el Ecuador la mlneria de oro tiene sus origenes en el area Portovelo-Zaruma y segun el 
estudio del Mlnisterlo de Energia y Minas (1998) comenzo en 1896. La pequeha mlneria o mlneria 
Informal, surge en la zona sur orlente del pais a finales de los ahos 70, debldo a la quiebra de grandes e 
Importantes compahfas extranjeras que explotaban oro en la zona de Zaruma y Portovelo (Sandoval, 
2001). En Nambija, con el descubrimlento del oro a principles de los 80, se Inicio la migraclon de mlneros 
segun afirma Sandoval (2001). En Bella Rica, campo mlnero enclavado en el este de la Provincia de 
Azuay, la explotaclon de oro se Inicio en 1983 (Appleton, 1996). El Mlnisterlo de Energfa y Minas (1998) 
reporta que el mercurio es usado en las areas mlneras del sur del Ecuador para extracclon del oro desde 
los concentrados medlante amalgamaclon. Consecuentemente, las aguas residuales del proceso de 
extracclon mezcladas con partfculas de amalgama son ellmlnadas a los cuerpos de agua y contamlnan 
su causes. 

El Proyecto de Desarrollo Mlnero y Control Amblental llevo a cabo el monitoreo amblental de las 
areas mlneras en el sur del Ecuador y reporta relativamente bajas concentraclones de mercurio en aguas 
de los rfos y sus sedimentos, en las cercanfas de las zonas mlneras (PRODEMINCA, 1998). 
Desafortunadamente, en ese estudio no se reporta el mecanismo de movlllzaclon y evoluclon del toxico 
en la zona de estudio. No se conoce de otros estudlos que hayan cuantificado la cantidad de mercurio 
presente en las aguas y sedimentos en esa zona mlnera. El objetivo de esta Investigaclon fue medir la 
contamlnaclon actual por mercurio de aguas y sedimentos en los rfos que reciben efiuentes de la mlnerfa 
de oro en los sectores de Nambija, Ponce Enrfquez y Portovelo y compararla con el estudio de 
PRODEMINCA. 

2. AREA DE ESTUDIO 

El estudio se desarrollo en la zona Sur del Ecuador entre 2°IT S y 4S0’S e Incluye en tres 
provinclas: Azuay, Zamora Chinchipe y El Oro. En cada provincia, se selecclono una localldad (ver Fig. 
1). El distrlto mlnero de Bella Rica (Ponce Enrfquez, Provincia del Azuay), esta ubicado a una altura 
entre 50 y 1600 m.s.n.m. y tiene cllma tropical humedo. Los yacimlentos mlnerales de Bella Rica 
contlenen pirlta, calcopirlta y arsenopirlta, mientras que el lecho rocoso esta conformado por rocas 
volcanicas andesfticas y basalticas (PRODEMINCA, 1998). El distrlto mlnero de Nambija (Provincia de 
Zamora Chinchipe) esta locallzado a una altura que osclla entre los 700 y 1400 m.s.n.m y goza de un 
cllma tropical medio seco. Posee rocas volcano-sedimentarlas con Intruslones de granodlorlta. Las 
mlnerallzaclones relaclonadas con estas Intruslones son yacimlentos Irregulares de sulfuro-skarn o 
sulfuro con vetas de cuarzo (PRODEMINCA, 1998). El distrlto mlnero Portovelo-Zaruma (Provincia del 
Oro) se encuentra entre los 1500 y 2100 m.s.n.m. y es una zona con cllma tropical humedo. En esta 
region el oro parece asoclado con vetas de cuarzo y las rocas son volcanicas dacfticas-andesiticas 
(PRODEMINCA, 1998). Los yacimlentos son complejos sistemas de vetas profundas en una zona de 
aproximadamente 50 km. de largo 10 km. de ancho y 1400 m verticalmente. 

3. MUESTREO Y PROCEDIMIENTOS ANALITICOS 
3.1 Muestreo 

Muestras de agua y sedimentos fueron tomadas de rfos y piscinas ubicados en el area de 
estudio (puntos de muestreo ver en Figs. 2 y 3). Las Investigaclones del campo Incluyeron, el 
posiclonamlento de los sitlos de muestreo en los cuerpos de agua utlllzando un sistema portatll de 
posiclonamlento global (GPS) y las medidas de los parametros de campo como pFI, temperatura y 
conductividad electrica. El muestreo de aguas Involucro la recolecclon de: I) muestras flltradas sin 
acldificar (0,45 para anallsis de aniones) y II) muestras flltradas acldificadas (para anallsis de Hg). Las 




Revista CIENCIA Volumen 12 Numero 2 de 2009 



173 



muestras fueron transportadas en un enfriador a los laboratories de ESPE para un analisis qufmico 
detallado. 

El muestreo de los sedimentos se realize tomando un peso aproximado de 0,5 kg desde un 
punto situado a 5 cm debajo de la superficie para cada sitio muestreado. Las muestras fueron luego 
almacenadas en contenedores plasticos y transportadas a los laboratories en otro enfriador. 




Figura 1. Area de estudio en la zona sur del Ecuador y localizacion de los puntos de muestreo. 

3.2 Extraccion de arsenico de los sedimentos 

Las muestras de sedimentos fueron secadas durante 12 boras en una estufa a una temperatura 
de 60°C. Los sedimentos secos fueron disgregados utilizando un mortero y pistilo y luego fueron 
tamizados con un cedazo de 425 pm de diametro. Las muestras tamizadas (250 mg) fueron sumergidas 
en 125 mL de HCI-4 N y colocadas en un agitador rotative durante 48 boras a 40 rpm. Posteriormente, 
los extractos fueron filtrados usando filtros de 0,45 mm y analizados por mercurio. 

3.3 Analisis qufmico 

Las muestras de agua acidificadas y los extractos diluidos fueron analizadas por mercurio 
usando un Generador Hidruros Perkin Elmer (HG) acoplado a un espectrometro de absorcion atomica 
Perkin Elmer, AA100, empleando el metodo de vapor en frfo, 3112 B del Standard Methods (APHA) 
modificado. Los resultados de los analisis quimicos fueron validados, comparando los contenidos de 
mercurio con las soluciones estandar compradas a Sigma Aldricb. 

4. RESULTADOS 

Con el proposito de investigar los procesos que pueden controlar la distribucion del mercurio en 
las aguas y sedimentos de rios y piscinas dentro del area de estudio, los datos estas fuentes ban sido 
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interpretados y comparados con resultados reportados por otros estudios. En general, se deternnina que 
el contenido de mercurlo en los sedlmentos y en las aguas de los rios depende su locallzaclon en las 
zonas mlneras. 

4.1 Mercuric en aguas superficiaies 

Ponce Enriquez 

En la Tabla 1 se Indica la altura, longitud, latitud y concentraclones de mercurlo total en los 
diferentes sitlos de muestreo del area de estudlo. De todas las muestras anallzadas en la zona de Ponce 
Enriquez (Bella Rica), el contenido de mercurlo mas elevado corresponds a la muestra M16 (Fig. 2A, 
Bella RIca-parte baja) con 47,88 pg/1. Las 15 muestras restantes presentan valores de Hgi Inferlores a 
los 10 (pg/L) y oscllan entre 1,90 y 8,90 pg/L. Estos valores son superlores al reportado en el estudlo 
realizado por PRODEMINCA en 1996-1998 (1 pg/L) 

Nambija 

En el distrlto mlnero Nambija, se colectaron 13 muestras de agua superficial. De las 13 muestras 
anallzadas, el valor mas alto de concentraclon de mercurlo total corresponds a la muestra M6 (29,84 
pg/1), tomada en la Quebrada Calixto (Fig. 2B). Las muestras restantes presentan valores de 
concentraclon Inferlores a 4 pg/L. En el estudlo de PRODEMINCA las concentraclones de mercurlo en la 
Quebrada Calixto se reportan ser 0.07-0.24 pg/L, mas bajas que aquellas obtenidas en la presents 
Investigaclon. 

Portoveio 

En el sector Portoveio, se colectaron 17 muestras de agua superficial, distrlbuldas en 16 sitlos 
(Tabla 1). Del total de muestras anallzadas, el valor mas alto de concentraclon de mercurlo total 
corresponds a la muestra M16 (Fig. 2C, Chanchas rio Calera) con 1812,50 pg/L, mientras que el 50% del 
resto de muestras presentan valores Inferlores a 1 pg/L (Tabla 1). 

4.2 Mercuric en sedimentos 

Ponce Enriquez 

En la Tabla 2 se Indica altitud, longitud, latitud y concentraclon de mercurlo total en sedlmentos 
de rios y piscinas dentro de la zona de estudlo. En el distrlto mlnero de Bella Rica, se observe que la 
concentraclon mas alta fue medida en la muestra M13 (440,76 pg/Kg, Fig. 3A) mientras que la 
concentraclon mas baja fue detectada es la de la muestra M12 con 3,16 pg/Kg. En la Figure 13 se 
observe que, aquellos sitlos de muestreo ubicados en la parte alta del Cerro Bella Rica tienen valores de 
concentraclon bajos y van Incrementandose a medida que se desclende el cerro. En la piscina de 
oxidaclon valores de concentraclon son medlos y altos. Del total de muestras, solo cinco muestras tienen 
concentraclones de mercurlo sobre la Norma de Calldad Amblental del Recurso Suelo (TULAS, 2003). 

Nambija 

En el distrlto mlnero de Nambija, se colectaron 16 muestras de sedlmentos. De las muestras de 
suelo sometidas a extracclon, la concentraclon mas alta corresponds la muestra M1 con 256,88 pg/Kg, 
tomada en el nacimlento del rIo Nambija (Fig. 3B), mientras que mas baja pertenece a la muestra M6 con 
0,48 pg/Kg. Como se puede ver en la Tabla 2 las concentraclones de mercurlo en las muestras se 
distrlbuyen casi uniformemente, con excepcion de dos picos bastante notorlos en las las muestras M1 y 
M13 (Fig. 3B). 
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Tabla 1. Caracterizacion de las aguas de rios y piscinas en tres regiones mineras del Sur del Ecuador 
(por localizacion de puntos de muestreo ver Figure 2). 



Muestra 


Sitio de muestreo 


Altitud 


Longitud 


Latitud 


Hgi 






(m.s.n.m.) 






(Hg/L) 


Ponce Enriquez (Bella Rica) 










Ml 


Bella Rica Parte Alta Bella Rica 


976 


9660134 


644083 


8,68 


M2 


Parte Alta 


911 


9659859 


644187 


5,94 


M3 


Bella Rica Parte Alta Bella Rica 


906 


9659835 


644184 


8,90 


M4 


Parte Media 


852 


9659857 


643845 


6,78 


M5 


Bella Rica Parte Media Bella Rica 


756 


9660416 


643439 


7,62 


M6 


Parte Media 


481 


9661404 


643159 


6,56 


M7 


Piscina de oxidacion 


433 


9661305 


642773 


6,56 


M8 


Piscina de oxidacion 


414 


9661255 


642714 


4,02 


M9 


Piscina de oxidacion 


419 


9661266 


642663 


7,20 


M10 


Piscina de oxidacion 


451 


9661225 


642756 


4,66 


Mil 


Piscina de oxidacion 


417 


9661334 


642774 


4,02 


M12 


Bella Rica Parte Baja 


43 


9661811 


640768 


4,02 


M13 


Bella Rica Parte Baja 


65 


9661888 


640734 


1,90 


M14 


Bella Rica Parte Baja 


62 


9661805 


640946 


5,72 


M15 


Bella Rica Parte Baja 


40 


9661 843 


639793 


5,50 


M16 


Bella Rica Parte Baja 


20 


9660302 


638704 


47,88 


Nambija 












M1 


Riachuelo 


1824 


9549830 


745334 


3,57 


M2 


Quebrada 1 


1425 


9551 935 


744789 


1,43 


M3 


Quebrada 2 (Camino de 


1454 


9552400 


744816 


0,19 


M4 


herradura) 


1277 


9553141 


744187 


0,14 


M5 


Quebrada 3 


1257 


9553528 


743993 


0,30 


M6 


Quebrada 4 


1190 


9554252 


743652 


29,84 


M7 


Quebrada Calixto 


1040 


9555210 


742000 


0,90 


M8 


Quebrada 7 y Rio Nambija 


1294 


9550808 


742316 


0,51 


M9 


Quebrada 1.1 


1305 


9550770 


742366 


3,29 


M10 


Playas de Nambija 


1122 


9553627 


741 742 


0,81 


M11 


Rio Nambija aguas abajo 


1090 


9554156 


741 729 


0,23 


M12 


Quebrada 1.2 


974 


9556705 


741809 


0,46 


M13 


Rio Nambija antes de Zamora 


854 


9562271 


739221 


0,57 




Gabarra 










Portovelo 












M1 


Aguas term ales 


539 


9590515 


655897 


8,59 


M2 


Rio Amarillo 


696 


9590947 


656011 


0,49 


M3 


Puente de madera 


688 


9591133 


656178 


3,25 


M4 


Quebrada baja Zaruma 


689 


9590856 


655892 


14,89 


M5 


El Qbrero 


635 


9588690 


653112 


0,28 


M6 


Rio Amarillo junto al rio Calera 


597 


9588153 


651796 


2,37 


M7 


Rio Calera cerca de confluencia 


602 


9588371 


651 844 


10,87 


M8 


Confluencia rios Amarillo y 


581 


9586847 


651 1 04 


6,24 


M9 


Calera 


543 


9586847 


650535 


0,03 


M10 


Cerca del rio Pindo 


538 


9583613 


650573 


0,03 


M11 


Rio Pindo 


595 


9589067 


651999 


0,03 
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M12 


Rio Calera cerca de chanchas 


637 


9591010 


651 832 


2,95 


M13 


Rio Calera cerca de puente Busa 


778 


9596966 


650283 


7,89 


M14 


Camino a Zaruma via Busa 


1015 


9596646 


650283 


0,20 


M15 


Rio Arcapamba 


623 


9589675 


652092 


24,57 


M16 


Canaleta cerca de chanchas 


605 


9588774 


652999 


1812,5 


M17 


Chanchas en rio Calera 
Sublendo a Zaruma 


1048 


9592455 


659320 


00,50 



639000 640000 641000 642000 643000 644000 





651000 654000 657000 660000 
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Figura 2 Localizacion de puntos de muestreo: a) Ponce Ennquez-Provincia del Azuay, b) Nambija- 
Provincia de Zamora Chinchipe, c) Portovelo-Provincia del Oro. 



Tabla 2 Mercuric total en sedimentos de rios y piscinas de tres reglones mlneras al Sur del Ecuador (per 
localizacion de puntos de muestreo ver Figura 3). 



Muestra 


Sitio de muestreo 


Altitud 

(m.s.n.m.) 


Longitud 


Latitud 


Hgr 

(tig/kg) 


Zona Minera Bella Rica 










M1 


Bella Rica Parte Alta 


976 


9660134 


644083 


10,34 


M2 


Bella Rica Parte Alta 


911 


9659859 


644187 


5,96 


M3 


Bella Rica Parte Media 


756 


9660416 


643439 


17,36 


M4 


Bella Rica Parte Media 


481 


9661404 


643159 


13,50 


M5 


Piscina de oxidacion 


433 


9661305 


642773 


257,48 


M6 


Piscina de oxidacion 


419 


9661281 


642735 


115,94 


M7 


Piscina de oxidacion 


419 


9661266 


642663 


140,84 


M8 


Piscina de oxidacion 


451 


9661225 


642756 


30,16 


M9 


Bella Rica Parte Baja 


43 


9661811 


640768 


8,94 


M10 


Bella Rica Parte Baja 


65 


9661888 


640734 


132,94 


M11 


Bella Rica Parte Baja 


62 


9661805 


640946 


4,74 


M12 


Bella Rica Parte Baja 


40 


9661843 


639793 


3,16 


M13 


Rio Siete 


20 


9660302 


638704 


440,76 


M14 


Rio Nueve 


22 


9660547 


638843 


10,88 


M15 


Chanchas Nelly 


53 


9661743 


640095 


27,72 



Zona Minera de Nambija 










Ml 


Nacimiento rio Nambija 


1824 


9549830 


745334 


256,8 


M2 


Gradas de Nambija 


1879 


9549863 


745591 


0,50 


M3 


Quebrada 1 


1425 


9551935 


744789 


62,22 


M4 


Quebrada 2 (Camino de herradura) 


1454 


9552400 


744816 


70,07 


M5 


Quebrada 3 


1277 


9553141 


744187 


49,25 


M6 


Quebrada 4 


1257 


9553528 


743993 


0,48 


M7 


Quebrada 5 


1190 


9554252 


743652 


67,81 


M8 


Quebrada 6 


1051 


9555157 


742213 


39,53 


M9 


Rio Campanas 


1040 


9555210 


742000 


53,11 


M10 


Quebrada 1 .1 


1294 


9550808 


742316 


23,66 


Mil 


Quebrada 1 .1 


1294 


9550808 


742316 


60,14 


M12 


Playas de Nambija 


1305 


9550770 


742366 


49,18 


M13 


Rio Nambija aguas abajo 


1122 


9553627 


741 742 


246,95 


M14 


Quebrada 1 .2 


1090 


9554156 


741 729 


23,44 


M15 


Rio Nambija antes de Zamora 


975 


9556705 


741 809 


35,72 


M16 


Gabarra 


854 


9562271 


739221 


73,81 


Zona Minera de Portovelo 










Ml 


Aguas termales 


539 


9590515 


655897 


19,22 


M2 


Rio Amarillo 


696 


9590947 


656011 


50,87 


M3 


Puente rio Amarillo 


700 


9591419 


656291 


72,67 


M4 


Puente de madera 


688 


9591133 


656178 


78,22 


M5 


Quebrada baja Zaruma 


689 


9590856 


655892 


85,77 


M6 


El Qbrero 


635 


9588690 


653112 


5,52 


M7 


Rio Amarillo junto al rio Calera 


597 


9588153 


651796 


2,41 


M8 


Rio Calera cerca de confluencia 


602 


9588371 


651 844 


11,02 


M9 


Confluencia rios Amarillo y Calera 


581 


9586847 


651104 


4,94 
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M10 

M11 

M12 

M13 

M14 

M15 

M16 

M17 



Cerca del rio Pindo 


543 


9586847 


650535 


20,79 


Rio Pindo 


538 


9583613 


650573 


6,83 


Rio Calera cerca de chanchas 


595 


9589067 


651999 


133,25 


Rio Calera cerca de puente Busa 


637 


9591010 


651832 


83,59 


Camino a Zaruma via Busa 


778 


9596966 


650283 


76,77 


Rio Arcapamba 


1015 


9596646 


650283 


3,93 


Pache 


623 


9589675 


652092 


1,18 


Chanchas en rio Calera 


605 


9588774 


652999 


398,21 



651000 654000 657000 660000 





651000 654000 657000 660000 
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Figura 3 Localizacion de puntos de muestreo: A) Ponce Enriquez-Provincia del Azuay, B) Nambija- 
Provincia de Zamora Chinchipe, C) Portovelo-Provincia del Oro. 



Portovelo 

En el sector Portovelo, se colectaron 17 muestras de sedimentos. Del total de muestras 
anallzadas por mercuric en los extractos del proceso de digestion, el contenido mas alto de mercuric en 
sedimento corresponde a la muestra M17 con 398,21 pg/Kg (Fig. 3C), mientras que el contenido mas 
bajo es 1,18 pg/Kg, en la muestra M16. De los resultados resumidos en la Tabla 2, unicamente 3 
muestras (M5, M12, M17) tienen contenidos de mercuric sobre los limites permisibles fijados por el 
TULAS (2003). 



5 DISCUSION 

5.1 Mercuric en ei agua superficiai de ios sectores mineros estudiados 

En Ponce Enriquez (Bella Rica), el contenido de mercuric mas elevado esta locallzado en las 
aguas del rio Slete con 47,88 pg/L (M16, Fig. 2k). El rio esta ubicado en la parte mas baja del lugar de 
estudio, por lo que las aguas de escorrentia que bajan desde las mlnas de Bella Rica pueden lavar los 
suelos y otros materlales solldos que contlenen trazas del mercuric utlllzado en la amalgamaclon del oro. 
En la piscina de oxidaclon ubicada en la parte media de Bella Rica, la concentraclon de mercuric en el 
agua esta en el range de 4,66 a 7,62 pg/l. Notar que el contenido de mercuric en el agua de la piscina no 
es muy elevado, a pesar de que las colas de Fig generadas en el proceso de amalgamlento del oro son 
descargadas hacia las aguas residuales, las que luego son transportadas a la piscina. Se cree que las 
trazas de mercuric, probablemente fueron transferldas a la fase de gas debldo a alta volatllldad del 
mercuric o fueron Inmovlllzadas por la reactividad maniflesta de los sedimentos. En la parte alta de Bella 
Rica (zona mlnera) los contenidos de mercuric disueltos en el agua son llgeramente superlores a los de 
la piscina de oxidaclon (6,78-8,68 pg/l). No obstante, de que en esta zona se lleva a cabo la 
amalgamaclon del oro con mercuric, el contenido del metal toxico disuelto en el agua es bajo. La baja 
solubllldad del mercuric (60 pg /L) (Stumm y Morgan, 1996) y la tendencia del metal por Interactuar 
rapidamente con las particulas solldas de suelos y sedimentos, hace que su contenido en las muestras 
de agua sea bajo. 

En Nambija, solamente una muestra de agua contlene mercuric sobre los 20 pg/L (M6, Fig. 2B). 
En este punto existe una significativa perturbaclon de sedimentos, lo que puede provocar llxivlaclon del 
mercuric hacia el agua de la Quebrada Calixto. El resto de muestras presentan valores bastante bajos. 
La baja concentraclon de mercuric en agua puede estar llgada a la escasa actividad mlnera existente en 
la zona y al reemplazo del mercuric por clanuro para la recuperaclon de oro desde los materlales 
rocosos. 



En el sector de Portovelo, la muestra de agua tomada en el rio Calera junto a una descarga de 
aguas residuales que provlenen de las canaletas, tiene el valor de mercuric mas alto del sector (1812,5 
pg/L) (M16, Fig. 2C). Este valor se explIca porque en las canaletas se acumulan finos que contlenen oro, 
los que posterlormente son amalgamados con mercuric. Las colas de la amalgamaclon pueden entonces 
fugar hacia la corriente acuosa y descargarse al rio Calera. Esta hipotesis se confirma al cuantificar la 
concentraclon de mercuric en una muestra de agua provenlente de una canaleta (muestra M15 = 24,57 
pg/l). El resto de muestras presentan concentraclones Inferlores a 15 pg/l, de las cuales alrededor del 
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50% exhiben valores de mercurio bajo los pernnisibles para consumo humano y use domestico y 
preservacion de la flora y fauna en aguas dulces (TULAS, 2003). 

5.2 Mercurio en ios sedimentos de ios sectores mineros estudiados 

En general los sedimentos contlenen particulas solldas activas que llgan e Inmovlllzan al mercurio en sus 
estructuras (Johannessen et al., 2005). Aunque existen desacuerdos en cuanto a los mecanismos de 
movllldad del metal en sedimentos y suelos (Rasmussen, 1994; Fitzgerald et al., 1998); sin embargo, las 
particulas solldas de los sedimentos (materia organica y fracclon Inorganica) retlenen cantidades 
significativas de Hg. El mercurio en los sedimentos esta principalmente en forma residual, como lones de 
Intercamblo, llgado a las substanclas humicas y rara vez comblnado con los carbonatos. 

En Ponce Enriquez, el contenido mas elevado de mercurio en sedimentos (440,76 pg/kg) esta 
locallzado en la orllla del rio Slete (M13, Fig. 3A). MIentras que valores entre 30,16 y 257, 48 pg/kg se 
tienen en los sedimentos de la laguna de oxidaclon ubicada en la parte Intermedia de Bella Rica. Las 
concentraclones de Fig encontradas en estos dos sitlos, confirman la tendencia de que el contenido del 
metal en los sedimentos es mucho mayor que en el agua. En la presente Investigaclon, la relaclon entre 
Hg en sedimentos y Hg en agua es mayor a 10. Sin embargo, esta es menor a la encontrada en la zona 
mlnera aurifera (200 a 600) de la cuenca del rio Tapajos, ubicada en la region amazonica brasilena 
(Telmer etal. 2006). 

La concentracion de mercurio en los sedimentos de las quebradas y rio Nambija es similarmente 
mayor a la de otros sitlos. Por ejempio en el nacimiento del rio Nambija, lugar donde se amalgama el oro 
con mercurio (M1, Fig. 3B), la concentracion es 256, 88 pg/kg mientras que en el punto de muestreo 
M15, situado a aproximadamente 6500 m aguas abajo, la concentracion desciende a 35,72 pg/kg. Este 
decrecimiento puede ser interpretado como que el mecanismo de inmovilizacion de mercurio, active la 
transferencia del mercurio disuelto en el agua del rio a los sedimentos, a medida que las aguas con 
trazas de mercurio avanzan. Esta afirmacion se puede confirmar, analizando los datos de concentracion 
de mercurio presentes en el agua del rio Nambija. El punto M1 tiene mayor concentracion (3,57 pg/L) 
comparado con la del punto M12 (0,46 pg/L) (ver Fig. 2B); sin embargo, se debe senalar que el aporte 
del agua proveniente de la Quebrada Calixto, puede provocar dilucion de los contenidos de mercurio y 
consecuentemente disminuir su concentracion en el agua en el rio Nambija. Ademas, es importante notar 
que en los puntos de muestreo M15 y M16 (Fig. 3B), antes de la confluencia del rio Nambija con el rio 
Zamora, existen contenidos significativos de mercurio en los sedimentos (35,72 y 73,81 pg/kg, Tabla 2). 
En esta zona no hay ninguna actividad mlnera cercana y la concentracion de mercurio disuelto en el 
agua en el mismo punto, es muy reducida (M13 = 0,57 pg/L, Tabla 1). Se cree que las significativas 
concentraclones de mercurio encontradas en los sedimentos del rio Nambija, son debido a su transports 
y redistribucion durante las crecidas del rio y las quebradas aledahas, las que arrastran abundantes 
sedimentos ricos en mercurio. 

En Portovelo las mas altas concentraclones de Hg, se presentan en los sedimentos del rio 
Calera (83,59 - 398,21 pg/kg, Tabla 2). Esto se explica por el hecho de que el rio recibe las aguas 
residuales provenientes del process de amalgamacion con mercurio. Investigadores chinos reportan 
similares concentraclones de mercurio en los sedimentos del rio Huangpu (70,52 pg/kg - 387,30 pg/kg), 
que igualmente recoge aguas residuales de minas de oro que usan amalgamiento con mercurio (Ding et 
al., 2005). De igual forma, los sedimentos del rio Amarillo contlenen significativas concentraclones de Hg 
(50,87- 72,67 pg/kg) porque este rio tambien es depositario de aguas residuales de las minas de oro de 
Zaruma, ubicadas aguas arriba. Comparando la concentracion de mercurio en los sedimentos con la 
concentracion del metal en las aguas de los rios del sector, se confirma que muy son superiores, 
consecuentemente los sedimentos son materiales activos para inmovilizar el mercurio que es 
transportado por las corrientes acuosas superficiales. 
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6 CONCLUSIONES 

Se puede concluir que las operaciones mineras que extraen oro con metodos tradicionales y con 
metodos sofisticados, son las responsables de las elevadas concentraclones de mercurlo. Sin embargo, 
la fuente domlnante de contamlnaclon no solo es la perdida de mercurlo en el proceso de amalgamlento, 
sino que la perturbaclon y movlllzaclon de grandes cantidades de sedimentos ricos en mercurlo y las 
Inundaclones provocadas por las operaciones mineras son tamblen fuente de contamlnaclon. Estos 
hallazgos camblan el enfoque que deberia darse a los esfuerzos de remedlaclon y prevenclon 
camblando del control al mercurlo en el sitio de uso, al control de la erosion de suelos y sedimentos. La 
mlnimizaclon o ellmlnaclon de las perdidas de mercurlo en el proceso de mlneria son Importantes para la 
salud de los habitantes y amblentes locales, pero el mantenimlento de los suelos y sedimentos 
confinados a su lugar seria una estrategla mas efectiva para mlnimizar los flujos de mercurlo a las 
microcuencas de los rios. 

■ En distrlto mlnero de Bella Rica (Ponce Enriquez), de las 16 muestras de agua anallzadas, el valor 
mas alto del contenido promedio de mercurlo total corresponde a la M16 (Bella RIca-parte baja) con 
47,88 pg/L. En cuanto a la concentraclon de Hgj en los sedimentos el range de concentraclon de 
Hgj pg/kg osclla entre 3,16 y 440,76 pg/kg. 

■ En Nambija de las 13 muestras anallzadas, el valor mas alto de concentraclon de Hgj en aguas, 
corresponde a la muestra M6 (29,84 pg/L). En lo que respecta a los sedimentos la concentraclon de 
mercurlo osclla entre 0,48 y 256,88 pg/kg. 

■ Para Portovelo de las 17 muestras de agua anallzadas, el valor mas alto de concentraclon promedio 
de Hgj corresponde a la M 16 (Chanchas rio Calera) con 1812,50 pg/L mientras que para los 
sedimentos el range osclla entre 1,18 y 398,21 pg/kg. 

■ En Ponce Enriquez, el contenido mas elevado de mercurlo en sedimentos (440,76 pg/kg) esta 
locallzado en la orllla del rio Slete (M13). Mientras que valores entre 115,94 y 257, 48 Dg/kg se 
tienen en los sedimentos de la laguna de oxidaclon ubicada en la parte Intermedia de Bella Rica. En 
la presente Investigaclon, la relaclon entre Hg en sedimentos y Hg en agua es mayor a 10. 

■ En Portovelo las mas altas concentraclones de Hg, se presentan en los sedimentos del rio Calera 
(83,59-398,21 pg/kg). Esto se explIca por el hecho de que el rio recibe las aguas residuales 
provenlentes del proceso de amalgamaclon con mercurlo. De Igual forma, los sedimentos del rio 
Amarillo contlenen significativas concentraclones de Hg (50,87- 72,67 pg/kg) porque este rfo tamblen 
es depositarlo de aguas residuales de las mlnas de oro de Zaruma, ubicadas aguas arriba. 

■ Probablemente las significativas concentraclones de mercurlo encontradas en los sedimentos del rfo 
Nambija, se deben a su transports y redlstrlbuclon durante las crecldas del rfo y las quebradas 
aledanas, las que arrastran abundantes sedimentos ricos en mercurlo. 
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RESUMEN 

En Ecuador son pocos los estudios de biodiversidad en paisajes fragmentados. Este trabajo, es 
pionero en proveer infornnacion sobre diversidad, abundancia y etnobiologia de aves que frecuentan los 
bosques riparios en tineas ganaderas y agricolas de Santo Domingo de los Tsachilas. Se registro a las 
especies mediante captura, observacion directa y rastros (sonidos). Se aplico encuestas a propietarios 
de tineas, restaurantes y cazadores para conocer el uso de las aves. Se registraron 199 individuos de 36 
especies de aves, (N=52, S=21 en tineas agricolas y N=147, S=25 en ganaderas), con 10 especies 
comunes entre sistemas productivos; se determine dos especies de aves en condiciones de amenaza de 
la UlCN y/o CITES, mientras que cinco tueron indicadoras del ambiente, nueve paraguas, cinco bandera 
y cuatro migratorias de la epoca. La abundancia, diversidad, riqueza promedio e indices ecologicos, 
entre los tipos productivos de tineas no presentaron diterencias signiticativas. Miembros de la comunidad 
ornitologica registrada no dependen estrictamente de la presencia de bosques (53%), trente a aves de 
arboledos (47%). El 42% de tinqueros le dieron uso y/o percepcion a alguna o varias especies de aves. 
El 14.3% de aves tue utilizada como alimento y el 4.1% como ornamento/mascotas. El 4.08%, 6.12% y 
20.41% de aves tuvieron la percepcion cultural, de beneticio para el ganado y como plagas de cultivos 
respectivamente. Las especies cazadas de aves, nunca tueron vendidas en el area de estudio. 

Palabras clave.- Bosque ripario, etnoornitologia, occidente ecuatoriano, intervencion antropica 

1.-INTRODUCCION 

El Ecuador es considerado como uno de los paises con mayor diversidad biologica del planeta 
(Mittermeier et al. 1997), siendo el primero en biodiversidad por unidad de superticie (Ministerio del 
Ambiente EC. et al., 2000); por esta razon, las investigaciones sobre biodiversidad, se ban visto 
plasmadas en importi antes obras (Albuja, 1999; Ridgely & Greentield; 2006; Tirira 2007). Sin embargo, 
este tipo de estudios, se limitan al estudio de la vida silvestre en areas protegidas, restando importancia 
a remanentes de zonas naturales presentes en areas alteradas como son los sitemas producctivos (Pozo 
et al., 2006; Pozo & Cardenas, 2008). 

Investigaciones toraneas en zonas con intervencion antropica, ban permitido conocer el papel 
que tienen en la conservacion de la vida salvaje, ciertos habitats de paisajes tragmentados, como cercas 
vivas, arboles disperses en potreros y bosques riberenos o riparios (Estrada et al., 1993; Maisonneuve & 
Rioux, 2001; Desebenes et al., 2003; Cardenas et al., 2004); comprobando que estos habitats son 
retugios, corredores, sitios de anidacion, reproduccion, descanso y alimentacion para la vida silvestre 
(Roth, 1976; Burel, 1996; Machtans et al., 1996; Maisonneuve & Rioux, 2001; Boutin et al., 2003; 
Desebenes et al., 2003). 
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A pesar de lo senalado, son pocas las investigaciones detalladas sobre la blodiversidad presente 
en palsajes fragmentados en toda America Latina y especificamente en el Ecuador (Pozo, 2004; 
Cardenas et al., 2004). Ultimamente, se ban desarrollado Investigaciones sobre el estado de 
conservaclon de especles sllvestres en areas productivas (Vlllacis et al., 2007; Pozo & Cardenas, 2008; 
Pozo et al., 2008; Eras, 2008; Davila 2009, Cardenas et al. en prensa). 

En Latinoamerica, se ban ejecutados varies estudlos sobre el uso dado por las personas a la 
flora y fauna (Mena, 1998). Vickers (1984), senala los estudlos realizados con las etnias Slona-Secoya, 
Cofan y HuaoranI (Vickers, 1981 ; Names & Vickers, 1982; Yost & Nelly, 1983). Las Investigaciones sobre 
los uses que comunidades Indigenas ecuatorlanas dan a los recursos naturales son muy comunes, como 
las realizadas en SInangue, Sucumbios (Mena, 1997a; Mena, 1997b), en Playa el Oro (Benitez & 
Matbeus, 1997), San Miguel (Benitez, et al., 1997), en Rio Negro Cbico (Gregory, 1997; Mena & Ruiz, 
1997) en Esmeraldas y en San Venancio, Imbabura (Gregory, 1997; Mena & Utreras, 1997). Ademas se 
ban desarrollado estudlos etnozoologicos explicitos para determlnar la Importancia economica de 
especles en diferentes etnias y comunidades (Mena, 1996; 1998). 

Como puede notarse, estudlos etnoblologicos, en comunidades productivas agropecuarlas 
ecuatorlanas no existen (Mena et al., 1997), dejando de lado colectividades no Indigenas, que usan los 
recursos disponibles de la naturaleza, estableclendo de este mode una Importancia y conocimlento 
antropologico, que podrian perderse sin los aportes de trabajos como el presente realizado en Santo 
Domingo de los Tacbllas - Ecuador (TSD). 

La presente Investigaclon es plonera, en lo que corresponde el estudio de la conservaclon de 
blodiversidad en habitats fragmentados de sistemas productivos; se pretende proveer Informaclon sobre 
la diversidad, abundancia, rlqueza y uses de aves sllvestres, alojados en los bosques riparlos de las 
fincas agricolas y ganaderas del bosque muy bumedo Pre-montano (bmbPM) del Noroccldente 
ecuatorlano. Ademas, espera proporclonar Informaclon util para sostener y promover producclones 
limplas 0 amigables con el amblente; contrlbuyendo de este mode con el mantenimlento de los bosques 
riparlos y, al control de la expansion de la frontera agricola. 

2.- MATERIALES Y METODOS 

Area de estudio 

La presente Investigaclon se realize en SDT la que se ubica entre los 500 a 700 m s.n.m., a 133 
km de Quito; entre la latitud S 0- 207 S 0- O'y la longitud W 79- 307 S 79- O' (IGM, 2007). SDT presenta 
una precipitacion media anual de 3 005 mm, bumedad relativa promedio de 85,8% y temperature media 
de 24,2-C con poca variacion durante el dia (INIAP, 2004). 

La zona de vida del area de estudio es el Bosque muy bumedo Pre-Montano (bmbPM) 
(Canadas, 1983), perteneciente al Piso Zoogeografico Tropical Noroccidental (Albuja et al., 1980), con 
dos formaciones vegetales: la sabana formada principalmente por arboles disperses en pastizales (Pozo 
& Cardenas, 2008) y el bosque bumedo siempre verde del occidente ecuatorlano (Sierra et al., 1999), 
cuyos remanentes sobreviven como bosques riparlos (Vlllacis et al., 2007; Pozo & Cardenas, 2008). 

Epoca de estudio 

Los dates fueron tornados en la epoca Iluviosa, puesto que la diversidad y abundancia de aves, 
son mayores que en la menos Iluviosa (Mena, 1997a; 1997b). La mayor rlqueza y abundancia de fauna 
en la epoca Iluviosa, se debe principalmente a la abundancia de alimento; en dicha epoca se ban 
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registrado especies de aves y manniferos que han estado ausentes en la menos Iluviosas (Gregory, 
1997; Mena, 1997a; 1997b; Mena & Ruiz, 1997). 

Diversidad de aves 

Al azar se seleccionaron 6 tineas (3 agricolas y 3 ganaderas), con remanentes de bosques 
natives riparios (Tabla 1). Las salidas de campo se realizaron en seis meses consecutivos (de junio a 
noviembre del 2007); los periodos de monitoreos fueron 48 h nnes-1 (Pozo & Cardenas, 2008; Cardenas 
et al., 2004). 



Tabla 1. Periodos de estudio y tineas monitoreadas, Santo Domingo de los Tsachilas, 2007 



Orden 


Nombre de Finca 


Tipo de 
Finca 


Parroquia o 
Recinto 


Fecha de Sal Ida 


1 


Esperanza I 


Agrtcola 


"Luz de America" 


8 junio 2007 


2 


Esperanza II 


Agrtcola 


"El Estuerzo" 


13julio 2007 


3 


El Rancho 


Ganadera 


"Luz de America" 


3 agosto 2007 


4 


Zoila Luz 


Agrtcola 


"Luz de America" 


20 agosto 2007 


5 


San Antonio 


Ganadera 


"km 35.5 " 


19 septiembre 2007 


6 


Luts 


Ganadera 


"El Estuerzo" 


16 noviembre 2007 



Las aves se colectaron con la ayuda de 4 redes de neblina de 6 m de longitud por 2.7nn de alto y 
de 35nnnn de mesh. Las redes se colocaron a lo largo de un sendero de 200 m de bosque ripario, y 
tueron templadas entre las 5:00 y 10:00h y entre las 14:00 hasta 18:00h. Los monitoreos y captures se 
realizaron en base de los patrones de actividad de las aves (Pozo & Cardenas, 2008). Un representante 
de cada especie de ave capturada tue totogratiado de pertil, dorso y vientre e inmediatamente liberado, 
ningun ave tue sacriticada, ya que se pretendio tundamentar el criterio de conservacion en los tinqueros. 

Se realizaron tambien recorridos de observacion por tinea, donde se reconocieron chillidos, 
gorgoteos, cantos y sonidos (Coopmans et al., 2004) los mismos que tueron guardados en cintas 
magnetotonicas, ademas se realizaron observaciones cuando las aves se posaron o levantaron vuelo 
con la ayuda de binoculares 10x40. Estos registros tueron bosquejados en una libreta de campo de 
escritura resistente al agua. 

Se obtuvieron registros indirectos, por medio de encuestas a los propietarios de las tineas 
productivas, de especies que no se pudieron capturar, ni registrar en las salidas de campo y sus 
monitoreos; de estas aves tambien se detallaron bosquejos; estos registros no se incluyeron en el 
analisis estadtstico pero si en el etnobiologico. 

Todas las especies registradas, se identiticaron con la Guta de Campo de las Aves del Ecuador 
(Ridgely & Greentield, 2006). Los registros auditivos, se identiticaron por medio de la guta de audio de 
Aves del Occidente Ecuatoriano (Coopmans et al., 2004). 

En 26 tineas se aplicaron encuestas sobre uses de aves para conocer los porcentajes de las 
especies mas importantes segun las utilidades y percepciones mas comunes dadas. 

Analisis de datos 

Indices de biodiversidad de Margalet, Simpson, Shannon y Hill tueros utilizados (BIODIVERSITY 
PRO, 1997; Moreno, 2001), el diseho experimental tue un DCA, para aceptar o rechazar la hipotesis se 
utilize las pruebas de Fisher con alta menor o igual a 0,05 (INFOSTAT, 2006). 
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3.- RESULTADOS 



Diversidad y abundancia de aves 

En ambos sistemas productivos, se registraron 199 individuos, de los cuales 52 (26 %) se 
encontraron en fincas agricolas y 147 (74 %) en fincas ganaderas. Estos individuos corresponden a 36 
especies de aves; de estas especies 21 fueron registradas en fincas agricolas y 25 en fincas ganaderas, 
siendo 10 comunes en ambos tipos de fincas (Porcentaje de similitud de Sorensen= 42,5 %). 



En los dos tipos de sistemas productivos se registraron 8 ordenes, 1 6 familias y 36 especies. Los 
ordenes mas abundantes en especies fueron Passeriformes (19), Apodiformes (5), Psittaciformes (3) y 
Columbiformes (3) (Tabla 2; Anexo A). 



Tabla 1. Fa milias y ordenes de aves por cada sistema productivo de Santo Do mingo, 2007. 



Ordenes 


Familias: 

Fincas Agrfcolas 


Fincas Ganaderas 


Apodiformes 


Trochilidae 


Trochilidae 


Ciconiiformes 


- 


Ardeidae 


Columbiformes 


Columbidae 


Columbidae 


Cuculiformes 


- 


Cuculidae 


Falconiformes 


Accipitridae 


- 





Emberizidae 






Fringillidae, 


Emberizidae 




Furnariidae 


Furnariidae 


Passeriformes 


Mimidae 


Hirundinidae 




Parulidae 


Thraupidae 




Thraupidae 

Tyrannidae 


Tyrannidae 



Piciformes 


Picidae 

Ramphastidae 


Ramphastidae 


Psittaciformes 


Psittacidae 


Psittacidae 



Tabla 2. Numero de familias y especies por ordenes de aves, en los dos sistemas productivos de Santo Domingo, 

2009. 



Ordenes 


Familias 


Especies 


Apodiformes 


1 


5 


Ciconiiformes 


1 


1 


Columbiformes 


1 


3 


Cuculiformes 


1 


1 


Falconiformes 


1 


1 


Passeriformes 


8 


19 


Piciformes 


2 


2 


Psittaciformes 


1 


3 


Total* 


16 


35 
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*No incluye a un individuo registrado 
como desconocido 

En las tineas agncolas se registraron 6 ordenes con 13 familias, mientras que en las ganaderas 
se registraron 11 familias de 6 ordenes. Se encontraron en comun 5 ordenes (Apodiformes, 
Columbiformes, Passeriformes, Piciformes y Psittaciformes) y 7 familias (Emberizidae, Furnariidae, 
Thraupidae, Tyrannidae, Psittacidae, Trochilidae y Ramphastidae) en ambos sistemas productivos (Tabla 
2 ). 



En las fincas agricolas se registraron 20 especies con 50 individuos, de las cuales la mas 
abundante fue Sporophila corvina seguida de Pteroglossus erythropygius y Celeus loricatus, las 17 
restantes tuvieron una abundancia menor a 3 individuos (Figura 1). 
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Especies 



□ Sporophila coruina 

m Pteroglossus erythropygius 

□ Celeus loricatus 

□ Fluvicola nengeta 

■ Sinallaxis brachyura 

□ Buteo nitidus 

□ Empidonax virescens 

□ Euphonia saturnata 
m Forpus coelestis 

□ Geothlypis auricularis 

□ Arremon aurantiirostris 

□ Carduelis xanthogastra 
m Columbia livia 

m Leucopternis occidentalis 

□ Mimus longicaudatus 

■ Myioborus miniatus 

□ Oryzoborus angolensis 

□ Phaetornis malaris 

□ Thraupis episcopus 

□ Tytira inquisitor 

□ Desconocido 



Figura 1. Ranking abundancia de especies de aves en bosques riparios de fincas agricolas, Santo Domingo, 2009. 



En las fincas ganaderas se registraron 148 individuos pertenecientes a 25 especies. Bubulcus 
ibis, Pteroglossus erythropygius, Fluvicola nengeta, Sporophila corvina, Claravis pretiosa y Oryzoborus 
angolensis, fueron las especies mas abundantes, las 19 restantes especies presentaron una abundancia 
inferior a 10 individuos (Figura 2). 



Las especies Empidonax virescens, Fluvicola nengeta, Forpus coelestis, Leucopternis 
occidentalis, Oryzoborus angolensis, Phaetornis malaris, Pteroglossus erythropygius, Synallaxis 
brachyura, Sporophila corvina y Thraupis episcopus, se encontraron tanto en las fincas agricolas como 
en las ganaderas. 



El ave cantante mas evidente en las fincas productivas fue Celeus loricatus (Carpintero Canelo), 
su voz es muy distintiva, compuesta por un chillido ruidoso y seguido de un largo borboteo, «iscuiial 
Klu-klu-klu-klu-klu-klu» (Ridgely & Greenfield, 2006). 
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Nuevos registros 



Durante este estudio se capture dos individuos de Phaetornis malaris, en dos diferentes tineas, 
una de tipo agrfcola y la otra ganadera. Con lo que se puede aseverar, que esta especie se encuentra 
tanto en las bajuras y estribaciones del este come en las bajuras del oeste, especialmente en SDT hasta 
los 1000 m s.n.m. (Pozo & Eras, 2007). 



35 n 



30 



25 - 




Especies 



□ Bubulcus ibis 

m Pterogiossus erythropygius 

□ Fluvicola nengeta 

□ Sporophila coruina 
m Claravis pretiosa 

□ Oryzoborus angolensis 

□ Thraupis episcopus 

□ Ramphocelus icteronotus 

m Neochelidon tibialis minima 
m Empidonax virescens 

□ Forpus coeiestis 

□ Riparia riparia 

m Tachyphonus rufus 
m Columbina buckleyi 

□ Crotophaga sulcirostris 

■ Leucopternis occidentalis 

□ Notiochelidon cyanoleucag 

□ Sinallaxis brachyura 

□ Amazilia amabilis 

□ Amazilia tzacatl 

□ Damophila julie 

□ Pachyramphus spodiurus 

□ Phaetornis malaris 

□ Phaetornis yaruqui 



Figura 2. Ranking de abundancia de especies de aves en bosques riparios de tineas ganaderas, Santo Domingo, 

2009 . 



Especies no capturadas 

Los propietarios de las tineas encuestadas, senalaron la presencia de 7 especies registradas y, 
de 14 mas que no se pudieron capturar ni observar en los monitoreos (Anexo C). 

Especies en riesgo 

Leucopternis occidentalis (Gavilan Dorsigrfs) tue la unica especie que se encuentra en la 
categorfa <En Peligro al igual que otras 15 especies en el Ecuador, lo que signitica, que se encuentra 
en alto riesgo de extincion en la siguiente o proxima decada; Pachyramphus spodiurus (Cabezon 
Pizarroso), se encuentra <AI Margen de Riesgo, por lo tanto, a pesar de que no este seriamente 
amenazada, su estado es preocupante, por lo que precise de cuidados tuturos (Ridgely & Greentield, 
2006). 
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Areas endemicas para las aves registradas 

Se deternnino nueve especies pertenecientes a tres de los centres endemicos (Ridgely & 
Greendfield, 2006); cuatro especies endemicas a las <Bajuras de Choc6>, siendo esta el area 
comprendida per el anillo de bosque pluvial del Choco (oeste de Colombia y noroeste de Ecuador); cinco 
especies del centre endemico <Bajuras Tumbesinas>, que se extiende desde el sur del Ecuador (Manabi 
y Los Rios), hasta fragmentos al norte de Esmeraldas y el noroeste del Peru. C. buckleyi y G. auricularis 
fueron especies circunscritas a las <Bajuras tumbesinas>, y ademas al <Rio Maran6n> (Tabla 3). 

Tabla 3 Especies registradas, pertenecientes a areas endemicas suscritas al territorio ecuat oriano 



Area endemica 


Nombre Vuigar 


Especie 




Loro Cachetirrosa 


Pionopssitta puicra 


Bajuras de Choco 


Ermitaho Bigotiblanco 


Phaethornis yaruqui 


Dios te de 


Pterogiossus erythropygius 




Tucan del Choco 


Ramphastos brevis 




Tortolita Ecuatoriana** 


Coiumbina buckieyi 


Bajuras 


Gavilan Dorsogris** 


Leucopternis occidentaies 


Tumbesinas* 


Periquito del Pacifico** 


Forpus coeiestis 




Urraca Coliblanca** 


Pachyramphus spodiurus 




Antifacito Lorinegro** 


Geothiypis auricuiaris 


Rio Marahon 


Tortolita Ecuatoriana 
Antifacito Lorinegro 


Coiumbina buckieyi 
Geothiypis auricuiaris 



* Area endemica a la que pertenece SDT, ** Especies endemicas de SDT 



Las areas endemicas senaladas y sus limites, corresponden a las nueve propuestas por Ridgely 
& Greenfield (2006), las que son adaptaciones a las propuestas por el Directorio Global de Areas 
Endemicas para las Aves (A Global Directory of Endemic Bird Areas, 1997) de Birdlife International. 

Aves indicadoras, paraguas y banderas 

Las especies consideradas como indicadoras del ambiente, fueron las aves endemicas del area 
de estudio (Gregory, 1997). SDT al encontrarse ubicado en el area endemica <Bajuras Tumbesinas> 
(Ridgely & Greenfield, 2006), present© 5 especies indicadoras del ambiente (Tabla 3). Este centre 
endemico tiene 59 especies endemicas en total. 

Accipiter superciliosus, Buteo nitidus, Leucopternis occidentales, L. plumbea y Pteroglossus 
erythropygius, por su condicion de rapaces, pueden ser consideradas tanto como paraguas como 
bandera (ANCON, 2002; Guerrero et al., 2006 Del Moral, 2007). 

Ademas, como especies bandera, en vista de su atractivo y carisma, frente a la comunidad 
biologica, se pueden considerar a: Phaethornis malaris, por su nuevo registro territorial; a Pionopssitta 
pulcra y Pionus chalcopterus, por la importancia frente a la comunidad de SDT; y a Ramphastos brevis y 
Pachyramphus spodiurus. 

Especies endemicas 

La unica especie endemica del Ecuador encontrada en las fincas productivas fue Pteroglossus 
erythropygius, conocida generalmente como “Arasari Piquipalido” o “Dios Te De”. Esta especie tiene un 
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ambito hogareno bastante amplio a traves de remanentes de bosques y arboledos del occidente 
ecuatoriano y pese a la decllnaclon de sus poblaclones a causa de la deforestaclon, muestra una 
supuesta tolerancia a la fragmentaclon y destrucclon de su habitat, per lo aparentemente no se halla bajo 
amenaza. 

Especies migratorias 

A pesar, de que el Ecuador al situarse sobre la linea equinocclal, este distanclado de los habitats 
de la mayoria de aves migratorias, clertas especies visitan al pais desde los cuatro puntos cardinales 
(RIdgely & Greenfield, 2006). En la presente Investigaclon se registraron tres migratorias boreales y una 
austral (Tabla 4). 

Tabla 4.Tipo, especies y condicion de criaderas de las aves migratorias registradas, Santo Domingo, 2009. 



Nombre cienti'fico 


Nombre Vulgar 


Tipo de Migratoria 


Criaderas en 
el Ecuador 


Bubulcus ibis 


Garceta Bueyera 


Boreal 


Si 


Empidonax virescens 


Mosquitero Verdoso 


Boreal local ista 


No 






Boreal (sep-nov; mar- 


Mn 


Riparia riparia 


Golondrina Ribereha 


abr) 


InU 


Notiocheiidon cyanoieuca 


Golondrina Azuliblanca 


Austral 


Si 



Las migratorias boreales, crian en el hemisferio norte y migran hacia el Ecuador durante el 
invierno del norte. Las migratorias australes, cria en el hemisferio sur y migran hacia el Ecuador. B. ibis al 
igual que N. cyanoleuca, a ciencia cierta cuenta con poblaclones migratorias como poblaclones criaderas 
en el Ecuador (RIdgely & Greenfield, 2006). 

Local izacion de las especies 



El 33 % de las aves inventariadas dependen totalmente de los bosques y remanentes del bosque 
muy humedo Pie Montano; el 14% no dependen plenamente, ya que habitan tanto dentro como fuera de 
areas boscosas; 38% (19 especies) dependen de los bosques fragmentados, fincas ganaderas y 
agricolas; 4% de especies habitan netamente en bosques riparios y a lo largo de los rios; 8% especies 
no viven ni dependen de los bosques y habitan en bosques secundarios y otros lugares; 8% son aves 
que habitan en zonas urbanas (poblados y casas) (RIdgely & Greenfield, 2006). 

En definitiva, las comunidades de aves, mostraron una menor proporcion de especies que 
dependen, facultativamente o estrictamente, de bosques y arboledos (23 especies, 47%); que de 
especies no dependientes en absolute de bosque (26 especies, 53%), estas especies incluyen a 
especies generalistas, de zonas agricolas, ganaderas, urbanas y fragmentadas (Figura 4). 
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Habitats 



Figura 4. Distribucion de las aves inventariadas per su localizacion (Adaptado de Ridgely & Greenfield, 2006), 

Santo Domingo, 2009 



Abundancia de aves 

Se hallaron en promedio 10.33 (± 1.98) especies y 7.43 (± 3.25) individuos por finca productiva. 
La mayoria de individuos presentes en los bosques riparios se concentraron en 9 especies de 7 families, 
que representaron el 72.36% del total de individuos (Tabia 5). La garza (Bubuicus ibis), ei arasari 
(Pteroglossus erythropygius) y el espiguero (Sporophila corvina) con 99 individuos representaron el 
49.76% del total, las otras 6 especies restantes de ia tabia 5 representaron 22.63%. Las restantes 26 
especies fueron poco abundantes (Figures 5 y 6). 



Tabia 5. Especies mas comunes de aves, nombre vulgar, familia, total de individuos y numero de fincas con la 
especie, Santo Domingo de los Tsachilas, Ecuador, 2009 



Nombre Vulgar 


Especie 


Famllla 


N 


% 


Fincas 
con la 
especie* 


Garceta Bueyera 


Bubuicus ibis 


Ardeidae 


30 


15,08 


1 


Arasarf Piquipalido 


Pteroglossus erythropygius 


Ramphastidae 


26 


13,07 


4 


Espiguero Variabie 


Sporophila corvina 


Emberizidae 


25 


12,56 


5 


Tirano de Ague 


Fluvicola nengeta 


Tyrannidae 


18 


9,05 


5 


Tortolita Azul 


Claravis pretiosa 


Columbidae 


11 


5,53 


1 


Semiiiero menor 


Oryzoborus angolensis 


Emberizidae 


11 


5,53 


3 


Tangara azulejada 


Thraupis episcopus 


Thraupidae 


10 


5,03 


3 


Tangara lolilimon 


Ramphocelus icteronotus 


Thraupidae 


7 


3,52 


2 


Carpintero Canelo 


Celeus loricatus 


Picidae 


6 


3,02 


3 


Total 






144 


72.36 





* Numero de fincas con ias especie de 6 en total 
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La especie mas abundante en las fincas agricolas fue el Espiguero Variable (Sporophlla corvina). 
En fincas ganaderas la Garza o Garceta Bueyera (Bubulcus Ibis) fue la especie mas abundante. El Dlos- 
te-de (Pteroglossus erythropyglus) en ambos tipos de fincas productivas fue la segunda ave en 
abundancla. 

Curvas de rango-abundancia de aves 

Los 3 grupos abundancla proporclonal (PI) de las especles de aves registradas por sistema 
productive se muestra en la tabla 6. PI de las fincas agricolas fue 0.048 (± 0.011) y la desviaclon 
estandar (ds) Igual a 0.050, mientras que de las fincas ganaderas fue 0.040 (± 0.010) y ds= 0.048 
respectivamente. 



Tabla 6. Ranges de abundancla de las aves registradas en las fincas productivas, Santo Domingo, 2007. 



Ranges de Abundancla 




Fincas Agrfcolas 
Maximo Mmimo 


Fincas Ganaderas 
Maximo Mmimo 


Abundante (A) 


x+ 16/34( c/s) 


1 


0,0712 


1 


0,0628 


Medianamente Abundante (MA) 




0,0712 


0,024 


0,0628 


0,0172 


Poco Abundante (PA) 


X- 16/34( c/s) 


0,024 


0 


0,0172 


0 



Las especles pertenecientes a cada range de abundancla, por sistema productive se muestran 
en la tabla 7. 



Tabla 7. Especles por ranges de abundancla y por sistema productivo, Santo Domingo, 2009. 



Ranges de 
Abundancla 


Abundante 

(A) 


Medianamente 
Abundante (MA) 


Poco Abundante 
(PA) 


Especles 

Fincas 

Agrfcolas 


Sporophila corvina 

Pteroglossus 
erythropygius 
Celeus loricatus 


Fluvicola nengeta 

Synallaxis brachyura 
Buteo nitidus 

Empidonax virescens 
Euphonia saturnata 
Forpus coelestis 

Geothlypis auricularis 


Arremon aurantiirostris 
Carduelis xanthogastra 
Columbia livia 

Leucopternis occidentalis 
Mimus longicaudatus 

Myioborus miniatus 

Oryzoborus angolensis 
Phaethornis malaris 

Thraupis episcopus 

Tytira inquisitor 


Especles 

Fincas 

Ganaderas 


Bubulcus ibis 

Pteroglossus 
erythropygius 
Fluvicola nengeta 

Sporophila corvina 

Claravis pretiosa 

Oryzoborus angolensis 


Thraupis episcopus 
Ramphocelus 
icteronotus 
Neochelidon tibialis 
Empidonax virescens 
Forpus coelestis 

Riparia riparia 

Tachyphonus rufus 


Columbine buckleyi 

Crotophaga sulcirostris 
Leucopternis occidentalis 
Notiochelidon cyanoleuca 
Synallaxis brachyura 

Amazilia amabilis 

Amazilia tzacatl 

Damophila julie 

Pachyramphus spodiurus 
Phaethornis malaris 

Phaethornis yaruqui 

Pionus chalcopterus 
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La abundancia proporcional (Pi) de las 21 especies registradas en las tineas agricolas indico que 
Sporophila corvina, Pteroglossus erythropygius y Celeus loricatus fueron las aves mas abundantes (entre 
Pi 1 y 0.071). Entre las especies medianamente abundantes (entre Pi 0.071 y 0.024), se encuentran 7 
Fluvicola nengeta, Synallaxis brachyura, Buteo nitidus, Empidonax virescens, Euphonia saturnata, 
Forpus coelestis y Geothlypis auricularis. Las aves menos o poco abundantes (entre Pi 0.024 a 0) tueron 
10 especies (7). 

La curva de las abundancias relativas (Pi) de las tineas ganaderas muestra que las unicas 
especies similares con las respectivas de las tineas agricolas en este range tueron Pteroglossus 
erythropygius y Sporophila corvina; con media abundancia Thraupis episcopus y Empidonax virescens, 
estuvieron en este range de abundancia en ambos tipos de tineas productivas. 

Porcentaje de similitud de especies de aves 

Se determine 43% de similitud de especies entre lo dos tipos de tineas; lo que indica 
cualitativamente que de cada 10 especies de los bosques riparios 4 son comunes entre ambos tipos de 
sistemas productivos. 

Entre tineas, los porcentajes de similitud, entre tineas agricolas indica que de 10 especies de 
aves, alrededor de 2 aves tueron comunes; mientras que entre tineas ganaderas 3 de cada 10 especies 
de aves se encontraron en comun. 

Conservacion de aves siivestres entre sistemas productivos 

La riqueza (p=0.4967), la diversidad (p=0.4501) y la abundancia de aves, no presentaron 
diterencias signiticativas entre tineas agricolas y ganaderas (F1,4= 1.59, p= 0.2757, Tabla 8). A pesar, de 
que la curva de acumulacion de especies (Figure 5) sugiere que las tineas ganaderas son las mas 
diversas. 

Tabia 8. Promedio (±Se), de Riqueza (S), abundancia (A/) y diversidad de aves registradas en las tineas productivas 
y de acuerdo al sistema productivo, Santo Domingo, 2009. 



Variables 


Fincas Productivas 
(n=6) 


Fincas Ganaderas 
(n=3) 


Fincas Agrtcolas 
(n=3) 


CV 

(%) 


Abundancia (A/) 


7.43 ± 3.25 


17.00 + 10.50a 


49.33 ± 23.10a 


52.37 


Riqueza (S) 


10.33 ±1.978 


7.43 + 3.25a 


9. 00 +3. 512a 


26.26 


Diversidad (/-/') 


0.872 ± 0.055 


0.92 ± 0.03a 


0.83 ±0.1 1a 


8.16 



Letras distintas indican diterencias signiticativas (p<= 0.05). 



La curva de acumulacion de especies de aves (Figure 5), indica que tanto en tineas agricolas 
como en ganaderas, existio un progresivo aumento de las especies contorme aumentaban los individuos 
inventariados. 

Las curves no se saturan al tinalizar el estudio, lo que indica que taltaron especies por ser 
registradas, y por ende al no llegar a ser horizontal taltan campahas de muestreo por sistema productivo. 
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AGRICOLAS 0 GANADERAS | 

Figura 5. Curva de acumulacion de especies de aves en los 
bosques riparios per sistema de productive, Santo Domingo, 2009 

De igual manera, no se observaron diferencias entre fincas ganaderas y agricolas, en los indices 
de diversidad Simpson (p= 0.4214) y su recfproco (p= 0.5647), al Igual que la abundancia reflejada a 
traves de la Equitatividad (p=0.3950) y la rlqueza expresada medlante el indIce de Margalef (p= 0.9307) 
(Tabla 9). 

Tabla 9. Promedio (±Se) de Indices de Simpson (D), Equitatividad (£), Margalef (Dwg) y Reciproco de Simpson (D ') 
de aves, d e acuerdo al sistema productive, Santo Domingo, Ecuador, 2007. 



Variables 


Fincas Agn'colas (n=3) 


Fincas Ganaderas (n=3) 


Indice de Simpson (D) 


0,07 ±0,22® 


0,11 ± 0,03a 


Indice de Equitatividad (£) 


0,94 ± 0,04® 


2,58 ± 1,69a 


Indice de Margalef (DMg) 


6,85 ± 1 ,34® 


6,87 ± 0,50a 


Recfproco de Simpson (D ') 


15,1 7 ±3,32® 


1 2,06 ± 5,22a 



Letras distintas indican diferencias significativas (p<= 0.05). 



En consecuencia estadisticamente, no existieron diferencias entre la conservaclon de los 
bosques riparios de ambos tipos de producclon, Indicando una Igual perturbaclon ecologica para las 
comunidades de especies de aves que habitan en dichos bosques. 

Etnobiologia de las aves 



El 42% de los finqueros (11), le dieron uso y/o percepcion a alguna o varlas aves del Inventarlo 
de este estudio, con lo que 24 especies de aves (49%), presentaron uno o varies usos/percepciones 
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(Figura 6). Ninguna especie de ave ha side comercializada per los propietarios de las fincas productivas 
despues de ser cazada. 




a) 



Aves con usos y/o Aves sin usos y/o 
percepciones percepciones 




Conocen usos y/o 
percepciones 



Figura 6. a) Especies de Aves registradas con usos y/o percepciones. b) Finqueros que conocian usos y/o 
percepciones a aves, Santo Domingo, 2007. 



Usos de ias aves 



Los usos determinados para ciertas aves se restringieron al de alimento y el de ornamento o 
mascotas tan solo 9 especies se enmarcaron con estas utilidades, representando el 22.4% del total de 
especies inventariadas y el 52.38% de especies con un uso y/o percepcion. 

Las especies de aves que poseen utilidad se resumen en la Tabla 10. Pionopssitta pulchra y 
Pionus chalcopterus, son especies susceptibles de ser cazadas como mascotas, ya que se conoce que 
los Tsachilas (Indios colorados) las persiguen por este motivo, sin embargo esta informacion no se 
obtuvo de los finqueros de la zona, por esta razon, no se los incluyo en el analisis etnobiologico 
de datos. 



Tabla 10. Especies de aves con utilidad para los finqueros de la zona. Santo Domingo, 2007. 



Nombre Vulgar 


Especie 


Usos 


Azulejo 


Thraupis episcopus 


Ornamento 


Loro alibronceado 


Pionus chaicopterus 


Alimento 


Loro/perico 


Pionopssitta puichra 


Alimento 


Paloma 


Coiumba subvinacea 


Alimento 


Paloma 


Coiumbia iivia 


Mascota 


Perdiz 


Geotrygon montana 


Alimento 


Periquito del pacifico 


Forpus coeiestis 


Alimento, Ornamento 


Tortolita Ecuatoriana 


Coiumbina buckieyi 


Alimento 


Tucan 


Ramphastos brevis 


Alimento 



Las especies destinadas como ornamento en SDT, fueron Thraupis episcopus, Columbia livia y 
Forpus coelestis. 
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Percepciones de las aves 

Veinte y dos especies tuvieron una o varias percepciones, representando el 42% del total de 
especies registradas y el 92% de las especies con uso y/o percepcion (Tabla 11 y 12). 

Tabla 1 1 . Percepciones, numero y porcentaje de especies de aves con percepciones y porcentaje del total de aves 
inventar iadas, Santo Domingo, 2007. 



Percepciones 


Numero de 
especies 


% de las especies con % del total de 
uso/percepcion (S=24) especies (S=49) 


Cultural 


7 


29,17 


4,08 


Beneticiosas para el ganado 


3 


12,5 


6,12 


Rapaces 


4 


16,67 


8,16 


Entermedades al ganado 


3 


12,5 


6,12 


Plagas 


10 


41,67 


20,41 



El Tirano Enmascarado o Dominico (Fluvicola nengeta), fue uno de los mas comunes visitantes 
de las areas construidas y frecuentemente habitadas por el hombre dentro de las tineas productivas; por 
esta razon, las personas de las relacionadas con las tineas se ban acostumbrado a su presencia, por su 
atractivo y por la talta de temor trente a la presencia del hombre, tanto ast, que se lo encuentra en 
parejas o solitario alimentandose en los pisos de dichas construcciones. 



Tabla 12. Aves con percepciones de tinqueros. Santo Domingo, 2009. 



Nombre Vulgar 


Especie 


Percepciones 


Aguila 


Accipiter superciliosus 


Rapaces 


Tangara azulejada 


Thraupis episcopus 


Plaga 


Espiguero Variable 


Sporophila corvina 


Plaga 


Eutonia coroninaranja 


Euphonia saturata 


Plaga 


Gallinazo 


Sarcoramphus papa 


Perjudicial para el ganado 
Beneticioso y perjudicial para el 


Garceta Bueyera 


Bubulcus ibis 


ganado 


Garrapatero 


Crotophaga suicirostris 


Beneticioso para el ganado 


Garza 


Butorides striatus 


Perjudicial para el ganado 


Gavilan Dorsigrts 


Leucopternis 

occidentaiis 


Rapaces 


Gavilan 


Leucopternis piumbea 


Rapaces 


Gavilan Gris 


Buteo nitidus 


Rapaces 


Golondrina 


Notiocheiidon 


Plaga 


Azuliblanza 


cyanoieuca 


Jilguero Ventriamarillo 


Cardueiis xanthogastra 


Plaga 


Loro alibronceado; 


Pionus chaicopterus 


Cultural 


Loro/perico 


Pionopssitta puichra 


Cultural 


Mosquitero verdoso 


Empidonax virescens 


Plaga 


Paloma 


Coiumba subvinacea 


Plaga 


Perdiz 


Geotrygon montana 


Plaga 


Ranchero 


Moiothrus bonariensis 


Beneticioso para el ganado 


Tangara tiliblanca 


Tachyphonus rufus 


Plaga 


Tortolita Ecuatoriana 


Coiumbina buckieyi 


Plaga 
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La importancia cultural de Plonus chalcopterus y Plonopssitta pulchra en la cosmovislon de los 
Tsachllas y clertos finqueros, fue un aspecto trascendental determlnado. Esto se puede deber a la gran 
abundancia y el constante contacto que la comunidad posee con ellas. Se conocio, que P. chalcopterus 
tiene un gran Importancia para los Tsachllas, en tal virtud, hasta se lo encuentra en uno de los 
monumentos mas Importantes del centre de la cludad (FIgura 7). 




Figura 7. Pionus chalcopterus en monumento a los Tsachllas (zona urbana de SDT). P. chalcopterus se observa 
en el antebrazo de un native (Foto: A. Eras, 2007). 

Metodo de cacen'a de las aves 

Los habitantes asentados y que concurren en las tineas productivas SDT, que han cazado aves 
silvestres, en su totalidad lo han realizado “sellandolos” (termino que los pobladores de SDT atribuyen al 
uso de escopetas, carabinas y armas de fuego en general, en actividades de caceria), sehalando que es 
lo mas practico si se busca alimento; sin embargo, se desconoce el metodo de caza para la especie 
Forpus coelestis, Plonopssitta pulchra y Pionus chalcopterus las mismas que se pueden destinar como 
mascotas. 

4.- DISCUSION 

Diversidad de aves 

Hasta ahora se desconoce estudios faunisticos en areas fragmentadas de SDT, por lo que este 
inventario permite conocer mas sobre el estado de conservacion de los remanentes del bosque muy 
humedo Premontano del occidente ecuatoriano (Anexo A). Este catalogo provee pautas para que 
parabiologos y cientificos hagan futures monitoreos en zonas fragmentadas, especialmente en SDT. 

Comparando este estudio con investigaciones hechas en areas protegidas, las cuales se han 
concentrado en varies localidades por periodos mas extensos y con la participacion de varies 
investigadores, la riqueza de STD (49 especies) es relativamente baja; En Sinangue (Reserve Ecologica 
Cayambe-Coca, Sucumbios) con identica zona de vida o formacion vegetal, se encontraron 278 
especies (Mena, 1997b); en Playa el Oro (Reserve Ecologica Cotacachi-Cayapas) se registraron 148 
especies (Benitez & Matheus, 1997), sin embargo, es logico que una zona alterada y fragmentada tenge 
menos diversidad de especies silvestres. 

En otros bosques riparios de fincas ganaderas como los de Matagalpa, Nicaragua, se 
encontraron 74 especies de aves (Perez et al., 2006), siendo mas ricos, comparadas con los bosques 
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riparios de fincas productivas en SDT. De igual manera, la diversidad se puede considerar baja con 
relacion a un estudio en Canas, Costa Rica, donde se registro 29 familias con 80 especies en la totalidad 
de un paisaje fragmentado (Cardenas et al., 2004), sin embargo, en la citada investigacion se 
caracterizaron las aves de varies habitats (fragmentos de bosque seco, bosques riparios, charrales, 
cercas vivas y potreros), siendo realmente la riqueza de los bosques riparios de 39 especies, lo que 
indica los bosques fragmentados del tropico ecuatoriano tienen aparentemente mayor biodiversidad, sin 
que esto signifique mayor o menor alteracion del area, lo que refleja la gran diversidad de aves presentes 
en Ecuador (Ridgely & Greendfield, 2006). 

La mayor o menor diversidad que presenta una tinea al memento de realizar los registros pudo 
estar influenciada por las condiciones climaticas, distancias recorridas, esfuerzo de monitoreo, 
facilidades del terrene y caracteristicas de los senderos, a pesar de que se realize las actividades de 
monitoreos lo mas homogeneamente posible, con el fin de que los dates representen exactamente la 
diversidad y abundancia de las mismas. Puesto, que las salidas de campo fueron planificadas 
previamente a su ejecucion, no se pudo escoger condiciones uniformes para todos los monitoreos, 
cumpliendo de este mode las fechas preestablecidas, siendo este uno de los factores que se infiere 
como el que mas pudo afectar la diferencia de registros. Es probable que en futures trabajos se muestre 
mayor diversidad basados en uniformizar los factores previamente sehalados. 

La precipitacion pluvial caracterfstica de la epoca de estudio, pudo haber afectando los 
resultados de diversidad, ya que existe una menor movilidad de animales frente la presencia de Iluvias. 
Las posibilidades de avistamiento de las aves, tambien pudo afectar negativamente los dates de 
diversidad y abundancia (Mena, 1997b). 

Los ordenes con mayor riqueza en los bosques riparios de las fincas productivas fueron 
Passeriformes (19), Apodiformes (5), los Psittaciformes (3) y Columbiformes (3); estos dates, son 
similares a los obtenidos, en estudio realizado en el bmhPM, con la unica diferencia de haberlos obtenido 
en una area protegida (Reserva Ecologica Cayambe-Coca), donde los ordenes mas abundantes fueron 
Passeriformes, Piciformes, Apodiformes y Psittaciformes (Mena, 1997b). Los ordenes que difieren entre 
estas dos investigaciones son Piciformes y Columbiformes. El orden Piciformes en el actual inventario 
ocupo uno de los ultimos lugares en diversidad de especies, lo que indica que puede ser uno de los mas 
afectados frente a la alteracion de los ecosistemas locales. El orden Columbiformes fue mas abundante 
frente a la citada publicacion, sin embargo, no se le puede atribuir ninguna relacion con la alteracion de 
los bosques natives, ya que su alta abundancia se debe principalmente a la preferencia que tienen sus 
especies por las zonas riberehas y los bosques riparios. 

Los Passeriformes, Apodiformes y Psittaciformes tuvieron la misma tendencia de abundancia en 
ambos estudios, por lo cual se supone que estos a pesar de enfrentar los mismos desaffos ambientales 
son los que mas se han adaptado a los cambios del medio, sin que esto signifique que dichas 
alteraciones no los hayan afectado, ya que, muestra clara de este trastorno, es la diferencia de riqueza 
que se presenta entre los paisajes agrfcolas-ganaderos (S= 49) de la misma zona de vida frente el area 
protegida (S=74)(Perez et al., 2006; Pozo & Eras, este estudio). 

No obstante, en estudios previos de avifauna del Ecuador, los ordenes Passeriformes y 
Apodiformes, son mas abundantes en el occidente ecuatoriano, independientemente de la zona de vida: 
en la Cordillera de Toisan, Imbabura (Bosque muy humedo Montano bajo) (Gregory & Vaca, 1994); en la 
Playa de Oro (Bosque muy humedo Tropical) y Rfo Negro Chico (Bosque humedo Tropical), Zona de 
Amortiguamiento de la Reserva Ecologica Cotacachi Cayapas, Esmeraldas (Benftez & Matheus, 
1997;Gregory, 1997); en San Venancio, Esmeraldas (Bosque humedo Montano bajo) (Gregory, 1997). 
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Esto no pernnite apreciar diferencias entre zonas alteradas y areas protegidas. Sin embargo, en las 
citadas publicaciones realizadas en areas protegidas, el orden Piciformes se present© en tercer lugar en 
abundancia, por lo tanto, se puede inferir que es un orden sensible frente a las alteraciones ambientales. 

Bubulcus ibis, Pteroglossus erythropygius, Claravis pretiosa y Thraupis episcopus se las observe 
en grupos, en tal virtud se las registro en numeros altos, por lo que se presentan en los inventarios con 
una alta frecuencia. El numero de individuos de estas especies fue mayor, a pesar del menor numero de 
observaciones realizadas. 

El encontrar nuevas distribuciones de especies de aves no es frecuente, sin embargo, no es 
insolito, ya que las distribuciones cambian con el tiempo (Ridgely & Greendfield, 2006). Es probable que 
P. malaris no fuera registrada con anterioridad, debido a que la mayoria de registros se ban realizado en 
areas protegidas, dejando de lado areas alteradas, siendo estos ultimos, sitios potenciales de 
conservacion y alojamiento de aves silvestres (Estrada et al., 1993; Pozo, 1994; Harvey et al., 2005; 
Pozo & Eras 2007). 

Leucopternis occidentalis y L. plumbea, se encuentran bajo amenazada en el Ecuador, debido a 
la presion de caceria, otras formas de persecucion directa por ser rapaces diurnas de aves domesticas. 
Ademas, al ser especies paraguas, son mas sensibles que otras especies, a la destruccion de su 
entorno, la perdida de habitat y a la caza excesiva (Pozo et al., 2006). 

De las 14 especies endemicas del Ecuador se inventario a Pteroglossus erythropygius, en ambos 
sistemas productivos con una categorfa de abundante, lo cual indica, que pese a sus poblaciones bajas y 
distribucion limitada, es una de las pocas especies endemicas del pais, que aparentemente no se 
encuentra bajo amenaza de conservacion (Ridgely & Greenfield, 2006); posiblemente, porque a traves 
del tiempo ha tolerado bien los disturbios de su ecosistema, lo que se confirma con su registro en ambos 
sistemas productivos. Los bosques riparios son indispensables para esta especie dentro de areas 
fragmentadas, ya que, si bien tolera la perdida de su habitat, no puede vivir sin areas boscosas (Ridgely 
& Greenfield, 2006). 

Las especies endemicas de SDT (centre endemico “Bajuras Tumbesinas”), pueden indicar la 
alteracion de la zona (Gregory, 1997). En tal virtud, se puede inferir que las fincas ganaderas fueron 
menos alteradas con relacion a las agrfcolas, ya que tuvieron mas especies endemicas (Columbina 
buckleyi, Leucopternis occidentales, Forpus coelestis y Pachyramphus spodiurus) frente a las fincas 
agrfcolas (L. occidentales, F. coelestis y Geothlypis auricularis). 

Para conocer, realmente cual es la alteracion que tienen los bosques riparios de SDT, como 
zona fragmentada, es necesario comparar la mayor dominancia de las especies de acuerdo a los valores 
de abundancia proporcional (Pi), con un inventario similar en bosques no alterados de la zona, 
determinando asf, la presencia de ciertas especies como indicadoras del ambiente (Benftez et al., 1997). 

Es la primera vez que en el Ecuador se registra a Columbia livia en los bosques riparios de las 
fincas productivas (Ridgely & Greenfield, 2006) lo que indica la proximidad, que las areas ocupadas por 
el hombre tienen con las fincas productivas y los bosques riparios, reflejando la magnitud de la 
fragmentacion del paisaje. 

Las curvas de dominancia-diversidad de aves son muy parecidas en ambos tipos de sistemas 
productivos; las mismas que presentan lo que Odum (1986) denomina curva tipo C, es decir un patron 
sigmoideo intermedio. Esto indica diferenciacion y superposicion de nichos; al ser un sistema abierto, las 
aves coexisten en condiciones de competencia parcial, mas no directa, puesto que sus adaptaciones 
promueven la diferenciacion de nichos, sin embargo, no se excluye habitats por competencia. La curva 
tipo C, es modelo mas adecuado para comunidades sin mayores alteraciones (Odum, 1986; Benftez et 
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al., 1997), por ende, se puede inferir que las comunidades biologicas de las fincas agncolas como 
ganaderas relativamente no se encuentran altamente perturbadas. Cuando las perturbaclones son 
sutlles, y tan solo reducen el acceso a los recursos allmentlclos o los sitlos para anidar, y dependlendo 
de las caracteristicas de las comunidades, la redlstrlbuclon de los organismos despues de la 
perturbaclon, no necesarlamente refleja un efecto negative (Rossetti & GIraudo, 2003), puesto que el 
nuevo patron Impuesto de distrlbuclon puede ser temporal. Despues de un tiempo los animales pueden 
volver a su distrlbuclon original para explotar el remanente de los recursos (Gill et al., 1996), a pesar de 
que sus poblaclones hayan decrecldo. 

El indIce de Simpson concede mayor Importancia a las especles comunes, Indicando de este 
mode la aglutinaclon de domlnancia de una o pocas especles, mientras que el indIce de Shannon 
concede mayor Importancia a las especles raras (Odum, 1986; Benitez & Matheus, 1997). Los indices de 
Shannon y Simpson, fueron significativamente Iguales para los dos tipos de fincas, reflejandose en que 
las especles raras o poco abundantes y las comunes tuvieran un patron de posiclon parecldo. 

Etnobiologi'a de las aves 

Bubulcus Ibis), fue una especle de ave con una percepcion muy polar por parte de los finqueros, 
ya que fue considerada como transmisora de enfermedades al ganado y como beneficlosa para el 
mismo, esta ultima, en virtud de que se allmenta de sus Insectos perjudiclales. La realldad es que B. Ibis 
se allmenta de Insectos espantados por las vacas, slendo realmente plagas de plantaclones de banano, 
ya que gustan de esta fruta (RIdgely & Greenfield, 2006). 

Pteroglossus erythropyglus), es una ave que se allmenta de frutas, lagartijas, huevos de aves y 
hasta crias, slendo esta una rapaz mas de SDT (RIdgely & Greenfield, 2006), mas no es perciblda de 
este mode por los finqueros. 

NInguna ave ha side cazada dentro de las fincas anallzadas. Esto se debe, a que cualquier 
mamffero susceptible de caza, proporclona mas came por bala que una ave, por lo que son preferldos 
por los cazadores (Gregory, 1 997). La no cacerfa de aves, muestra que la preslon de caceria a diferencia 
de la perdida de habitat, no es un problema que afecte la diversidad ornitologica de SDT, a pesar de que 
clertos finqueros perciben notorlamente que los Columbldae son susceptibles de ser cazadas como 
allmento. Probablemente no se cazan aves en este sector, porque las pavas grandes (Cracldae), en el 
occldente ecuatorlano nunca han sido abundantes, conoclendo solo la existencia de Penelope 
purpurascens en dicha area del Ecuador (RIdgely & Greenfield, 2006). 

Los bosques riparlos dentro de las fincas productivas, presentaron relativamente poca 
diversidad, pero excelentes Indiclos de conservaclon, que en teorfa no se esperaban, pues se encontro a 
una de las 14 especles de aves endemicas del Ecuador, en los dos sistemas productivos estudlados 
(Pteroglossus erythropyglus); al Igual que Columblna buckleyl, Forpus coelestls, Leucopternis 
occldentales y Pachyramphus spodlurus, como Indicadoras de buen amblente; como bandera a P. 
spodlurus, Phaetornis malarls, Plonopsitta pulcra, Plonus chalcopterus y Ramphastos brevis tanto como 
bandera y paraguas a Accipiter superclllosus, Buteo nitidus, L. occldentales, L. plumbea y P. 
erythropyglus (Gregory, 1997; Mena & Ruiz, 1997; Mena & Utreras, 1997; Mena, 1997; ANCON, 2000; 
Guerrero et al., 2006; Del Moral, 2007; UlCN, 2006). 

La conservaclon de una comunidad blologica, depende de la estructura poblaclonal y especifica 
de las especles que la conforman; el encontrar las sehaladas especles "focales" (paraguas, clave, 
Indicadoras y bandera), en vista de su caracter articulador de la comunidad blologica, su gran 
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importancia dentro del ecosistema (al mantener el equillbrio de un ampllo range de plantas y anlmales), 
su atractivo frente a la colectividad de seres humanos y el alto grade de amenaza que enfrentan, pernnite 
sostener que las alteraclones de la comunldad blologica de los basques riparlos, no ban side de gran 
magnltud. 

Esto brinda buenas perspectivas de conservaclon, ya que emprendlendo un arduo trabajo con 
los finqueros, se aumentaria y mantendria los remanentes de bosques natives dentro de sus tineas, con 
lo que se puede ampllar las subpoblaclones de especles sllvestres; que en base del concepto de 
metapoblaclon (Suarez, 1998), dichas poblaclones pequenas (de cada especle), podrian ser manejadas 
como una sola poblaclon, cuyo tamano colectivo puede ser suflclentemente grande como para pernnitir 
su persistencia a largo plazo (HanskI & Gilpin, 1991; Suarez, 1998). 

Esta propuesta de manejo, bajo el slstema de “metapoblaclon” (HanskI & Gilpin, 1991; Suarez, 
1998), propone, Incrementar y mantener los remanentes boscosos entre tineas adyacentes, consecutivas 
y cercanas, para que las poblaclones pequenas locales (subpoblaclones) de estos bosques remanentes, 
Interactuen a traves de los Individuos que se movlllzan entre dichas poblaclones, lo que Implica que 
todas estas subpoblaclones sean manejadas como una sola poblaclon. 

Para aumentar la conservaclon de especles sllvestres de los bosques riparlos, trente a la real 
alteraclon que sutre el palsaje, se debe tomentar esquemas de protecclon y restauraclon de los 
corredores riparlos (todas aquellas areas boscosas tluviales conectadas, que disminuyen el alslamlento 
entre Islas habitats) mantenlendo la conectividad blologica entre las tineas, en las cuales se promueva en 
los productores la regeneraclon e Incremento del area torestal (Gordlllo et al., 2006), creando “habitats 
temporales” y “permanentes” para estas especles (Burel 1996; Harvey et al., 2005). 

No obstante, el creer clegamente en uno de los dos polos que han caracterizado a la comunldad 
agroblologica, es declr, o blen la producclon de allmentos sin que Importe la ecologia, o blen la 
conservaclon de la naturaleza sin que Importe la allmentaclon humana, tomenta el actual esquema 
Inmanejable de producclon, donde la posible y blen tundamentada sostenibllldad no es tomada en 
cuenta, slendo la talta de consensos, el primer problema causante de la Imparable tragmentaclon y 
perdida de bosques natives (Obs. pers.). Esta sostenibllldad, se entlende como la utlllzaclon de la 
diversidad blologica de un mode y a un ritmo que no ocaslone su disminuclon en el largo plazo, 
mantenlendo las posibllldades de satlstacer las necesidades y aspiraclones de las generaclones actuales 
y tuturas (Mena et al., 1 997). 

Una alternativa trente al estllo de desarrollo -Irraclonal y depredador-, como el Implementado en 
la Amazonia y Costa ecuatorlana causa Irreversible de la tragmentaclon de los bosques natives (Suarez, 
1998), es el tomento del manejo de vIda sllvestres en cautiverlo y semicautiverlo, en areas sin alteraclon, 
en las cuales con una considerable cantidad de terrene se puede produclr de manera extensiva, 
explotando al menclonado recurso al ritmo de su reproducclon natural (Tapia, 1998). El agroturlsmo es 
otra alternativa para mejorar la calldad de vIda de las personas principalmente rurales (Pozo & Cardenas, 
2008). El agroecoturlsmo puede nacer en SDT, en busqueda de opciones para una mayor productividad 
de las tierras, con clerto caracter turistico, en vista de la gran cantidad de movimlentos “ecoturisticos” 
que giran a sus alrededores (Arrieta, 2000). 

Con el agroturlsmo se mantlene los habitats naturales y la comunldad blologica que la reside, en 
perspective de los grandes beneticlos economicos y soclales derivados de su uso sostenible (Pozo & 
Cardenas, 2008). Esta alternativa, comprende el turlsmo en areas agropecuarlas, en las cuales se 
aprovecha el medio amblente rural, ocupado por una socledad autoctona, que muestra y comparte no 
solo su Idlosincrasla y tecnicas agricolas, sino tamblen su entorno natural protegido, sin embargo, solo 
se puede promover en base de la protecclon y uso sostenible de los recursos, con tundamentos de 
consclencia ecologica en la Idlosincrasla de las personas (Arrieta, 2000) 
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La inexistencia de un mercado y forma de comercializacion de aves, sugiere que no se las altera 
Incldentalmente. En tal virtud, una actividad Importante y positiva apllcable a SDT, que se esta 
desarrollando en el Ecuador desde hace unas pocas decadas es el ecoturlsmo y el agroecoturlsmo, 
enfocada principalmente en las aves, pues es el grupo de entre los organismos vivos mas popular 
globalmente (RIdgely & Greenfield, 2006), slendo muy Importante en el area continental y mas aun en 
Galapagos; por esto, el aviturlsmo es una de las alternativas que Impactaran positivamente los 
ecosistemas no solo de areas protegidas sino de palsajes fragmentados, como el de SDT 

5. - CONCLUSIONES 

Se reconocleron 36 especles de aves (21 en fincas agricolas y 25 fincas en ganaderas), con 199 
Individuos (26% en agricolas y 74% en ganaderas), perteneclentes a 8 ordenes y 16 famlllas. Los 
ordenes con mayor rlqueza fueron Passeriformes, Apodiformes, Psittaciformes y Columbiformes. Las 
especles mas abundante fueron: el Espiguero Variable (Sporophlla corvina) (12.6%) en agrfcolas y, 
Garza (Bubulcus Ibis) (15.1%) en ganaderas. Phaethornis malarls (Collbrf ermitano) habita en el 
occldente ecuatorlano como nuevo registro. El 47% de aves dependen, facultativamente o estrictamente 
de bosques y arboledos, mientras que el 53% no dependen en absolute de bosque. 

Los promedlos de rlqueza (S), diversidad (H') y la abundancia (N), al Igual que los indices 
ecologicos de diversidad Simpson, su reciproco, equitatividad y Margalef, no presentaron diferenclas 
estadisticas entre los bosques riparlos de los sistemas productivos propuestos. Esto muestra Iguales 
niveles de perturbaclon ecologica para las comunidades de especles que habitan en los bosques. 

Existio un evidente y progresivo aumento de especles conforms aumentaban los Individuos 
Inventarlados. Sin embargo, se observe y en base de las curvas de acumulaclon de especles, que 
faltaron especles por ser registradas, y por ende campanas de muestreo por sistema productive. 

El 42% de finqueros le dieron uso y/o percepcion a alguna o varlas especles de aves. El 49 % de 
las aves, presentaron uno o varies usos/percepciones. NInguna especle de ave es comerclallzada por los 
propletarlos de las fincas productivas. Los uses de las aves se restringleron al de allmento (14.3%) y el 
de ornamento o mascotas (6.1 1%). El 4.1% de aves fueron percibldas con Importancia cultural, 6.1% con 
beneficlos para el ganado e Igual porcentaje como perjudiclales para el ganado. El 20.4% de aves eran 
percibldas como plagas de cultivos y el 8.2% como rapaces. 

El mantener los remanentes de bosques riparlos, dentro de los palsajes productivos, ha 
significado un gran aporte para la conservaclon de especles sllvestres en SDT, frente a la perdida de sus 
habitats, tanto asi, entre aves y mamfferos (Pozo & Eras, en prep.) que se registraron 2 especles 
“Vulnerables”, 4 “CasI amenazadas”, 1 “En Pellgro”, 1 “Al Margen de RIesgo” y 13 especles dentro de los 
apendices del CITES. Las perturbaclones, existentes en la comunidad blologica de los bosques riparlos, 
no han sido de gran magnitud, en vista de los registros de 1 ave endemica del Ecuador, 8 especles 
Indicadoras de un buen ecosistema, 12 especles paraguas y 17 especles bandera. 

6. - RECOMENDACIONES 

Realizar estudlos de regeneraclon y aumento de bosques riparlos y areas boscosas dentro de las 
fincas productivas, para reevaluar la diversidad y abundancia de especles, y las metapoblaclones fruto 
de estos Incrementos. 
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Realizar estudios similares entre niveles de intensificacion de sistemas productivos, y en otros 
habitats de las tineas come: cercas vivas, setos vivos, arboles disperses en potreros. 

Aumentar los estudios interdisciplinarios entre agropecuarios, biologos y parabiologos, para 
mejorar y/o crear un esquema mas equilibrado de produccion, basado en la sostenibilidad (utilizacion de 
la diversidad biologica de un mode y a un ritmo que no ocasione su disminucion en el largo plazo). 

Los registros de aves podrian aumentar, y ampliarse a las demas clases de vertebrados, si se 
amplian las campahas de muestreo, los dias por campaha, y el personal encargado. 

Los valores de diversidad y abundancia de especies silvestres se pueden incrementar, y 
representar un ejempio para las demas zonas productivas ecuatorianas, en base de aumentar y 
mantener los remanentes de bosques natives dentro de las tineas productivas, con lo que las 
poblaciones podrian tener un tamaho colectivo suticientemente grande como para persistir 
indetinidamente. Creando de este mode corredores riparios, los que mantendran la conectividad 
biologica entre las tineas productivas, tundamentando su sostenibilidad y sustentabilidad en: 1) el 
agroturismo, aprovechando la gran cantidad de movimientos “ecoturisticos” de SDT y, 2) el aviturismo, 
actividad altamente impulsada en las ultimas decadas en el Ecuador, en virtud de que las aves son los 
organismos vivos mas populares. 
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Anexo A. Lista de aves de los bosques riparios de fincas agricolas y ganaderas, Santo Domingo, 2009. 



Especie 


Nombre Vulgar 


Fincas 


Registro 


Agri. 


Gan. 


APODIFORMES 


Trochilidae 


Amazilia amabilis* 


Amazilia Pechiazul 




P 


D 


Amazilia tzacatl* 


Amazilia Colirrufa 




P 


D 


Androdon aequatorialis* 


Colibri piquidentado 






I 


Damophila julie* 


Colibri ventrivioleta 




P 


D 


Phaetornis malaris* 


Ermitano Piquigrande 


P 


P 


D 


Phaetornis yaruqui* 


Ermitano Bigotiblanco 




P 


D 


CICONIIFORMES 


Ardeidae 


Bubulcus ibis 


Garceta Bueyera 




A 


D 


Butorides striatus 


Garza 






I 


ICathartidae I 


I Sarcoramphus papa* 


Gallinazo 






I I 


COLUMBIFORMES 


Columbidae 


Ciaravis pretiosa* 


Tortolas;Tortolita Azul 




A 


D 


Coiumba subvinacea* 


Paloma 






I 


Coiumbia iivia 


Paloma 


P 




D 


Coiumbina buckieyi 


Tortola; Tortolita Ecuatoriana 




P 


D 


CUCULIFORMES 










Cuculidae 










Crotophaga suicirostris 


Garrapatero 




P 


D 


FALCONIFORMES 


Accipitridae 


Accipiter superciiiosus* 


Aguila 






I 


Buteo nitidus 


Gavilan Gris 


M 




D 


Leucopternis occidentaiis 


Gavilan Dorsigris 


P 


P 


D 


Leucopternis piumbea 


Gavilan 






I 


PASSERIFORMES 


Icteridae 


Cacicus ceia 


Cacique 






I 


Moiothrus bonariensis 


Rachero 






I 


Emberizidae 


Arremon aurantiirostris* 


Salton Piquinaranja 


P 




D 


Oryzoborus angoiensis* 


Negro yermo; Semillero menor 


P 


A 


D 


Sporophiia corvina 


Espiguero Variable 


A 


A 


D 


Fringillidae 


Cardueiis xanthogastra 


Jilguero Ventriamarillo 


P 




D 


Furnariidae 


Synaiiaxis brachyura* 


Carpintero; Colaespina 

pizarrosa 


A 


P 


D 
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Continua en la siguiente paginal 

* Especies que dependen estrictamente o facultativamente de bosques, Agri = Fincas 
Agncolas, Gan = Fincas Ganaderas, A = Abundante, M = Medianamente Abundante, P= 
Poco Abundante, D = Directo (Capture, Observacion, Rastro), I = Indirecto (Encuestas) 
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Lista de aves de los bosques riparios de tineas agricolas y ganaderas, Santo Domingo, 2009. 
(Conti nuacion) 



Especie 


Nombre Vulgar 


Fincas 


Registro 


Agri. 


Gan. 


PASSERIFORMES 


Flirundinidae 


Neochelidon tibialis minima* 


Golondrina 




M 


D 


Notiochelidon cyanoleuca 


Golondrina Azuliblanca 




P 


D 


Progne chalybea 


Martin 






I 


Riparia riparia 


Golondrina riberena 




P 


D 


Mimidae 


Mimus longicaudatus 


Cucube 


P 




D 


Parulidae 


Geothlypis auricularis* 


Antitacito lorinegro 


M 




D 


Myioborus miniatus* 


Candelita goliplomiza 


P 




D 


Thraupidae 


Euphonia saturnata 


Eutonia coroninaranja 


M 




D 


Ramphocelus icteronotus* 


Tangara Amarilla o lolilimon 




M 


D 


Tachyphonus rufus 


Tangara tiliblanca 




P 


D 


Thraupis episcopus 


Azulejo; Tangara azulejada 


P 


A 


D 


Tyrannidae 


Colonia colonus* 


Viuditas 






I 


Empidonax virescens 


Amarillo; Mosquitero Verdoso 


M 


M 


D 


Fluvicola nengeta 


Dominico; Tirano Enmascarado 


A 


A 


D 


Pachyramphus spodiurus 


Cabezon Pizarroso 




P 


D 


Tityra inquisidor* 


Titira coroninegra 


P 




D 


PICIFORMES 


Picidae 


Ceieus loricatus* 


Carpintero Canelo 


A 




D 


Piculus rubiginosus* 


Carpintero 






I 


Ramphastidae 


Pteroglossus erythropygius* 


Dios te de. Arasari Piquipalido 


A 


A 


D 


Ramphastos brevis* 


Tucan 






I 


PSITTACIFORMES 


Psittacidae 


Forpus coelestis 


Periquito del pacitico 


M 


M 


D 


Geotrygon montana* 


Perdiz 






I 


Pionopsitta pulchra 


Loro 0 perico 






I 


Pionus chalcopterus 


Loro alibronceado 




P 


D 



* Especies que dependen estrictamente o tacultativamente de bosques, Agri = Fincas 
Agricolas, Gan = Fincas Ganaderas, A = Abundante, M = Medianamente Abundante, P= 
Poco Abundante, D = Directo (Capture, Observacion, Rastro), I = Indirecto (Encuestas) 
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RESUMEN 



Este articulo discute el desarrollo de modelos analiticos, usando la tecnica de secciones 
divididas en fibras, para describir el comportamiento estructural no lineal de conexiones de hormigon 
pila-losa de cubierta sujetas a cargas ciclicas laterales. Los modelos son validados en parte usando 
resultados experimentales obtenidos de la literatura tecnica sobre el tema. En el nivel de seccion, las 
comparaciones estan enfocadas principalmente a los comportamientos esfuerzo-deformacion y 
momento-curvatura, mientras que al nivel de elemento estructural la comparacion de memento vs. 
desplazamiento lateral es usada para arrojar luz sobre el comportamiento global de la conexion pila-losa 
de cubierta. Los resultados indican que los modelos desarrollados son capaces de predecir con precision 
los comportamientos experimentales local y global, aun incluyendo la ocurrencia de varies estados de 
daho. Finalmente, se conduce un estudio para explorar el impacto del uso de chaquetas de polimero 
reforzado con fibras de carbono (CFRP), en la zona de la rotula plastica de la pila, en el comportamiento 
dclico de la conexion. La aplicacion de chaquetas CFRP puede retrasar la separacion del recubrimiento 
de hormigon y per tanto mejorar la resistencia y la ductilidad de la conexion. 



ABSTRACT 



This paper discusses the development of analytical models, using the fiber section 
technique, to describe the nonlinear structural behavior of concrete pile-wharf connections subjected to 
cyclic lateral loading. The models are validated in part using experimental results obtained from the 
literature. At the section level, comparisons are primarily focused on the stress-strain and moment- 
curvature behaviors, while at the element level comparisons of lateral force and moment vs. displacement 
are utilized to shed light on the global behavior of a pile-wharf connection. Results indicate that the 
developed models are able to accurately predict the global and local experimental behavior, even 
including the occurrence of various damage states. Finally, a study is conducted to explore the impact of 
using Carbon Fiber 
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Reinforced Polymer (CFRP) jackets at the plastic hinge zone of the pile on the cyclic behavior of 
the connection. The application of CFRP jackets can delay spalling of the cover concrete and thus 
enhances the strength and ductility of the connection. 

1 . INTRODUCCION 

La respuesta sfsmica de estructuras portuarias tales come muelles apoyados sobre pilas puede 
asumirse que esta fuertemente relacionada a un par de factores: la configuracion geometrica de la 
estructura y las condiciones del suelo de cimentacion. Dos zonas crfticas son 

reconocidas en las pilas - una dentro del suelo de soporte y otra en la conexion pila-losa de 
cubierta. Esta ultima zona deberia merecer una consideracion especial porque es usualmente accesible 
y por lo tanto podrfa ser relativamente facil de reforzar. Debido al incremento del trafico comercial entre 
el continente americano en general y Asia, hay una considerable demanda para expandir puertos y 
construir otros nuevos. Por lo tanto, una especial atencion es requerida en la planificacion y construccion 
de nuevas instalaciones portuarias, tal como podrfa ser el Puerto de Manta, en la figura 1 (a) 

Entre las conexiones pila-losa de cubierta, un subtipo comun comprende una pila de hormigon 
presforzado con una losa de cubierta fundida in-situ. La pila y la losa de cubierta son conectadas por 
medio de barras “pasador” (dowel bars), con la pila embebida de 5 a 7.5 cms. dentro de la losa, como se 
puede ver en la Figura 1(b). El trabajo pionero sobre esta materia fue realizado en la Universidad de 
Canterbury, Nueva Zelanda (Joen et al. 1988). Su especimen (40 cms. de seccion octogonal con cuatro 
D20 barras pasador) fue considerado el de peer comportamiento, porque la rotacion plastica concentre el 
daho de forma no deseada en una amplia grieta cerca de la interface pila-cubierta de pila, pero los 
investigadores sugirieron que este comportamiento podrfa mejorar con barras pasador adicionales. Un 
estudio de 1997 conducido en la Universidad de California San Diego (Silva et al. 1997) ilustro que una 
conexion pila-cubierta a escala natural (pila tipo Caltrans Class 70 ton. con 6-#6 barras pasador 
embebidas en la cubierta) puede ser susceptible de una reduccion significativa en su capacidad a flexion 
debido a una gran extension de la perdida de recubrimiento de hormigon cuando es sometida a cargas 
laterales cfclicas junto con carga axial variable. Otra conexion pila-cubierta que usaba barras pasador 
longitudinales y barras con cabeza en forma de T actuando como enlace en la region nodal fue luego 
probada bajo carga lateral cfclica (sin carga axial) para el Puerto de Los Angeles (POLA) (Sritharan and 
Priestley 1998). Elios concluyeron que este tipo de barras aseguraban la ductilidad de la conexion. 
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Figura 1.(a) Puerto en Manta y (b) Esquema de la conexion pila-losa de cubierta 



Mas recientemente, un estudio conducido en la Universidad de Washington (Roeder et 
al. 2001) involucrando el ensayo de algunas conexiones pila-cubierta, con pilas apenas 
embebidas dentro de la cubierta, indico que definitivamente tales conexiones podrian sostener 
daho significative bajo cargas laterales reversibles. En un estudio experimental y analftico 
patrocinado por POLA, dos conexiones pila-cubierta a escala completa (60 cms. de seccion 
octogonal) fueron probadas bajo carga lateral ciclica (Restrepo et al. 2007). Un ensayo 
represento la conexion “no-sismica” pila-cubierta (4#9 barras pasador embebidas 45 cms. 
dentro de la cubierta), mientras que el segundo ensayo represento una conexion “sismica” (8- 
#10 barras pasador embebidas 75 cms. dentro de la cubierta). Ambos ensayos mostraron que 
este tipo de conexiones pueden tener respuestas predecibles y que los desplazamientos 
laterales correspondientes a los limites de deformacion requeridos por el Codigo para Diseho 
Sismico, Modernizacion y Reparacion de Muelles para Contenedores del Puerto de Los 
Angeles (Puerto de Los Angeles 2004) pueden ser predichos con precision. Finalmente, cuatro 
especimenes a escala completa fueron recientemente probados en la Universidad de 
Washington para evaluar y comparar el rendimiento de diferentes detalles de conexion (Jellin 
2008). Estos detalles que incluyeron la separacion parcial de las barras pasador del hormigon 
circundante, la colocacion de un cojfn semi-rfgido entre la pila y la cubierta, y la adicion de un 
abrigo de espuma sintetica en el perfmetro de la corta longitud embebida de la pila parecieron 
reducir la extension del daho en la interfase pila-cubierta (a pesar de alguna perdida en la 
rigidez de la conexion). 

Todos los especfmenes descritos mas arriba experimentaron daho en forma de 
separacion de los recubrimientos de hormigon tanto en la pila como en la cubierta y la 
exposicion de los refuerzos longitudinal y lateral a “derivas” entre el 2% y el 5%, aun si la 
capacidad flexural fue aceptable. Ademas, a pesar de todos los estudios experimentales 
descritos previamente, aun hay una considerable falta de herramientas analfticas y modelos 
que podrfan ser usados para investigar y predecir el comportamiento de conexiones pila- 
cubierta bajo carga sfsmica. Este artfculo se enfocara por tanto en el desarrollo y calibracion de 
modelos analfticos precisos para conexiones de hormigon armado pila-cubierta que puedan ser 
usados en futures investigaciones. Los modelos desarrollados podran tambien ser utilizados 
para explorer el comportamiento de conexiones cuando son reforzadas con chaquetas de 
polfmeros reforzadas con fibras de carbono (CFRP). 

2 . MODELO ANAUTICO - DESCRIPCION GENERAL 



Los componentes fundamentales del modelo analftico desarrollado son la pila, la 
cubierta y su interfase. La pila es modelada en una posicion sobre la cubierta porque este es el 
arreglo experimental tfpico. Sus diferentes secciones, si bien siempre octogonales, modelan la 
variable presencia y caracterfsticas del hormigon (incluyendo los efectos de confinamiento), del 
acero presforzado, y de las barras pasador, como se puede ver esquematicamente en la Figura 
2. La cubierta propiamente dicha es globalmente considerada como un cuerpo rfgido, pero su 
influencia sobre el comportamiento de la pila sf se modela gracias a las caracterfsticas de la 
interfase pila-cubierta. Debido al hecho de que los elementos de enlace son las barras pasador, 
su efecto en el comportamiento plastico de la interfase (principalmente debido a la penetracion 
de la cedencia del acero) es tornado en cuenta cuidadosamente. El software usado es 
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OpenSees (Mazzoni et al. 2009) debido a su conveniente implementacion de la tecnica de 
secciones analizadas per medio de fibres y a su biblioteca de materiales que incluye un amplio 
rango de modelos no lineales de materiales que ban side desarrollados especificamente para 
aplicaciones sismicas. 



Axial load 




5 ^ 



Figura 2. Caracterfsticas generates del modelo y vista de una seccion tfpica de pila 

3. MODELO DE PILA 

La tecnica de modelacion de secciones por medio de fibres, que es a la fecha el mas 
eficiente mode de analizar elementos de hormigon presforzado con comportamiento plastico, 
es usada para modeler las secciones de la pila. Un papel clave dentro de este tipo de analisis 
es el del modelamiento de los materiales. Cuatro modelos basicos son requeridos: hormigon no 
confinado para el recubrimiento de la pila, hormigon confinado para el nucleo, acero 
presforzado para los tendones, y acero comun para las barras pasador que conectan la pila a 
la cubierta. 

Modelamiento del Hormigon 

El modelo para hormigon ConcreteOf de OpenSees fue usado para representar el 
comportamiento esfuerzo vs. deformacion del hormigon; esta basado en el modelo uniaxial de 
Kent-Scott_Park con rigidez de descarga/recarga linealmente degradada de acuerdo al trabajo 
de Karsan y Jirsa (1969). En este modelo, la resistencia a tension del hormigon no es 
considerada. Por otro lado el comportamiento a compresion es definido por los valores de 
esfuerzo y deformacion en los puntos denominados en este artfculo “pico” y “ultimo”. La curva 
hasta la resistencia pico es parabolica y luego lineal. Las curvas para hormigon no confinado y 
hormigon confinado (en la zona de la rotula plastica) usadas en este estudio pueden ser vistas 
en Fig. 3, asf como una curva para el hormigon reforzado con CFRP (la que sera explicada 
mas tarde). En el hormigon no confinado la deformacion a la resistencia pico es calculada con 
las ecuaciones de Thorenfeldt et al. (1987), porque cubren un amplio rango de resistencias del 
hormigon hasta valores de 18 ksi. La resistencia del hormigon se asume como cero en el punto 
de deformacion ultima, el cual fue de nuevo calculado usando las ecuaciones de Thorenfeldt et 
al. (1987). 

En la vecindad de la conexion pila-cubierta, el nucleo de la pila esta tfpicamente 
confinado con una relativamente alta cantidad de refuerzo transversal en forma de espirales. 
Considerar el efecto de tal confinamiento en el comportamiento constitutivo del hormigon en el 
nucleo es crucial. Los datos de entrada para el modelo de material ConcreteOI fueron por tanto 
modificados para considerar el efecto del confinamiento. Esfuerzo y deformacion en el punto 
pico se determinaron basandose en el modelo presentado por Mander et al. (1988), pero el 
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esfuerzo residual del hormigon (20% del esfuerzo pico) y deformacion (aproximadamente 10 
veces la deformacion en el esfuerzo pico) fueron considerados tomando en cuenta valores 
tfpicos para un modo de falla de deslizamiento friccional ultimo (Mazzoni et al. 2009). 




Strain [in/in] 

Figura 3. Modelos de hormigon no confinado, confinado y reforzado con CFRP 

Modelamiento del Acero 

Como se ilustra en la Figura 2, la pila esta tfpicamente reforzada con acero regular (en 
la forma de barras pasador que cruzan la interface pila-cubierta) y acero presforzado. El 
modelo uniaxial de material Steel02 disponible en la librerfa de materiales del programa 
OpenSees fue usado para describir el comportamiento de ambos tipos de refuerzo. En este, el 
comportamiento esfuerzo-deformacion del acero de refuerzo es descrito por el modelo no lineal 
de Menegotto y Pinto (1973). Para elementos tales como pilas, los cuales tienen niveles de 
compresion axial bajos a moderados, la curva monotonica esfuerzo-deformacion provee una 
respuesta razonable en el range a tension pero no en el de compresion (Priestley et al. 1996). 
Afortunadamente la contribucion del acero a compresion no es tan importante como la del 
hormigon, por lo que el modelo de Menegotto-Pinto es usado en su forma estandar. Se 
considera, sin embargo, que el acero presforzado tiene un efecto mfnimo sobre el 
comportamiento de la pila cerca de la conexion, principalmente debido a la limitada 
transferencia de esfuerzos entre el acero presforzado y el hormigon circundante en esa region. 
Sin embargo vale la pena decir que el acero presforzado juega un importante papel en la 
reduccion del agrietamiento a lo largo de la longitud de la pila. 

4. MODELO DE LA INTERFASE PILA-CUBIERTA 

El cuidadoso modelamiento de la rotacion inelastica en el extreme de la pila es de la 
mas extrema importancia porque esta asociada con la extension del daho estructural y con las 
estrategias de refuerzo. Ignorar esto puede sobrestimar la contribucion de los mecanismos de 
flexion regulares y por lo tanto sobrestimar el daho estructural (Zhao y Sritharan 2007). La 
rotacion inelastica del extreme de la pila es debida a dos fenomenos: plastificacion del 
hormigon y deslizamiento del refuerzo. Debido a la coexistencia de estos dos fenomenos, la 
rotacion inelastica del extreme de la pila sera modelada con un resorte rotacional con dos 
componentes (acero y hormigon) que los caractericen. Este resorte usara un elemento 
analizado con fibras pero de longitud unitaria (zero-length section element), que por su 
generalidad permite la captura del efecto de penetracion de deformacion sin importar la forma 
de la seccion transversal y la direccion de la carga lateral (Zhao and Sritharan 2007). 

Resorte Rotacional 

La rotacion en una grieta por flexion es igual al deslizamiento de la varilla de acero 
dividido por la distancia entre esta y la punta de la grieta. Esta altura depende de las 
caracterfsticas del hormigon plastico (con resistencia reducida) (Oehlers et al. 2009). El modelo 
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para el hormigon presentado per Oehlers et al. (2009) fue considerado un metodo conveniente 
para replicar el comportamlento plastico del hormigon en la Interfase. Como se llustra en la Fig. 
4(a), bajo carga lateral excesiva y debldo a la preslon en la base de la plla, una cuha a 
compreslon se desarrolla. Como se llustra en la Fig. 4(a), la fuerzas paralela y perpendicular a 
la cara Incllnada (el piano de cortante) pueden ser calculadas por equillbrio y luego 
relaclonadas usando la clasica teoria de Mohr-Coulomb para pianos de falla por cortante. La 
resistencia a compreslon longitudinal efectiva de la zona de resistencia reduclda del hormigon 

al comlenzo de la falla por cortante fricclonal es 

_ c + cos a(sin a + m cos a) 
sin a(cos a - m sin a) 



( 1 ) 



Donde: c es el termino cohesive, es la preslon de confinamlento, m es la pendlente 
del piano de falla en la ecuaclon de Mohr-Coulomb, y ^en el angulo entre la vertical y el piano 
Incllnado de la cuha de falla 



a = tan 




( 2 ) 



La rama Iniclal de la curva esfuerzo-deformaclon para el hormigon sigue la forma 
parabolica usual. Sus valores pIco son la resistencia a compreslon no confinada fco y la 
deformaclon a la resistencia no confinada pIco eco. SI hay un nuevo Incremento en 
deformaclon, aun a pesar de que la resistencia del hormigon no haya alcanzado su maximo 
confinado, las cuhas por cortante fricclonal ocurren. Notese que la profundidad de la zona de 
resistencia reduclda (plastica) del hormigon puede aumentar hasta un Ifmite que depende del 
desllzamlento total de la cuha. 

La metodologfa descrita mas arriba fue Implementada a traves del modelo ConcreteOt 
y puede ser vista en la Fig. 4(b). 



(b) 





Figura 4. (a) y (b) Hormigon plastico (Oehlers et al. 2009) y (c) esfuerzos de adherencia para barras de 
refuerzo (Haskett et al. 2009) 



El cuidadoso modelamiento del desllzamlento del acero es el factor mas importante en 
la descripcion de la penetracion de cedencia, como tambien es explicado por Zhao y Sritharan 
(2007). El modelo para acero presentado por Haskett et al. (2009) es usado en este artfculo por 
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su base matematica, en contraste con la usual experimental, respecto a la descripcion de la 
penetraclon de cedencia de las barras de refuerzo, por lo que puede ser facllmente adaptado 
para diferentes condiclones de adherencla. Este modelo establece que si las barras de refuerzo 
permanecen elasticas y no se fracturan antes de fallar por adherencla: 
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Donde: 
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r T 

es la fuerza en las barras pasador, ”“es la resistencia de adherencla maxima, 

db es el diametro de las barras pasador, es el desllzamlento despues del cual el esfuerzo 

de adherencla permanece en cero, es el desllzamlento de las barras pasador, y E es el 
modulo de elasticldad del acero. La Fig. 4(c) muestra dos curvas para diferentes condiclones 
de adherencla. Estas ecuaclones son luego modificadas para tomar en cuenta las reglones de 
comportamlento del acero elastica y de endurecimlento por deformaclon y se llmitan por la 
resistencia a la fractura de las barras pasador. Este Ifmite es a menudo usado para conexiones 
plla-cublerta (Zhao y Sritharan 2007). FInalmente estas ecuaclones pueden ser modeladas a 
traves del modelo de material Hysteretic del programa OpenSees por su simpllcldad respecto al 
modelamlento del comportamlento cfcllco y sobre todo porque en un elemento “zero-length 
section”, que como se ha Indicado supone una longitud unitarla, el “desplazamlento” del 
elemento (en este case, desllzamlento de la barra de acero), es Igual a la deformaclon de la 
secclon. 

Reforzamiento con CFRP 

Parte del objetivo de este trabajo fue explorar la posibllldad de usar chaquetas de 
CFRP para reduclr el daho sufrido por las pllas en la region adyacente a la conexion plla- 
cublerta cuando son sometidas a cargas cfcllcas. Con este fin, el modelo analftico propuesto 
por Salldl et al. (2005) fue adoptado y usado para predeclr las propledades del hormigon si es 
confinado con una chaqueta de CFRP. En este artfculo una chaqueta de CFRP (30 cms. de 
longitud) fue disehada para un esfuerzo de confinamlento de 300 psi y una deformaclon de 
diseho ultima de 0.004 en la direcclon de las fibras de carbono (Caltrans 1996). El modulo de 
elasticldad del material CFRP fue asumido como 29,200 ksl. El modelo da entre otros 
resultados un 0.5% de volumen de CFRP con respecto al del hormigon, en otras palabras un 
total de 14 capas de CFRP. 

Validacion del Modelo 

El modelo analftico presentado en este artfculo fue comprobado con los resultados 
experlmentales obtenidos en la Universidad de Washington en su Especimen 9 (Jellln 2008). 
La plla de hormigon presforzado tiene una secclon octogonal de 60 cms. de diametro y una 
longitud de 260 cms. desde la Interfase hasta el punto de apllcaclon de la carga lateral. La losa 
de cublerta de hormigon fundido en sitio fue simulada con un bloque rectangular macizo de 
hormigon: 235 cms. de longitud, 132 cms. de ancho y 74 cms. de alto. 




108 



Revista CIENCIA Vol. 12 N.- 2 de 2009 



Descripcion de los materiales usados en el Especimen 9 

El hormigon usado en las pllas tuvo 8000 psi de resistencia de diseno a compreslon. 
Las barras con cabeza en T fueron fijadas dentro de tubos Incluldos en la plla y fundidas dentro 
de la cublerta para crear una conexion resistente a momento. Estas barras fueron ocho #10, 
ASTM 706 (Fy = 70 ksl), de 193 cms. de longitud. La plla tamblen tenia refuerzo en espiral 
W11, con un espaclamlento centro a centra variable entre 2.5 y 7.5 cms. El presfuerzo fue 
alcanzado usando veinte y dos, 1.25 cms en diametro, 270 ksl, tendones de baja relajaclon, 
esforzados a 31 kIps/tendon. La plla fue sujeta a una carga axial de 450 kips y a una historla de 
carga lateral consistente en ciclos entre 0 y 9% de “deriva”. 

Graficos Momento-Curvatura 

En la secclon crftica (cerca del extreme de la plla, en la Intersecclon con la cublerta), 
las curvas analftica y experimental se muestran en la Fig. 5. 




Notese que los resultados Incluyendo el reforzamlento con CFRP tamblen ban side 
Incluldos en la figura y muestran que el confinamlento con CFRP produce una mejora de la 
resistencia del hormigon. Con el fin de simpllficar el anallsis, la secclon octogonal fue modelada 
como circular; ademas la resistencia a tension del hormigon, como ya se ha Indicado, no se 
considero. Estas simpllficaclones se cree que son las principales razones para las llgeras pero 
perceptibles diferenclas en resistencia y pendlente Iniclal entre los resultados experimental y 
analftico. Tamblen es Importante notar como se llustra en la Fig. 5, que dado que el 
recubrimlento de hormigon (7.5 cms.) representa una relativamente gran porclon de la secclon 
de la plla, juega un papel Importante en la definiclon de la capacldad a momento de la secclon. 
De aquf que, luego del desprendimlento del recubrimlento, tamblen una relativamente gran 
reducclon en la capacldad a momento fue observada. 

Graficos Desplazamiento Lateral vs. Momento en la Base de la Pila 

Las relaclones entre el desplazamiento lateral a nivel del punto de apllcaclon de la 
carga lateral y el momento resistente en la union con la losa de cublerta se muestran en la Fig. 
6 para dos modelos analfticos y se pueden comparer con los resultados experlmentales. El 
primer modelo analftico (“Analytical As-Bullt”) considero la conexion plla-cublerta como se 
describe anterlormente, Incluyendo el resorte rotaclonal (pero no el refuerzo de CFRP). El 
segundo modelo analftico (“Analytical Retrofitted”) muestra los efectos de Incluir una chaqueta 
de CFRP con una longitud de 30 cms. en el extreme de la plla en su union con la losa de 
cublerta. 
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Displacement [in] 

Figura 6. Graficos Desplazamiento Lateral vs. Memento. 

Los resultados del modelo “Analytical As-Built” estan en muy buen acuerdo con los 
resultados experimentales, excepto por la zona de deterioro extreme de resistencia. Per otro 
lado el confinamiento adicional con CFRP basicamente produce un moderado pero continue 
incremento en ia capacidad pico a memento en comparacion con ei modeio sin este refuerzo. 
Este comportamiento se extiende a casi ei range compieto de despiazamientos ciciicos. Debe 
notarse que ia chaqueta de CFRP se asumio que se extendia dentro de ia porcion embebida 
de ia piia. 

5. CONCLUSIONES 

Este articuio desarroiia un modeio anaiitico de fibras para expiicar ei comportamiento 
de conexiones piia-iosa en ias que ia piia esta embebida dentro de ia cubierta de 2.5 a 7.5 
cms. Materiaies taies come hormigon no confinado y confinado, asi come acero reguiar y 
presforzado, se presentan, en adicion a modeios para hormigon piastico y para ei 
desiizamiento de barras de acero (ios cuaies fueron necesarios para describir con precision ia 
rotacion en ia interfase piia-cubierta). Los resuitados anaiiticos fueron comparados con 
diferentes resuitados experimentaies, enfatizando ei comportamiento iocai y giobai a traves de 
curves momento vs. curvature y momento vs. despiazamiento, respectivamente; ia coincidencia 
entre eiios se considera exceiente. Y finaimente un modeio basico para refuerzo con CFRP fue 
introducido para expiorar ei efecto de mejorar ia resistencia dei hormigon en ia conexion piia- 
iosa de cubierta. 
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ANEXOS 
Programa 1 

# 

# Geometria general y creacion de archives de dates 

# 

# define GEOMETRY 

set LCel 1 20; # celumn length 

set Weight 450.; # 550 limit superstructure weight [kips] 

set PCel 45.; # 5% celumn capacity (deberia ser igual a Weight) 

# set Mass [expr $Weight/386.4j; # nedal mass 

# nedal ceerdinates: 

nede 10 0; # nede#, X, Y 

nede 2 0 7.5; 

nede 3 0 37.5; 

nede 4 0 59.0; 

nede 5 0 82.5; 

nede 6 0 103; 

nede 7 0 112.5; 

nede 8 0 120; 

nede 9 0 -3; 

nede 10 0-3;# prueba flexien 

nede 11 0 -3; # prueba flexien+cempresien 

nede 12 0-3; 

# Single peint censtraints - Beundary Cenditiens 
fix 12 10 0; 

fix 9 1 1 1 ; # nede DX DY RZ 

fix 1 0 1 0 0; # prueba flexien 

fix 1 1 1 0 0; # prueba cempresien 

# equalDOF $rNedeTag $nNedeTag $def1 (master, slave, def) 
equalDOF 12 10 2 3; # prueba memente (resiste:steel-i-cencrete) 
equalDOF 12 112 3; # prueba memente-rfaxial (resiste:cencrete) 

# nedal masses: 

# mass 8 $Mass 1 e-9 0.; # nede#, Mx My Mz, Mass=Weight/g, neglect retatienal inertia at nedes 

# define geemetric transfermatien: perferms a linear geemetric transfermatien ef beam stiffness and 

resisting ferce frem the basic system te the glebal-ceerdinate system 
set CelTransfTag 1 ; # asseciate a tag te celumn transfermatien 

geemTransf PDelta $CelTransfTag ; 

# element cennectivity: 

set numIntgrPts 2; # number ef 

integratien peints fer ferce-based element 
#element dispBeamCelumn 1 9 2 $numlntgrPts $SecTag1 $CelTransfTag; 
element nenlinearBeamCelumn 1 1 2 $numlntgrPts $SecTag1 $CelTransfTag; 
element nenlinearBeamCelumn 2 2 3 $numlntgrPts $SecTag2 $CelTransfTag; 
element nenlinearBeamCelumn 3 3 4 $numlntgrPts $SecTag3 $CelTransfTag; 
element nenlinearBeamCelumn 4 4 5 $numlntgrPts $SecTag4 $CelTransfTag; 
element nenlinearBeamCelumn 5 5 6 $numlntgrPts $SecTag5 $CelTransfTag; 
element nenlinearBeamCelumn 6 6 7 $numlntgrPts $SecTag5 $CelTransfTag; 
element nenlinearBeamCelumn 7 7 8 $numlntgrPts $SecTag6 $CelTransfTag; 

#element zereLengthSectien 8 9 1 $SecTag7 -erient 0 10-100; #vect eje leng glebceerd, vect en 
seccien el. glebceerd 

#element zereLengthSectien 9 9 1 $SecTag8 -erient 0 10-100; 
element zereLengthSectien 8 9 10 $SecTag7 -erient 0 10-100;# prueba memente 
(resisteisteel-rcencrete) 
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element zeroLengthSection 9 9 11 $SecTag8 -orient 0 10-100;# prueba memento-i-faxial 
(resiste:concrete) 

element nenlinearBeamColumn 10 12 1 $numlntgrPts $SecTag9 $ColTransfTag; 

# Define RECORDERS 

set dataDir Data; # set up name of data directory - simple 

file mkdir $dataDir; # create data directery 

recorder Node -file Data/DFree.out -time -nede 6 -dof 1 disp; # displacements of free nodes 

recorder Node -file Data/RBaseM.out -time -node 9 -dof 3 reaction; # support reaction 

recorder Node -file Data/RBaseV.out -time -node 12 -dof 1 reaction; # support reaction 

#recorder Element -file Data/RBaseELE.out -time -ele 1 sectien 1 ferce; # support reaction 

recorder Element -file Data/SteeIRS.out -time -ele 8 sectien fiber 8.04 -1 .16 $IDbond9 stressStrain; 
recorder Element -file Data/CoverRS.out -time -ele 9 section fiber -$ro 0 $IDzero1 1 stressStrain; 
recorder Element -file Data/CoreRS.out -time -ele 9 sectien fiber -$rc 0 $IDzero10 stressStrain; 
recorder Element -file Data/SteellOB.out -time -ele 10 sectien 1 fiber 8.04 -1 .16 $IDreinf5 stressStrain; 
recorder Element -file Data/Coverl OB.out -time -ele 10 section 1 fiber -$re 0 $IDcencCover4 stressStrain; 
recorder Element -file Data/Corel OB.out -time -ele 10 sectien 1 fiber -$rc 0 $IDconcCore1 stressStrain; 
recorder Element -file Data/SteellOT.out -time -ele 10 section 2 fiber 8.04 -1 .16 $IDreinf5 stressStrain; 
recorder Element -file Data/CeverlOT.eut -time -ele 10 section 2 fiber -$ro 0 $IDconcCover4 stressStrain; 
recorder Element -file Data/CerelOT.out -time -ele 10 sectien 2 fiber -$rc 0 $IDcencCore1 stressStrain; 
recorder Element -file Data/Steell .out -time -ele 1 section 1 fiber 8.04 -1 .16 $IDreinf5 stressStrain; 
recorder Element -file Data/Ceverl .eut -time -ele 1 section 1 fiber -$ro 0 $IDcencCover4 stressStrain; 
recorder Element -file Data/Corel .out -time -ele 1 section 1 fiber -$rc 0 $IDconcCore1 stressStrain; 

# define GRAVITY 

pattern Plain 1 Linear { 

load 8 0 -$Weight 0 



# Gravity-analysis parameters - lead 
set Tol 1 .Oe-8; # 

constraints Plain; # 

numberer Plain; # 

system SparseGeneral -piv; 
analysis 

#test NormDispIncr $Tol 6 ; 

an iteration step 
test Energyincr 1.0e-9 10 
algerithm Newten; 

at every iteratien 

set NstepGravity 10; # 

set DGravity [expr 1 ./$NstepGravity]; 
integrator LoadControl $DGravity; # 
analysis Static; # 

analyze $NstepGravity; # 



-controlled static analysis 
convergence tolerance fer test 
how it handles boundary conditions 

renumber dot's to minimize band-width (optimization), if you want to 

# how to store and solve the system of equations in the 

# determine if convergence has been achieved at the end ef 



# use Newton's solution algorithm: updates tangent stiffness 

apply gravity in 1 0 steps 

# first lead increment; 

determine the next time step for an analysis 
define type of analysis static or transient 
apply gravity 



# maintain censtant gravity leads and reset time to zero 

loadConst -time 0.0 



puts "Medel Built' 
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Programa 2 

# 

# Build a section 

# Silvia Mazzoni & Frank McKenna, 2006 

# Modified by Pablo Caiza 1 1 -2008 

# 



# SET UP 



wipe; # clear memory of all past model definitions 

model BasicBuilder -ndm 2 -ndf 3; # Define the model builder, ndm=#dimension, ndf=#dofs 
set dataDir Data; # set up name of data directory - simple 

file mkdir $dataDir; # create data directory 

source LibUnits.tcl; # define units 



# MATERIAL parameters 

# material ID tag 



set IDconcCorel 1 ; 


set IDconcCore2 2; 


set IDconcCore3 3; 


set IDconcCover4 4; 


set IDreinf5 


5; 


set IDreinf6 


6; 


set IDreinf7 


7; 


set IDreinf8 


8; 


set IDbond9 


9; 


set IDzerolO 


10; 


set IDzerol 1 


11; 


set IDcfrp12 


12; 


set IDcfrp13 


13; 



# Concrete for cover reinforced with CFRP 
#cfrp12 

uniaxialMaterial Hysteretic$IDcfrp1 2 1.145 0.0002 0.5725 0.0006 0.0 0.001 -11.45-0.002 -12.51 - 
0.0215 0 -0.0248 0.0 0.0 0.0 0.0 0.0 

#Concrete for core reinforced with CFRP 

# cfrpIS 

uniaxialMaterial Hysteretic $IDcfrp1 3 1.613 0.0005 0.8065 0.0015 0.0 0.0025 -16.13-0.005 -17.0- 
0.0215 -3.23-0.035 0.0 0.0 0.0 0.0 0.0 



# nominal concrete compressive strength 

setfc [expr -1 1*$ksi]; # ksi (-rTension, -Compression) 

set Ec [expr 57*$ksi*sqrt(-$fc/$psi)]; # Concrete Elastic Modulus 

#concCore1 (pitch lin) OJO! VALORES MODIFICADOS 

setfcll -16.13; # -12.28 ksi, maximum stress 

setepsll -0.00501; # -0.01 15 strain at maximum stress 

set fc21 -3.23; # -7.69 -2.46 ultimate stress 

set eps21 -0.044; # -0.029 strain at ultimate stress 



# concCore2 
set fc12 
set eps12 
set fc22 
set eps22 



(pitch 1 .5in) 
-14.55; 

-0.00437; 

-2.91; 

-0.0224; 



#-10.99 ksi, maximum stress 

# -0.00900 strain at maximum stress 
# -6.06 -2.20 ultimate stress 

# -0.0291 strain at ultimate stress 



# concCore3 
set fc13 
set eps13 
set fc23 
set eps23 



(pitch 3in) 

-12.80; 

-0.00364 ; 
-2.56; 

-0.0141; 



# -9.53 ksi, maximum stress 

# -0.0061 strain at maximum stress 
# -4.26 -1 .91 ultimate stress 

# -0.0148 strain at ultimate stress 



# concCover4 (unconfined) 
setfc14 $fc; 



# UNCONFINED concrete, maximum stress 
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set eps14 


-0.0028; 


concrete 




set fc24 


0.0 


set eps24 


-0.01; 


set lambda 0.1 ; 





$Ec 

# tensile-strength properties 

set ft14 [expr 0.0075*sqrt(-$fc*1000)]; 

set epst14 -0.01 ; 

set ft24 [expr0.414*$ft14]; 

set beta 4; 

set epst24 0.02; 

set Ets4 70; 



# 0.0031 strain at maximum strength of unconfined 

# [expr 0.1*$fc14] ultimate stress 

# 0.01 strain at ultimate stress 

# ratio between unloading slope at $eps24 and initial slope 



# tensile strength -rtension 



# reinfS 

set Fy5 [expr 71 .0*$ksi]; # STEEL yield 

set Es5 [expr 29000. *$ksi[; 

set Bs5 0.005; 

set R05 18; 

set cR15 0.925; 

set cR25 0.15; 



stress 

# modulus of steel 

# strain-hardening ratio 

# control the transition from elastic to 

# control the transition from elastic to 

# control the transition from elastic to 



plastic branches 
plastic branches 
plastic branches 



# re inf 6 

set Fy6 [expr 243.5*$ksi];# STEEL yield stress 230 248.8 

set Es6 [expr 27890. *$ksi[; # modulus of steel 27000 

set Bs6 0.01174; # strain-hardening ratio 0.008 0.00772 

set R06 7.344; # 7.39 control the transition from elastic to plastic branches 

set cRI 6 0.925; # 1 .0852 control the transition from elastic to plastic branches 

set cR26 0.1 5; # 0.00772 control the transition from elastic to plastic 

branches 

set inistress6 21 .5; # initial stress [ksi]21 .5 10.7 98 202.614 



# reinf7 

set Fy7 [expr 243.5*$ksi];# STEEL yield stress 100 235 
set Es7 [expr 27890. *$ksi]; # modulus of steel 

set Bs7 0.01174; # strain-hardening ratio 

set R07 7.344; # control the transition from elastic to plastic branches 

set cRI 7 0.925; # control the transition from elastic to plastic branches 

set cR27 0.15; # control the transition from elastic to plastic branches 

set inistress7 128.7; # initial stress [ksi[128.7 64.4 98 202.614 170 



# reinf8 

set Fy8 [expr 243.5*$ksi];# STEEL yield stress 100 230 
set Es8 [expr 27890. *$ksi]; # modulus of steel 

set Bs8 0.01174; # strain-hardening ratio 

set R08 7.344; # control the transition from elastic to plastic branches 

set cRI 8 0.925; # control the transition from elastic to plastic branches 

set cR28 0.15; # control the transition from elastic to plastic branches 

set inistress8 171.6; # initial stress [ksi[ 138 98 202.614 170 

uniaxialMaterial ConcreteOI $IDconcCore1 $fc11 $eps11 $fc21 $eps21 ; 
uniaxialMaterial ConcreteOI $IDconcCore2 $fc12 $eps12 $fc22 $eps22 ; 
uniaxialMaterial ConcreteOI $IDconcCore3 $fc13 $eps13 $fc23 $eps23 ; 
uniaxialMaterial ConcreteOI $IDconcCover4 $fc14 $eps14 $fc24 $eps24 ; 

#uniaxialMaterial Concrete02 $IDconcCore1 $fc11 $eps11 $fc21 $eps21 $lambda $ft14 $Ets4; 
#uniaxialMaterial Concrete02 $IDconcCover4 $fc14 $eps14 $fc24 $eps24 $lambda $ft14 $Ets4; 
#uniaxialMaterial Concrete03 $IDconcCore1 $fc11 $eps11 $fc21 $eps21 $lambda $ft1 4 $epst1 4 $ft24 
$beta $epst24 ; 

#uniaxialMaterial Concrete03 $IDconcCover4 $fc14 $eps14 $fc24 $eps24 $lambda $ft14 $epst14 $ft24 
$beta $epst24; 



uniaxialMaterial Steel02 $IDreinf5 $Fy5 $Es5 $Bs5 $R05 $cR15 $cR25; 
# uniaxialMaterial ElasticPP 5 27000 0.00245; 
uniaxialMaterial ElasticPP 6 27890 0.0008 -0.0004 -0.0008; 
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uniaxialMaterial Steel02 $IDreinf7 $Fy7 $Es7 $Bs7 $R07 $cR17 $cR27 0 1 0 1 $inistress7; 
uniaxialMaterial Steel02 $IDreinf8 $Fy8 $Es8 $Bs8 $R08 $cR18 $cR28 0 1 0 1 $inistress8; 

# Steel bond slip behavior 
#factor=2.5 Esec=145[ksi] 

#uniaxialMaterial Hysteretic $IDbond9 71 0.009 98 0.366 100 0.4 -71.0-0.029 -83-0.7945 -90.0- 

1.56 0.0 0.0 0.0 0.0 0.0 
#factor=1.25 Esec=250[ksi] 

#uniaxialMaterial Hysteretic $IDbond9 71 0.017 98 0.522 100 0.57 -71.0 -0.029 -83 -0.7945 -90.0- 

1.56 0.0 0.0 0.0 0.0 0.0 

# factor=1 .7 Esec=1 80[ksi] 

#uniaxialMaterial Hysteretic $IDbond9 71 0.013 98 0.593 100 0.6 -71.0-0.029 -83-0.7945 -90.0- 

1.56 0.0 0.0 0.0 0.0 0.0 

# actual values (august 2009) 

uniaxialMaterial Hysteretic $IDbcnd9 71 0.013 98 0.262 100 0.3 -71.0-0.029 -83 -0.7945 -90.0- 

1.56 0.0 0.0 0.0 0.0 0.0 

# actual values (finales de august 2009) 

#uniaxialMaterial Hysteretic $IDbcnd9 71 0.0024 98 0.262 -71.0-0.0024 -98-0.262 0.0 0.0 0.0 0.0 

0.0 

#Ccncrete fcr the bend slip zerc-length element 

#uniaxialMaterial CcncreteOI $IDzerc10 -12.28 -0.138 -2.46 -.3528; # cere: strains*12=displ. 
#uniaxialMaterial CcncreteOI $IDzerc11 -8 -0.037 0.0 -0.12; # ccver:strains*12=displ. 

#uniaxialMaterial CcncreteOI $IDzerc10 -11 -0.0196-7.54-0.07; # cere: ccncrete crushing using shear. 
#uniaxialMaterial CcncreteOI $IDzerc1 1 -11 -0.0196 -7.54-0.07; # ccver:ccncrete crushing using shear. 
uniaxialMaterial CcncreteOI $IDzerc10 -11 -0.0028-8.15-0.07; # cere: -7.54 ccncrete crushing using 
shear. 

uniaxialMaterial CcncreteOI $IDzerc1 1 -11 -0.0028-8.15-0.07; # ccver:-7. 54 ccncrete crushing using 
shear. 

#secticn Uniaxial 8 12 Mz; 

# secticn GEOMETRY 

set DSec 24; 

#data fcr reinf5 
set ccverSec5 [expr 4*$in]; 
set numBarsSec5 8; 
set barAreaSec5 1 .27; 

#data fcr reinf6 
set ccverSec6 [expr 3*$in]; 
set numBarsSec6 22; 
set barAreaSec6 0.153; 

#data fcr reinf7 
set ccverSec7 [expr 3*$in]; 
set numBarsSec7 22; 
set barAreaSec7 0.153; 

#data fcr reinf8 

set ccverSec8 [expr 3*$in]; # Cclumn cever tc prestressed steel NA. 

set numBarsSec8 22; # number cf unifcrmly-distributed Icngitudinal strands 

set barAreaSec8 0.1 53; # area cf Icngitudinal-reinfcrcement strands 

# Generate a circular reinfcrced ccncrete secticn 

# with cne layer cf steel evenly distributed arcund the perimeter and a cenfined cere. 

# cenfined cere. 

# by: Michael H. Scett, 2003 

# Ncte: 

# Zerc degrees is aicng secticn y-axis 

# 

set ri 0.0; # inner radius cf the secticn, cniy fcr hcllcw sectiens 

set rc [expr $DSec/2]; # cverall (cuter) radius cf the secticn 

set nfCcreR 19; # number cf radial divisiens in the cere (number cf "rings") 



# Cclumn Diameter [in] 

# Cclumn cever tc reinfcrcing steel NA. 

# number cf unifcrmly-distributed Icngitudinal bars 

# area cf Icngitudinal-reinfcrcement bars 



# Cclumn cever tc prestressed steel NA. 

# number cf unifcrmly-distributed Icngitudinal strands 
# area cf Icngitudinal-reinfcrcement strands 



# Cclumn cever tc prestressed steel NA. 

# number cf unifcrmly-distributed Icngitudinal strands 
# area cf Icngitudinal-reinfcrcement strands 
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set nfCoreT 22; 
set nfCoverR 6; 
set nfCoverT 22; 



# number of theta divisions in the core (number of "wedges" 
# number of radiai divisions in the cover 
# number of theta divisions in the cover 



# Define the fiber section 1 (1 pitch + normai reinf) 
set SecTagI 1 ; # iD 

section fiberSec $SecTag1 { 

set rc [expr $ro-$coverSec6]; # radius press. (Core 

radius) 

patch circ $IDconcCore1 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; # Define the core 

patch 

patch circ $IDconcCover4 $nfCoverT $nfCoverR 0 0 $rc $ro 0 360; # Define the covei 

patch 

set theta [expr 360.0/$numBarsSec6]; # angie increment betwee 

bars 

iayer circ $IDreinf6 $numBarsSec6 $barAreaSec6 0 0 $rc $theta 360; # Define the prestress 
iayer 



# Define the core 



# Define the cover 



# angie increment between 



fiber -8.04 1.16 $barAreaSec5 $IDreinf5 
fiber -6.14 5.32 $barAreaSec5 $IDreinf5 
fiber 0.00 8.13 $barAreaSec5 $IDreinf5; 
fiber 6.14 5.32 $barAreaSec5 $IDreinf5; 
fiber 8.04 -1 .16 $barAreaSec5 $IDreinf5 
fiber 6.14 -5.32 $barAreaSec5 $IDreinf5 
fiber -2.29 -7.80 $barAreaSec5 $IDreinf5 
fiber -6.14 -5.32 $barAreaSec5 $IDreinf5 



# south + 

# south + 

# south + 

# south + 

# south + 

# south + 

# south + 

# south + 



# Define the fiber section 2 (1 .5 pitch + normai reinf) 
setSecTag2 2; # iD 

section fiberSec $SecTag2 { 

set rc [expr $ro-$coverSec7]; # radius press. (Core 

radius) 

patch circ $IDconcCore2 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; # Define the core 

patch 

patch circ $IDconcCover4 $nfCoverT $nfCoverR 0 0 $rc $ro 0 360; # Define the covei 

patch 

set theta [expr 360.0/$numBarsSec7]; # angie increment betwee 

bars 

iayer circ $IDreinf7 $numBarsSec7 $barAreaSec7 0 0 $rc $theta 360; # Define the prestress 
iayer 



# Define the core 



# Define the cover 



# angie increment between 



fiber -8.04 1.16 $barAreaSec5 $IDreinf5 
fiber -6.14 5.32 $barAreaSec5 $IDreinf5 
fiber 0.00 8.13 $barAreaSec5 $IDreinf5; 
fiber 6.14 5.32 $barAreaSec5 $IDreinf5; 
fiber 8.04 -1 .16 $barAreaSec5 $IDreinf5 
fiber 6.14 -5.32 $barAreaSec5 $IDreinf5 
fiber -2.29 -7.80 $barAreaSec5 $IDreinf5 
fiber -6.14 -5.32 $barAreaSec5 $IDreinf5 



# south -I- 

# south -I- 

# south + 

# south + 

# south -I- 

# south -I- 

# south -I- 

# south -I- 



# Define the fiber section 3 (3 pitch + normai reinf) 
setSecTag3 3; # iD 

section fiberSec $SecTag3 [ 

set rc [expr $ro-$coverSec8j; # radius press. (Core 

radius) 

patch circ $IDconcCore3 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; # Define the core 

patch 

patch circ $IDconcCover4 $nfCoverT $nfCoverR 0 0 $rc $ro 0 360; # Define the covei 

patch 

set theta [expr 360.0/$numBarsSec7]; # angie increment betwee 

bars 

iayer circ $IDreinf8 $numBarsSec8 $barAreaSec8 0 0 $rc $theta 360; # Define the prestress 
iayer 



# Define the core 



# Define the cover 



# angie increment between 




118 



Revista CIENCIA Vol. 12 N.- 2 de 2009 



# Define the fiber section 4 (3 pitch oniy) 

setSecTag4 4; # iD 

section fiberSec $SecTag4 { 

set rc [expr $ro-$coverSec8]; # radius press. (Core 

radius) 

patch circ $IDconcCore3 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; # Define the core 

patch 

patch circ $IDconcCover4 $nfCoverT $nfCoverR 0 0 $rc $ro 0 360; # Define the cover 

patch 

set theta [expr 360.0/$numBarsSec7]; # angie increment between 

bars 

iayer circ $IDreinf8 $numBarsSec8 $barAreaSec8 0 0 $rc $theta 360; # Define the prestress 
iayer 



# Define the fiber section 5 (1 .5 pitch oniy) 

setSecTagSS; # iD 

section fiberSec $SecTag5 { 

set rc [expr $ro-$coverSec7]; # radius press. (Core 

radius) 

patch circ $IDconcCore2 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; # Define the core 

patch 

patch circ $IDconcCover4 $nfCoverT $nfCoverR 0 0 $rc $ro 0 360; # Define the cover 

patch 

set theta [expr 360.0/$numBarsSec7]; # angie increment between 

bars 

iayer circ $IDreinf7 $numBarsSec7 $barAreaSec7 0 0 $rc $theta 360; # Define the prestress 
iayer 



# Define the fiber section 6 (1 pitch oniy) 

setSecTag6 6; # iD 

section fiberSec $SecTag6 [ 

set rc [expr $ro-$coverSec6j; # radius press. (Core 

radius) 

patch circ $IDconcCore1 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; # Define the core 

patch 

patch circ $IDconcCover4 $nfCoverT $nfCoverR 0 0 $rc $ro 0 360; # Define the cover 

patch 

set theta [expr 360.0/$numBarsSec7]; # angie increment between 

bars 

iayer circ $IDreinf6 $numBarsSec6 $barAreaSec6 0 0 $rc $theta 360; # Define the prestress 
iayer 



# Define the fiber section 7 (Bond Siip-Fiexion) 

setSecTag7 7; # iD 

section fiberSec $SecTag7 [ 

# set rc [expr $ro-$coverSec6j; 
radius) 

# patch circ $IDzero1 0 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; 

# patch circ $IDzero1 1 $nfCoverT $nfCoverR 0 0 $rc $ro 0 360; 



fiber -8.04 1.16 $barAreaSec5 $IDbond9; # 

fiber -6.14 5.32 $barAreaSec5 $IDbond9; # 

fiber 0.00 8.13 $barAreaSec5 $IDbond9; # 

fiber 6.14 5.32 $barAreaSec5 $IDbond9; # 

fiber 8.04 -1.16 $barAreaSec5 $IDbond9; # 

fiber 6.14 -5.32 $barAreaSec5 $IDbond9; # 

fiber -2.29 -7.80 $barAreaSec5 $IDbond9; # 

fiber -6.14 -5.32 $barAreaSec5 $IDbond9; # 



# radius press. (Core 

# Define the core patch 

# Define the cover patch 
south + 

south + 
south -I- 
south -I- 
south + 
south + 
south -I- 
south -I- 



} 



# Define the fiber section 8 (Bond Siip-Compression) 
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setSecTag8 8; # ID 

section fiberSec $SecTag8 { 

set rc [expr $ro-$coverSec6]; 
radius) 

patch circ $IDzero10 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; 
patch circ $IDzero1 1 $nfCoverT $nfCoverR 0 0 $rc $ro 0 360; 

} 

# Define the fiber section 9 (Stiffened End) 
set SecTag9 9; 
set ro9 [expr $ro+1]; 

set nfCoverR9 8; # ID 

section fiberSec $SecTag9 { 

set rc [expr $ro-$coverSec6j; # radius press. (Core 

radius) 

# patch circ $IDconcCore1 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; # Define the core 

patch 

# patch circ $IDconcCover4 $nfCoverT $nfCoverR9 0 0 $rc $ro9 0 360; # Define the core 

patch 

patch circ $IDcfrp13 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; # Define the core patch 

patch circ $IDcfrp12 $nfCoverT $nfCoverR9 0 0 $rc $ro9 0 360; # Define the core 

patch 

# set theta [expr 360. 0/$numBarsSec6]; # angle increment between 

bars 

# layer circ $IDreinf6 $numBarsSec6 $barAreaSec6 0 0 $rc $theta 360; # Define the prestress 
layer 

# set rc1 [expr $ro-$coverSec5j; 

# setthetaf [expr 360. 0/$numBarsSec5]; 

# layer circ $IDreinf5 $numBarsSec5 $barAreaSec5 0 0 $rc1 $theta1 360; # Define the 

reinforcing layer 

fiber -8.04 1.16 $barAreaSec5 $IDreinf5 
fiber -6.14 5.32 $barAreaSec5 $IDreinf5 
fiber 0.00 8.13 $barAreaSec5 $IDreinf5 
fiber 6.14 5.32 $barAreaSec5 $IDreinf5 
fiber 8.04 -1.16 $barAreaSec5 $IDreinf5 
fiber 6.14 -5.32 $barAreaSec5 $IDreinf5 
fiber -2.29 -7.80 $barAreaSec5 $IDreinf5 
fiber -6.14 -5.32 $barAreaSec5 $IDreinf5 

} 

# Define the fiber section 10 (CFRP: 1 pitch -i- normal reinf) 

set SecTaglO 10; # ID 

section fiberSec $SecTag1 0 [ 

set rc [expr $ro-$coverSec6j; # radius press. (Core 

radius) 

patch circ $IDcfrp13 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; # Define the core patch 

patch circ $IDcfrp12 $nfCoverT $nfCoverR 0 0 $rc $ro 0 360; # Define the cover patch 

set theta [expr 360.0/$numBarsSec6j; # angle increment between 

bars 

layer circ $IDreinf6 $numBarsSec6 $barAreaSec6 0 0 $rc $theta 360; # Define the prestress 
layer 

fiber -8.04 1.16 $barAreaSec5 $IDreinf5 
fiber -6.14 5.32 $barAreaSec5 $IDreinf5 
fiber 0.00 8.13 $barAreaSec5 $IDreinf5 
fiber 6.14 5.32 $barAreaSec5 $IDreinf5 
fiber 8.04 -1.16 $barAreaSec5 $IDreinf5 
fiber 6.14 -5.32 $barAreaSec5 $IDreinf5 
fiber -2.29 -7.80 $barAreaSec5 $IDreinf5 
fiber -6.14 -5.32 $barAreaSec5 $IDreinf5 

} 

# Define the fiber section 1 1 (CFRP: 1 .5 pitch -i- normal reinf) 

set SecTagI 111; # ID 

section fiberSec $SecTag1 1 [ 



# south -I- 

# south -I- 

# south -I- 

# south -I- 

# south -I- 

# south -I- 

# south -I- 



# south -I- 

# south -I- 

# south -I- 

# south -I- 

# south -I- 

# south -I- 

# south -I- 



# radius press. (Core 

# Define the core patch 

# Define the cover patch 
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# radius press. (Core 



set rc [expr $ro-$coverSec7]; 
radius) 

patch circ $IDcfrp1 3 $nfCoreT $nfCoreR 0 0 $ri $rc 0 360; 
patch circ $IDcfrp1 2 $nfCoverT $nfCoverR 0 0 $rc $ro 0 360; 
set theta [expr 360.0/$numBarsSec7]; 
bars 

iayer circ $IDreinf7 $numBarsSec7 $barAreaSec7 0 0 $rc $theta 360; # Define the prestress 
iayer 



# Define the core patch 

# Define the cover patch 

# angie increment between 



fiber -8.04 1.16 $barAreaSec5 $IDreinf5 
fiber -6.14 5.32 $barAreaSec5 $IDreinf5 
fiber 0.00 8.13 $barAreaSec5 $IDreinf5 
fiber 6.14 5.32 $barAreaSec5 $IDreinf5 
fiber 8.04 -1.16 $barAreaSec5 $IDreinf5 
fiber 6.14 -5.32 $barAreaSec5 $IDreinf5 
fiber -2.29 -7.80 $barAreaSec5 $IDreinf5 
fiber -6.14 -5.32 $barAreaSec5 $IDreinf5 



# south 

# south 

# south 

# south 

# south 

# south 

# south -I- 

# south -I- 



set Dmax [expr 0.005*1 03]; # maximum dispiacement of pushover, push to 5% drift, 

set Diner [expr 0.0005*103]; # dispiacement increment for pushover, you want this to be very smali, 

but not too smaii to siow down the anaiysis 

set Hioad [expr $Weight]; # define the iaterai ioad as a proportion of the weight so that the pseudo time 
equais the iaterai-ioad coefficient when using iinear ioad pattern 
pattern Piain 200 Linear {; # define ioad pattern - generaiized 
ioad 6 $Hioad 0.0 0.0 0.0 0.0 0.0 

} 
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Programa 3 

set constraintsType Plain; # default; 

constraints $constraintsType 

numberer Plain 

#system BandGeneral 

system SparseGeneral -piv; # Overkill, but may need the pivoting! 
set Tol 1 .e-9; # Convergence Test: tolerance 

set maxNumIter 1 0; # Convergence Test: maximum number of iterations that will be performed 

before "failure to converge" is returned 

set printFlag 0; # Convergence Test: flag used to print information on convergence (optional) # 1 : 

print information on each step; 
set TestType Energyincr; # Convergence-test type 
test $TestType $Tol $maxNumlter $printFlag; 

set algorithmType Newton 
algorithm $algorithmType; 

#cycle 1 

foreach Aux1 {4 4} { 

set Nsteps [expr $Aux1/2]; # $Aux1 number of pushover analysis steps 

foreach Dincrl {.00721 -.00721 -.008755 .008755} ( 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nsteps]; # this will return zero if no convergence problems were encountered 

if {$ok 1= 0} ( 

puts "Analysis 1 failsl" 

}; # end if ok 10 



} 

#cycle 2 

foreach Auxi {4 4} ( 

set Nsteps [expr $Aux1/2]; # $Aux1 number of pushover analysis steps 

foreach Dincrl {.009785 -.009785 -.01 133 .01 133} [ 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nstepsj; # this will return zero if no convergence problems were encountered 

if [$ok 1= 0} [ 

puts "Analysis 2 failsl" 

}; # end if ok 10 

} 

} 

#cycle 3 

foreach Auxi [4 4} [ 

set Nsteps [expr $Aux1/2j; # $Aux1 number of pushover analysis steps 

foreach Dincrl (.012875 -.012875 -.01545 .01545} [ 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nstepsj; # this will return zero if no convergence problems were encountered 

if [$ok 1= 0} [ 

puts "Analysis 3 failsl" 

}; # end if ok 10 

} 

} 

#cycle 4 

foreach Auxi [12 12} [ 

set Nsteps [expr $Aux1/2j; # $Aux1 number of pushover analysis steps 

foreach Dincrl {.01 1 502 -.01 1 502 -.012188 .01 21 88} [ 
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integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nsteps]; # this will return zero if no convergence problems were encountered 

if {$ok != 0} { 

puts "Analysis 4 fails!" 

}; #endifok!0 

} 

} 

#cycle 5 

foreach Auxi {60 60 60 60 60 60} { 

set Nsteps [expr $Aux1/2j; # $Aux1 number of pushover analysis steps 

foreach Dincrl {.010334 -.01 0334 -.01 091 8 .01 091 8} { 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nsteps]; # this will return zero if no convergence problems were encountered 

if [$ok != 0} { 

puts "Analysis 5 fails!" 

}; # end if ok 10 

1 

1 

#cycle 6 

foreach Auxi [100 100 100 100} { 

set Nsteps [expr $Aux1/2]; # $Aux1 number of pushover analysis steps 

foreach Dincrl (.009579 -.009579 -.010321 .010321} { 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nstepsj; # this will return zero if no convergence problems were encountered 

if [$ok != 0} { 

puts "Analysis 6 fails!" 

}; # end if ok 10 

} 

} 

#cycle 7 

foreach Auxi [128 128 128 128} { 

set Nsteps [expr $Aux1/2]; # $Aux1 number of pushover analysis steps 

foreach Dincrl (.010332 -.01 0332 -.01 1 1 37 .01 1 1 37} { 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nstepsj; # this will return zero if no convergence problems were encountered 

if ($ok != 0} { 

puts "Analysis 7 fails!" 

}; # end if ok 10 

} 

} 

#cycle 8 

foreach Auxi (180 180} { 

set Nsteps [expr $Aux1/2]; # $Aux1 number of pushover analysis steps 

foreach Dincrl {.010346 -.01 0346 -.01 1 273 .01 1 273} { 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nstepsj; # this will return zero if no convergence problems were encountered 

if ($ok != 0} { 

puts "Analysis 8 fails!" 

}; # end if ok 10 

} 

} 

#cycle 9 

foreach Auxi (228 228} { 

set Nsteps [expr $Aux1/2]; # $Aux1 number of pushover analysis steps 
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foreach Dincrl {.010372 -.01 0372 -.01 1 1 94 .01 1 1 94} { 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nsteps]; # this will return zero if no convergence problems were encountered 

if [$ok != 0} { 

puts "Analysis 9 fails!" 

}; #endifok!0 

} 

} 

#cycle 10 

foreach Auxi [272 272} { 

set Nsteps [expr $Aux1/2j; # $Aux1 number of pushover analysis steps 

foreach Dincrl (0.010429 -0.010429 -0.010997 0.010997} ( 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nstepsj; # this will return zero if no convergence problems were encountered 

if [$ok 1= 0} [ 

puts "Analysis 10 fails!" 

}; #endifok!0 

} 

} 

#cycle 1 1 

foreach Auxi [488 488} [ 

set Nsteps [expr $Aux1/2j; # $Aux1 number of pushover analysis steps 

foreach Dincrl {.010338 -.010338 -.010389 .010389} [ 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nstepsj; # this will return zero if no convergence problems were encountered 

if [$ok != 0} [ 

puts "Analysis 1 1 fails!" 

}; #endifok!0 

} 

} 

#cycle 12 

foreach Auxi [400 400} [ 

set Nsteps [expr $Aux1/2j; # $Aux1 number of pushover analysis steps 

foreach Dincrl {0.020487 -0.020487 -0.020023 0.020023} [ 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nstepsj; # this will return zero if no convergence problems were encountered 

if [$ok != 0} [ 

puts "Analysis 12 fails!" 

}; #endifok!0 

} 

} 

#cycle 13 

foreach Auxi [560 560} [ 

set Nsteps [expr $Aux1/2j; # $Aux1 number of pushover analysis steps 

foreach Dincrl {.020589 -.020589 -.019978 .019978} [ 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nstepsj; # this will return zero if no convergence problems were encountered 

if [$ok != 0} [ 

puts "Analysis 13 fails!" 

}; #endifok!0 

} 

} 

#cycle 14 

foreach Auxi [720 720} [ 
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set Nsteps [expr $Aux1/2]; # $Aux1 number of pushover analysis steps 

foreach Dincrl {.02046 -.02046 -.019713 .019713} { 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nsteps]; # this will return zero if no convergence problems were encountered 

if [$ok != 0} { 

puts "Analysis 14 fails!" 

}; # end if ok 10 

1 

1 

#cycle 15 

foreach Auxi [840 840} { 

set Nsteps [expr $Aux1/2]; # $Aux1 number of pushover analysis steps 

foreach Dincrl {.020551 -.020551 -.019725 .019725} [ 
integrator DisplacementControl 6 1 $Dincr1 
analysis Static 

set ok [analyze $Nstepsj; # this will return zero if no convergence problems were encountered 

if [$ok != 0} [ 

puts "Analysis 15 fails!" 

}; #endifok!0 

} 

} 

puts "Pushover Done. Control Disp=[nodeDisp 6 1j" 




